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ABSTRACT 


Recent  studies  of  Cr  in  GaAs  have  established  that 
significant  Cr  diffusion  may  occur  during  crystal  growth  or 
subsequent  processing;  however,  most  work,  has  been  done  without 
concern  for  phase  equilibria  and,  as  a  result,  inconsistent  and 
invalid  results  are  common  in  the  literature.  Some  invest¬ 
igators  were  aware  of  the  need  for  phase  equilibria  information 
and  stated  that  their  work  was  limited  because  the  Ga-As-Cr 
ternary  system  had  not  been  determined. 

This  research  program  concerns  the  study  of  phase 
equilibria  in  the  Ga-As-Cr  system  and  applying  this  knowledge 
to  the  design  and  implementation  of  well-defined  experiments  to 
study  the  solubility  and  diffusion  of  Cr  in  GaAs. 

The  Ga-As-Cr  phase  diagram  (for  arsenic  a/o  £  50)  was 
determined  f "  m  room  temperature  to  1300°C  by  differential 
thermal  analysis  and  x-ray  diffraction.  Of  particular 
significance  to  solubility  and  diffusion  studies  are:  1)  Cr 
and  GaAs  are  not  in  equilibrium  with  each  other  under  any 
conditions;  2)  two  and  three  condensed  phase  regions  are 
present  which  would  be  suitable  Cr  sources  for  solubility  and 
diffusion  studies;  and  3)  a  Ga-rich  liquid  is  in  equilibrium 
with  GaAs  which  has  a  very  high  (up  to  10  a/o)  Cr  solubility. 

Well-defined  solubility  studies  of  Cr  in  GaAs  have  been 
performed  resulting  in  solubility  values  for  various  regions  in 
the  ternary  phase  diagram  and  distribution  coefficients.  The 
Cr  solubility  in  GaAs  for  the  GaAs ( s )-CrAs ( s )-Ga ( 4 ) 
region  ranges  from  1.3  x  10-*-®  to  1.2  x  10^  cm“3  from  700 


to  1000  JC ;  the  solubility  in  the  GaAs ( s ) -Ga ( z )  region  is 
lower.  Results  from  these  studies  confirm  the  applicability  of 
ternary  phase  equilibria  to  these  solubility  studies  and  also 
confirm  the  phase  diagram  determinations.  The  results  of  Cr 
solubility  studies  in  doped  GaAs  confirm  the  model  of  Cr  as  a 
deep  acceptor  on  a  Ga  site,  and  capacitance-voltage 
measurements  confirm  that  the  bulk  Cr,  measured  by  secondary 
ion  mass  spectrometry  (SIMS)  is  in  solution  and  electrically 
active. 

Three  types  of  diffusion  experiments,  indiffusion, 
outdif fusion,  and  "junction  diffusion”,  have  been  performed, 
again  utilizing  the  phase  equilibria  information.  The  Cr 
concentration-depth  profiles,  determined  by  SIMS,  have  been 
fitted  to  the  "substitutional/interstitial  dissociative 
diffusion"  model.  Three  regimes  of  the  Cr  diffusion  have  been 
identified  -  vacancy  diffusion  controlled,  vacancy  regeneration 
controlled,  and  interstitial  diffusion  controlled  -  and  values 
for  several  key  parameters  have  been  determined.  Diffusion 
studies  were  also  done  in  doped  GaAs  and  in  GaAs  from  different 
suppliers,  and  the  different  diffusion  results  can  be  explained 
in  terms  of  the  diffusion  model  and  the  defect  chemistry  of 
GaAs . 

In  addition,  some  of  the  anomalous  behavior  of  Cr  in  GaAs 
observed  in  this  work  and  in  other  studies  can  be  explained  in 
terms  of  the  phase  equilibria  information  and  the  solubility 
and  diffusion  results:  very  fast  Cr  diffusion  observed  in  some 
cases  is  accomplished  via  fast  interstitials;  large  Cr  build-up 
seen  at  the  GaAs  surface  after  a  diffusion  anneal  is  mostly  a 


result  of  Ga ( ^ j  formation;  and  a  depletion  of  Cr  near  the 
surface  after  Cr  indiffusion  is  due  to  the  surface  accumulation 


of  manganese. 

Finally,  the  results  from  this  research  have  been  used  in 
correlating  and  explaining  many  of  the  confusing  results  in  the 
literature  concerning  Cr  in  GaAs.  From  this  work,  a  much 
clearer  picture  of  Cr  in  GaAs  has  emerged. 
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1.  INTRODUCTION 


The  use  of  chromium  in  gallium  arsenide  in  microelectronic 
devices  is  widespread.  Chromium  is  intentionally  introduced 
into  the  GaAs  substrate  material  during  the  initial  growth  of 
the  bulk  single  crystal  in  order  to  compensate  the  electrically 
active  residual  defects  and  impurities,  and  render  the 
substrate  semi-insulating.  However,  there  have  been  problems 
with  the  incorporation  of  chromium  in  the  active  region  of  the 
devices,  either  during  the  epitaxial  growth  on  the  substrates, 
or  during  the  subsequent  process  annealing  of  the  wafers.  This 
has  often  resulted  in  poor  device  performance  and  has  prompted 
some  device  manufacturers  to  explore  other  possibilities  for 
semi-insulating  substrates,  particularly,  the  use  of 
ultra-pure  "undoped  semi- insulating  GaAs,"  grown  by  the  Liquid 
Encapsulated  Czochralski  (LEC)  method  (1-3).  However,  the 
majority  (60  -  80%)  of  GaAs  devices  are  still  being  made  using 
Cr-doped  GaAs  substrates  (4). 

Because  of  this  problem  with  Cr-doped  GaAs,  many  studies 
have  been  done  recently  on  the  redistribution  of  chromium  in 
gallium  arsenide.  In  many  of  these  studies,  solubility  and 
diffusivity  values  have  been  cited  or  used  in  attempts  to 
explain  the  behavior.  However,  in  none  of  these  studies  had 
ternary  phase  equilibria  been  considered,  simply  because  there 
was  no  information  on  appropriate  regions  of  the  Ga-As-Cr  phase 
diagram.  Furthermore,  only  a  small  percentage  of  the 
researchers  even  mentioned  the  necessity  of  considering  phase 
equilibria  information  in  reporting  or  using  other's  solubility 
and  diffusivity  information,  or  in  performing  their 


experiments.  As  a  result,  many  diffusivity  and  solubility 
values  used  in  explaining " uheir  results  are  not  valid. 

The  goals  in  this  work  have  been  the  following: 

1.  Experimentally  determine  the  regions  of  the  Ga-As-Cr 
ternary  phase  diagram  that  are  appropriate  for  solubility  and 
diffusion  studies  of  Cr  in  GaAs; 

2.  Use  this  knowledge  to  perform  valid  and  well-defined 
solubility  and  diffusion  experiments  of  Cr  in  GaAs; 

3.  Explain  other  phenomena,  such  as  surface  degradation 
and  the  anomalous  surface  build-up  of  Cr,  using  the  phase 
equilibria  information; 

4.  Correlate  and  explain  the  present  literature  on 
the  solubility  and  diffusivity  of  Cr  in  GaAs. 

The  outline  for  this  thesis  is  as  follows: 

In  chapter  2,  I  discuss  the  need  to  consider  ternary 
phase  equilibria,  and  to  know  the  ternary  phase  diagram  in 
conducting  solubility  and  diffusion  experiments  of  an  impurity 
in  a  compound  semiconductor  and  in  reporting  and  using  these 
values.  Also,  the  concept  of  "equilibrium’’  is  discussed  in 
relation  to  these  studies. 

In  chapter  3,  the  Ga-As-Cr  phase  diagram  determination 
is  presented,  along  with  its  implications  and  applications  to 
the  solubility  and  diffusion  studies. 

Chapter  4  is  devoted  to  the  solubility  studies  of  Cr  in 
GaAs.  Included  is  a  review  of  previous  work  concerning  the 
solubility  of  Cr  in  GaAs  and  how  earlier  results  were  in  error 
or  were  misleading  because  phase  equilibria  was  not 


considered.  The  results  of  the  present  work,  including 
solubility  values  of  Cr  in  undoped,  and  doped,  GaAs,  are 
presented,  as  well  as  results  which  confirm  both  the  validity 
of  applying  the  phase  equilibria  cons ideration  to  solubility 
studies,  and  of  the  phase  diagram  determinations.  In  addition, 
the  use  of  Secondary  Ion  Mass  Spectrometry  (SIMS)  is  discussed. 

Chapter  5  is  devoted  to  the  substitutional/interstitial 
dissociative  diffusion  mechanism,  including  a  literature 
review,  the  governing  equations  for  the  different  cases,  and 
its  application  to  real  systems. 

In  chapter  6,  I  discuss  my  work  on  the  diffusion  of  Cr  in 
GaAs.  Three  types  of  diffusion  experiments  -  indiffusion, 
outdif f us  ion ,  and  "junction"  diffusion  -  have  been  done,  all 
using  the  phase  diagram  and  solubility  results.  From  these 
studies,  a  previously  proposed  (5)  diffusion  mechanism  was 
confirmed,  expanded,  and  quantified.  In  addition,  what  has 
been  learned  about  GaAs  itself,  and  other  impurities  and 
defects,  from  these  studies  is  presented.  Previous  work 
concerning  the  diffusion  of  Cr  in  GaAs  is  also  discussed  and 
analyzed . 

In  chapter  7,  studies  of  various  surface  phenomena  in  GaAs 
are  reported.  These  include:  surface  build-up  of  Cr  after 
annealing;  surface  degradation  of  GaAs;  Cr  depletion  near  the 
surface,  even  after  indiffusion;  and  high  Cr  levels  in  LPE 
layers.  These  are  discussed  in  relation  to  the  previous 
findings. 

In  chapter  8,  a  summary  and  conclusions  of  this  work  will 
be  presented,  as  well  as  suggestions  for  future  work. 
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2.  THE  NEED  FOR  THE  TERNARY  PHASE  DIAGRAM 


2 . 1  Introduction 

The  study  of  a  dopant  in  a  compound  semiconductor  is  more 
complex  than  that  of  a  dopant  in  an  elemental  semiconductor 
because  of  the  extra  degree  of  freedom  (6).  The  number  of 
degrees  of  freedom  is  the  number  of  intensive  variables  that 
must  be  controlled  or  measured  in  order  for  an  equilibrium 
system  to  be  properly  specified  or  defined.  Once  the  correct 
number  of  degrees  of  freedom  in  a  specific  system  are 
specified,  then  the  system  is  "well-defined".  The  failure  to 
recognize  this  has  led  to  confusion,  and  to  erroneous  and 
invalid  conclusions,  in  the  past.  In  order  to  know  the  number 
of  degrees  of  freedom,  and  to  establish  if  the  system  is  at,  or 
even  near,  equilibrium,  the  ternary  phase  diagram  must  be 
known. 

The  importance  of  knowing  the  appropriate  ternary  phase 
diagram  in  diffusion  and  solubility  studies  of  an  impurity  in  a 
compound  semiconductor  has,  in  the  past,  only  occasionally  been 
espoused.  Examples  of  this  are:  Shih,  Allen,  and  Pearson's 
(6)  article  on  "Diffusion  in  Compound  Semiconductors  in  Terms 
of  the  Ternary  Phase  Diagram";  Tuck's  book  (7)  on  Diffusion  in 
Semiconductors ;  and  Casey's  review  article  (8)  on  "Diffusion 
in  the  III-V  Compounds  Semiconductors".  In  all  these  cases, 
the  example  of  Zn  in  GaAs  is  used,  mainly  because  it  is  the 
only  system  that  has  been  extensively  studied  in  this  manner, 
due  to  its  importance  in  GaAs  electronic  devices,  and  to  the 
fact  that  the  Ga-As-Zn  system  is  relatively  simple  and 
well-known  (9).  Casey  and  Panish's  article  (10)  on  Zn  in  GaAs 


is  another  excellent  example  of  ternary  phase  considerations  in 
diffusion  studies. 


A  good  review  of  ternary  phase  diagrams  in  general  is 
Reisman's  Phase  Equ i 1 ibr ia  (11).  (Additional  references  on 
determining  phase  diagrams  are  be  given  in  the  succeeding 
chapter . ) 

In  the  present  chapter,  the  basic  concepts  regarding  the 
relation  of  ternary  phase  equilibria  to  solubility  and 
diffusion  are  presented  first,  including  the  use  of  the  phase 
rule.  Then  the  need  for  the  ternary  phase  diagram  for 
establishing  well-defined,  and  for  citing  and  using  solubility 
and  diffusion  results,  is  discussed.  Finally,  the  need  for  the 
ternary  phase  diagram  for  equilibrium  conditions  will  be 
presented,  along  with  some  notes  regarding  equilibrium, 
quasi-equilibrium,  the  determination  of  the  final  equilibrium 
state,  and  suitable  sources  for  solubility  and  diffusion 
studies  based  on  the  above  considerations. 

2 . 2  Basic  Concepts 

The  number  of  degrees  of  freedom  needed  to  define  a  system 
is  stipulated  by  the  phase  rule: 

f  =  c  -  p  +  2  (2-1) 

where  f  is  the  number  of  degrees  of  freedom,  c  is  the  number  of 
components  (e.g.,  three  for  a  ternary  system,  two  for  a  binary 
system),  and  p  is  the  number  of  phases  present.  A  simple 
derivation  of  this,  taken  from  R.E.  Dickerson  (12),  follows. 

It  is  a  matter  of  experience  that  the  variables 
temperature,  pressure,  and  the  mole  fractions  of  each  component 
in  a  phase  fix  all  other  macroscopic  intensive  properties  of  a 


system,  including  mechanical,  thermal,  and  chemical 
properties.  This  means  that  for  p  phases,  2  +  p  c  intensive 
variables  are  sufficient  to  define  all  the  intensive  properties 
of  the  system  (assuming  temperature  and  pressure  are  constant 
throughout  the  system).  The  number  of  independent  intensive 
variables  is  the  total  number  of  variables  minus  the  number  of 
restraining  conditions.  There  are  two  types  of  restraining 
equations:  1)  the  sum  of  the  mole  fractions  in  each  phase  must 
equal  unity  -  there  are  p  such  equations;  and  2)  the  chemical 
potential  of  each  component  must  be  the  same  in  all  phases  - 
for  each  component,  there  are  p-1  equations  (  j  =  j 
=  ...  =  jP),  therefore  there  are  c(p-l)  equations  of  this 
type.  The  total  number  of  restraining  equations  is  thus 
p+p  c-c.  Therefore,  the  number  of  independent  variables  is: 
f  =  2+p  c  -  (p+p  c-c)  =  c+2-p,  as  earlier  stated. 

For  ternary  phase  diagrams,  with  3  components  by 
def inition : 

f  =  3+2 -p  =  5-p  (2-2  ) 

If  the  system's  total  composition  lies  in  a  three  phase  region 
(two  condensed  phases  plus  the  ever  present  vapor  phase),  then 
the  degrees  of  freedom  are  two.  Such  a  system  may  be  defined 
by  specifying  the  temperature  and  the  vapor  pressure.  Or 
alternately,  the  temperature  and  the  composition  of  one  of  the 
condensed  phases  may  be  specified.  For  a  four  phase  system 
(three  condensed  phases  plus  the  vapor  phase),  there  is  only 
one  degree  of  freedom.  In  this  case,  specifying  only  the 
temperature  completely  defines  the  system.  In  addition,  the 
overall  composition  of  the  system  (prescribed  amounts  of 


component  A,  B,  and  C)  lies  in  a  tie  triangle  region  in  a 
ternary  phase  diagram,  which  by  definition  is  a  region  in  which 
3  condensed  phases  are  present  and  enclosed  by  three  tie-lines 
(hence  a  tie-triangle).  As  long  as  the  system  composition  lies 
within  such  a  region  at  a  given  temperature,  the  system  is 
well-defined;  thus,  there  are  no  remaining  degrees  of  freedom. 

Once  a  system  is  well-defined,  all  intensive  variables  are 
fixed,  such  as  the  chemical  potentials  of  each  component  in 
each  phase,  the  partial  pressure  of  each  component,  the  vacancy 
concentration  in  each  phase,  and  such  related  properties  as  the 
solubility  and  diffusivity.  Even  if  the  net  composition  of 
the  system  is  changed,  as  long  as  the  composition  remains  in 
the  same  tie  triangle  region,  the  values  of  all  the  variables 
will  not  change.  Similarly,  if  the  system  does  change  from  one 
region  in  the  phase  diagram  to  another  (e.g.,  from  a 
3-condensed  phase  region  to  a  two  phase  region,  or  to  a 
different  3  phase  region),  or  from  one  position  in  a  two 
condensed  phase  region  to  another,  then  the  variables  mentioned 
above  are  expected  to  change. 

2 . 3  The  Need  for  The  Ternary  Phase  Diagram  for  Well-defined 

Conditions 

It  is  obvious,  then,  that  to  establish  well-defined 
conditions  for  studying  the  behavior  of  an  impurity  in  a 
compound  semiconductor,  one  needs  the  appropriate  ternary  phase 
diagram.  Without  this  knowledge,  one  would  not  know:  1)  if 
there  were  any  unspecified  degrees  of  freedom,  whether  it  is 
necessary  to  control,  for  example,  the  exact  composition  of  a 
liquid  phase  in  the  system;  2)  if  changing  such  quantities  as 
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the  amounts  of  the  components  or  the  volume  of  the  system 
influence  the  results;  or  3)  whether  the  results  of  a 
particular  experiment  could  be  compared  to  another  experiment. 
This  last  point  is  an  important  consideration  in  presenting 
results  and  using  results  of  others.  To  do  these  properly,  one 
must  have  knowledge  of  the  ternary  phase  diagram  since  the 
exact  conditions  of  the  experiment  may  alter  the  solubility, 
diffusivity,  and  partial  pressure,  and  usually  do.  The  person 
citing  results  (e.g.  the  "solubility  of  C  in  AB" )  should  report 
the  conditions  that  led  to  that  result,  and  if  possible,  the 
exact  location  in  the  phase  diagram.  Usually  the  experimental 
conditions  stipulate  at  least  the  first  part  of  this.  If  those 
results  are  then  used  to  explain  some  other  experimental 
result,  such  as  using  solubility  results  in  explaining 
diffusion  behavior,  one  must  make  sure  that  the  respective 
diffusion  and  solubilty  values  correspond  to  the  same  system, 
as  defined  by  the  phase  rule.  This  is  commonly  not  done,  as  I 
will  discuss  in  upcoming  chapters.  In  both  cases,  the  ternary 
phase  diagram  is  required. 

2 . 4  The  Need  for  the  Ternary  Phase  Diagram  for  Equilibrium 

Conditions 

Another  use  of  the  ternary  phase  diagram  in  these  studies 
is  to  establish  whether  the  system  is  at  or  near  equilibrium. 
The  best  explaination  of  this  is  by  example.  If  a  dopant  is 
not  in  equilibrium  with  the  semiconductor  material  for  the 
specified  conditions,  and  the  solubility  and  diffusion 
experiments  are  performed  by  using  the  pure  elemental  dopant  as 
the  source  material  (as  is  usually  done),  then  during  the 


coarse  of  the  experiment,  either  the  source,  or  the 
semiconductor  material,  or  both,  depending  on  the  position  in 
the  ternary  phase  diagram  indicated  by  the  overall  composition, 
will  undergo  a  change  in  phase.  In  that  case,  the  results  are 
not  valid.  This  is  especially  true  for  the  case  when  the 
entire  semiconductor  material  itself  is  changing  phase  and  the 
desired  result  is  the  solubility  of  a  specific  dopant  in  that 
semiconductor  material.  But  it  is  also  true  for  the  case  when 
the  source  material  is  changing  because:  1)  the  activity  of  the 
dopant  will  be  changing;  and  2)  in  changing  phase,  the  dopant 
usually  gains  components  from  the  semiconductor,  and  hence  part 
of  the  semiconductor,  usually  near  the  surface  and  thus 
critical  to  the  diffusion  and/or  solubility,  many  change 
phase.  To  know  if  the  initial  phases  are  at  equilibrium  and 
stable  with  respect  to  each  other  requires  a  knowledge  of  the 
phase  diagram.  As  will  be  discussed  later,  the  lack  of 
appropriate  phase  equilibria  was  the  major  limitation  in  the 
experiments  of  chromium  in  gallium  arsenide  performed  by  Tuck 
and  Adegboyega  (5),  as  was  indicated  by  Tuck  himself. 

2 . 5  Equilibrium  and  Quasi-equilibrium 

A  comment  should  be  made  regarding  equilibrium  and  the 
diffusion  process.  Diffusion  is  necessarily  a  non-equilibr : um 
process.  It  is,  in  fact,  a  process  in  which  equilibrium  is 
approached;  i.e.,  the  component  concentration  approaches  the 
solubility  level  in  the  host  material  for  the  given  source 
conditions.  However,  as  long  as  the  starting  materials  are 
different  from  the  equilibrium  phases  only  by  that  the  amounts 
of  the  components  ("dopants")  in  them  are  different,  one  may 


say  that  the  system  is  at  all  times  in  "quasi-equilibrium" 

(6).  In  such  cases,  the  chemical  potentials  of  the  diffusing 
species  at  the  surface  may  be  assumed  to  remain  constant 
throughout  the  experiment,  and  the  diffusion  experiment  is 
considered  valid.  This  is  not  the  case  when  the  starting 
materials  are  different  phases  from  the  equilibrium  phases  (as 
discussed  in  the  preceding  section).  In  that  case,  the 
chemical  potentials  may  undergo  a  significant  change. 

Another  comment  is  appropriate  regarding  the  equilibrium 
of  lattice  defects,  such  as  vacancies  and  interstitials, 
especially  deep  within  the  sample.  The  equilibrium 
concentrations  of  lattice  defects  may  take  a  relatively  long 
time  to  attain  equilibrium,  and  this  must  be  considered. 

Also,  equilibrium  and  phase  diagram  considerations  may 
affect  some  aspects  of  the  solubility  and  diffusion  behavior 
more  than  other  aspects.  For  example,  the  solubility  of  a 
component  in  a  phase  is  very  dependent  on  phase 
considerations,  as  has  been  discussed;  however,  the 
diffusivity  of  an  impurity  by  an  interstitial  mechanism  may 
not  be  as  dependent  on  the  component  chemical  potentials,  as 
dictated  by  other  phases  in  the  system,  and  thus  may  not  be  as 
dependent  on  the  phase  diagram  considerations. 

2 . 6  Determining  the  Final,  Equilibrium  State 

A  point  should  be  made  concerning  equilibrium  in  regards 
to  the  final,  equilibrium  state  of  the  system.  To  determine 
the  final  state  of  a  particular  experimental  system,  one 
usually  j  s t  determines  the  overall  composition  of  the  system 
(a/o  Ga,  As,  and  Cr,  for  example)  from  the  starting  material, 


finds  that  spot  on  the  isothermal  section  of  the  ternary  phase 
diagram  for  the  appropriate  temperature,  and  establishes  what 
phases  are  stable  for  that  composition.  However,  one  must 
consider  the  vapor  phase  since  conventional  isothermal  sections 
in  ternary  diagrams  only  consider  the  condensed  phases.  One 
must  account  for  the  loss  of  components  from  the  condensed 
phases  to  establish  the  equilibrium  vapor  pressure.  In  many 
cases  this  will  not  have  any  effect  on  the  final  state,  but  in 
some  cases  it  could.  An  example  of  a  case  where  it  will  not 
have  an  effect  is  in  an  indiffusion  performed  with  plenty  of 
the  source  material  and  a  large  diffusion  sample,  and  the 
overall  system  composition  lies  well  within  the  interior  of  a 
tie  triangle.  In  this  case,  losing  a  small  amount  of  one  or 
more  of  the  components  to  the  vapor  phase  will  probably  not 
shift  the  system  to  a  different  region  of  the  phase  diagram. 
(However,  it  may,  as  we  will  see  later,  affect  how  much  of  each 
phase  is  present. )  An  example  of  a  case  where  losing 
components  from  the  condensed  phase  to  establish  the 
equilibrium  vapor  pressure  may  have  a  large  effect  is  during  an 
outdiffusion  anneal  of  a  doped  semiconductor.  Here,  losing  the 
same  amount  of  a  component  as  in  the  above  case  could  very  well 
put  the  system  in  a  new  region  of  the  phase  diagram  just  due  to 
the  geometry  of  the  phase  diagram  and  the  relatively  small 
solubilities  of  dopants  in  semiconductors. 

In  all  the  considerations  mentioned  in  this  chapter,  it 


has  been  assumed  that  the  system  in  question  is  a  closed 
system.  Only  in  that  case  will  equilibrium  be  achieved.  It 
may  take  a  while,  and  it  may  require  that  some  of  the 
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components  leave  the  condensed  phases  and  go  into  the  vapor 
phase,  which  can  be  taken  into  account  without  much  trouble. 

But  if  the  system  is  open,  as  in  the  case  of  open  flow  tubes 
commonly  used  in  semiconductor  processing  and,  in  some  cases, 
in  studying  solubility  and  diffusion,  then  the  system  never 
reaches  a  final  equilibrium  state.  The  components  enter  the 
vapor  phase  to  try  to  achieve  equilibrium;  however,  they  are 
removed  before  equilibriun  is  established.  As  a  result,  the 
system  becomes  very  complex  to  study.  Therefore,  in  our 
studies,  a  closed  ampoule  is  used  in  all  cases.  (It  is 
possible,  though,  to  qualitatively  explain  some  phenomena 
reported  by  others  who  use  open  tubes. ) 

2 . 7  Suitable  Sources 

It  is  obvious  from  the  before-mentioned  considerations 
regarding  equilibrium  conditions  what  source  materials  should 
be  used  for  valid  solubility  and  diffusion  studies.  Because 
you  want  the  material  to  be  in  equilbrium  -  or  actually  "in 
quasi-equilibrium"-  with  the  phase  you  are  diffusing  the 
component  into,  the  source  material  should  be  made  up  of  phases 
that  all  lie  in  the  same  region  (e.g.,  a  tie-triangle  region) 
as  the  substrate  phase.  In  practice,  one  anneals  a  mixture 
which  has  an  overall  composition  in  the  region  so  that  the 
proper  equilibrium  phases  are  formed.  Then  the  actual 
solubility  or  diffusion  experiment  is  done  by  annealing  this 
source  with  the  undoped  substrate. 

2 . 8  Summary 

In  summary,  to  have  valid  solubility  and  diffusion 
experiments,  one  must  have  a  well-defined  equilibrium,  or 
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quasi-equilibrium,  system.  To  know  if  this  is  the  case,  one 
must  have  a  knowledge  of  the  appropriate  phase  diagram.  In  the 
case  of  an  impurity  in  a  compound  semiconductor,  this  means 
knowing  the  ternary  phase  diagram.  Further,  in  citing  results, 
or  using  other's  results,  one  must  be  aware  of  the  ternary 
phase  equilibrium  considerations. 


As  discussed  in  chapter  2,  knowledge  of  the  appropriate 
ternary  phase  equilibria  is  necessary  to  design  and  interpret 
solubility  and  diffusion  studies  of  an  impurity  in  a  compound 
semiconductor.  Since  the  Ga-As-Cr  ternary  phase  diagram  had 
not  previously  been  determined,  this  was  the  first  goal  of  my 
research. 

A  logical  starting  point  in  determining  a  ternary  phase 
diagram  is  the  knowledge  of  the  respective  binary  systems, 
including  a  complete  description  of  the  compounds,  with  their 
melting  points,  x-ray  diffraction  patterns  (for  identification 
purposes),  and,  if  possible,  the  Gibb’s  formation  energies  of 
the  compounds.  With  this  information,  one  can  predict  possible 
ternary  sections.  The  true  sections  may  then  be  established  by 
annealing,  quenching,  and  x-raying  appropriate  ternary 
compositions  at  appropriate  temperatures.  The  melting  points 
of  the  binary  compounds  establish  the  temperatures  where 
two  or  three  condensed  phase  regions  may  be  present,  and 
provide  clues  to  which  compounds  are  the  "dominant"  ones,  with 
tie  lines  radiating  from  them  to  other  compounds  and  elements; 
these  dominant  compounds  are  usually  high  melting  point 
compounds  which  are  centrally  located  along  the  binary  tie 
lines.  The  Ga-As-Cr  system  was  particularly  difficult  because 
there  are  nine  reported  intermediate  binary  compounds,  with 
incomplete  and  confusing  information  concerning  many  of  them. 

The  Ga-As  system,  with  only  one  intermediate  compound, 

GaAs ,  is  well  known  and  described  in  a  good  review  article 


(13),  and  the  Cr-Ga  system,  with  4  intermediate  compounds,  has 
been  extensively  studied  (14-16).  However,  the  Cr-As  system 
had  not  been  adequately  studied  prior  to  the  present  work,  and 
there  is  confusing  information  in  the  literature.  In  1960, 
Yuzuri  (17)  proposed  a  tentative  Cr-As  phase  diagram  for  the 
composition  region  above  50  a/o  Cr,  with  the  following  binary 
CrAs  compounds:  CrAs,  Cr3As2,  Cr2As,  and  Cr5As.  The  melting 
point  of  each  compound  was  only  estimated.  In  1967,  Boiler, 
Wolfsgruber,  and  Nowotny  (18)  reported  that  the  Cr5As  phase  was 
actually  Cr4As.  In  1970,  Baurecht,  Boiler,  and  Nowotny  (19) 
proposed  the  existence  of  Cr4As3;  however,  it  is  not  clear  from 
their  description  whether  this  compound  replaced  Cr3As2  or  is 
an  additional  compound.  (The  latter  statement  is  implied  in 
their  paper;  however,  the  fact  that  their  x-ray  pattern  for 
Cr^Asj  is  identical  to  Yuzuri 's  Cr3As2  pattern  implies  that  the 
former  statement  is  correct.)  All  references  mention  CrAs,  but 
there  are  uncertainties  regarding  its  x-ray  pattern.  Nowotny 
and  Arstad  (20)  published  an  x-ray  pattern  in  1937;  however, 
many  pairs  of  reflections  were  grouped  together  into  single 
d-spacings  and  29  values.  (From  calculations  using  their 
lattice  spacing  values,  I  have  found  that  these  pairs  split 
into  reflections  as  much  as  1°  apart.  This  led  to  many 
problems  in  identifying  the  compound. )  There  are  also 
inconsistencies  in  the  reported  melting  points  of  some  of  the 
compounds.  Yuzuri  (17)  estimated  the  melting  point  of  CrAs  to 
be  *800°C;  however  Baurecht,  Boiler,  and  Nowotny  (21),  in  their 
Cr-As-C  phase  diagram,  show  it  to  be  a  solid  at  1000°C,  thus,  a 
melting  point  of  greater  than  1000°C.  With  regard  to  the 
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compound  Cr4As3,  it  was  shown  to  be  present  in  the  Cr-As-C 
phase  diagram  at  1000°C  (21),  but  not  present  in  the  Cr-As-Ge 
phase  diagram  at  800°C  (18).  This  implies  a  peritectic 
decomposition  of  this  compound,  contrary  to  all  other  studies. 
The  compound  Cr4As  was  shown  to  be  present  at  800 °C  in  the 
Cr-As-Ge  system,  but  not  at  1000 °C  in  the  Cr-As-C  system.  This 
would  imply  a  melting  point  between  800  and  1000°C;  however, 
Yuzuri  estimated  its  melting  point  to  be  greater  than  1200°C. 
The  existence  of  Cr2As  and  its  x-ray  pattern  have  not  been 
disputed  (22).  However,  its  melting  point  has  only  been 
estimated,  and  low  and  high  temperature  forms  have  been 
reported,  with  no  transition  temperature  given.  In  addition,  a 
high  pressure  CrAs2  phase  on  the  arsenic  side  of  the  phase 
diagram  has  been  reported  (23). 

The  only  study  on  the  Ga-As-Cr  system,  besides  determining 
the  distribution  coefficient  of  Cr  in  GaAs(S)  versus  Cr  in 
Ga ( i )  (24-29),  was  the  determination  of  the  arsenic 
solubility  in  Cr-doped  Ga  melt,  described  by  Andre  and  LeDuc 
(30).  The  arsenic  and  chromium  concentrations  both  ranged  from 
2  to  10  a/o  over  a  temperature  range  of  600  to  900°C.  They  did 
not  say  whether  10  a/o  was  the  limiting  Cr  concentration  in  the 
Ga  melt,  only  that  it  was  at  least  that  high. 

My  initial  objective  was  to  confirm  the  information 
concerning  the  Cr-Ga  system,  and  to  synthesize  the  Cr-As 
compounds  and  experimentally  determine  necessary  information 
concerning  them,  such  as  their  melting  points  and  x-ray 
diffraction  patterns.  Then  using  the  annealing,  quenching,  and 
x-ray  method  and  differential  thermal  analysis,  I  could 


determine  the  important  features  of  the  Ga-As-Cr  ternary  phase 
diagram.  Since  I  was  mainly  interested  in  the  Cr  and  Ga  rich 
portion  of  the  system,  and  because  of  the  difficulty  in 
annealing  arsenic-rich  compositions,  I  limited  my  work  to  the 
arsenic-poor  portion  of  the  phase  diagram,  with  a/o  As  £  50. 

In  the  next  section,  the  experimental  techniques  and 
procedures  used  in  the  phase  equilibria  studies  are  discussed, 
followed  by  the  results  of  these  studies.  Then  the 
implications  and  applications  of  these  findings  to  solubility 
and  diffusion  studies  of  chromium  in  gallium  arsenide  are 
presented . 

3 . 2  Experimental  Techniques 

The  experimental  techniques  used  in  the  phase  equilibria 
studies  -  x-ray  diffraction,  differential  thermal  analysis, and 
sample  synthesis  and  annealing  -  will  be  discussed  in  the 
following  three  subsections,  in  which  the  basic  principles  of 
each  will  be  discussed  in  addition  to  the  particular 
experimental  methods  and  operations  employed  in  these  studies. 
Then  the  specific  experimental  procedures  concerning  the  Cr-Ga 
and  Cr-As  binary  systems,  and  the  Ga-As-Cr  ternary  system  will 
be  described. 

3.2.1  X-ray  Diffraction 

Unknown  compounds  may  be  identified  by  x-ray  diffraction 
( XRD ) .  The  basic  principle  of  this  method  is  the  behavior  of 
diffracted  x-ray  beams  from  a  crystal  in  accord  with  Bragg's 
law  Ubeam  =  2dj1]cj_sin  Q).  The  diffraction  angles,  a, 
depend  on  the  d-spacings,  and  hence  on  the  size  and  shape  of 
the  unit  cell  of  the  analyzed  sample.  In  addition,  the 


relative  intensity  of  each  diffracted  beam  depends  on  the 
contents  of  the  unit  cell.  If  a  powder  sample  is  analyzed,  the 
proper  Bragg  angle  will  be  achieved  for  all  possible 
crystallographic  planes.  The  sample  is  scanned  over  a  range  of 
0's,  and  the  intensity  of  the  diffracted  beam  at  each  0  is 
measured.  Because  of  the  above  considerations  regarding  unit 
cell  size,  shape,  and  contents,  each  element  or  compound  will 
have  a  characteristic  diffraction  pattern,  or  "finger  print". 
Therefore,  if  one  knows  the  x-ray  pattern  for  a  particular 
element  or  compound,  one  may  identify  that  material  in  an 
"unknown"  sample  by  matching  the  positions  of  the  peaks  and 
relative  intensities  of  the  obtained  diffraction  pattern  from 
the  "unknown"  sample  to  those  in  the  "known"  or  "standard" 
diffraction  pattern.  If  two  or  more  compounds  are  present  in 
the  powder  sample,  then  both  diffraction  patterns  will  appear 
in  the  obtained  pattern,  superimposed  on  each  other.  Hence,  an 
unknown  sample  composed  of  one,  two,  or  more  phases  may  be 
analyzed,  and  the  constituents  identified  by  obtaining  its 
diffraction  pattern  and  matching  the  peaks  to  standard  patterns 
of  the  different  phases.  (This  assumes,  of  course,  that  the 
pattern  of  each  compound  or  phase,  is  already  known.) 

In  this  work,  a  Norelco  x-ray  diffraction  apparatus,  using 
a  Cu  x-ray  source  (X  =  1.54178  A)  was  used.  The  beam  was 
operated  at  36  KeV  and  14  yamps.  The  samples  were  treated  as 
follows:  they  were  crushed  into  powders  with  a  ceramic  mortar 
and  pestle;  cemented  to  a  microscope  slide  using  a  mixture  of 
acetone  and  colloidion  (=>50-50);  and  slowly  heated  on  a  hot 
plate.  (I  discovered  that  this  mixture  dissolves  gallium,  so 
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in  the  cases  that  Ga  was  suspected  to  be  present  in  the  sample, 
i.e.,  whenever  the  composition  was  in  the  Ga  corner  or  side  of 
the  phase  diagram,  and  when  the  sample  was  not  easy  to  crush, 
indicating  a  liquid  phase  present,  I  would  also  do  an 
additional  XRD  analysis  using  Duco  cement  to  hold  the  powder  on 
the  slide).  The  usual  scan  rate  was  2°/minute,  and  in  all 
cases,  the  scan  was  done  from  2s  =  10°  to  85°. 

A  good  review  of  powder  x-ray  techniques  is  found  in 
Schwartz  and  Cohen's  Diffraction  from  Materials  (31). 

3.2.2  Differential  Thermal  Analysis 

Differential  Thermal  Analysis  (DTA)  was  used  to  determine 
phase  transition  temperatures  (melting  points  of  compounds, 
liquidus  and  solidus  temperatures  of  two  and  three  phase 
regions,  etc.)  for  various  binary  and  ternary  compositions. 

In  DTA  analysis,  the  sample  in  question,  along  with  a 
reference  material,  is  heated  slowly  to  a  maximum  temperature, 
and  then  cooled  slowly.  The  temperatures  of  the  sample  and  the 
reference  material  during  the  heating/cooling  cycle  are  each 
measured  by  separate  thermocouples.  The  output  is  the 
difference  between  those  two  readings  versus  the  system 
temperature  (measured  by  a  third  thermocouple),  or  sometimes 
just  the  temperature  of  the  reference  material.  As  long  as  no 
phase  transitions  are  occurring,  and  the  heat  capacities  of  the 
two  materials  (sample  reference)  are  similar,  the  output  should 
yield  a  straight  horizontal  line.  If  there  is  a  slope,  it  can 
be  electrically  adjusted  to  make  up  for  any  difference  in  heat 
capacities  or  masses.  But  when  a  phase  transition  occurs  in 
the  sample  (such  as  solid*l iquid ) ,  then  the  temperature  of  the 
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sample  relative  to  the  reference  changes.  For  example,  during 

melting,  the  temperature  of  the  sample  remains  constant  until  0 

the  entire  phase  turns  to  liquid,  while  the  temperature  of  the  ' 

reference  material  continues  to  increase.  The  result  in  the  : 

output  in  this  case  is  a  downward  pointing  spike,  and  the 

temperature  (x-axis)  at  which  this  begins  is  the  melting 

point.  Similarly,  upon  cooling,  an  upward  spike  or  peak  occurs 

upon  solidifying.  Because  of  superheating  and  supercooling, 

one  usually  does  both  heating  and  cooling  runs,  and  the  actual  • 

phase  transformation  temperature  is  in  between,  usually  much 

closer  to  the  heating  peaks  since  superheating  is  usually  much  , 

0 

less  than  supercooling.  Any  phase  transition  in  which  the  heat 
capacity  of  the  material  changes,  or  which  is  endothermic  or 
exothermic,  will  cause  a  change  in  slope  of  the  DTA  curve  or  a 

0 

peak,  thus  enabling  one  to  measure  the  phase  transformation 
temperatures,  and  hence  to  obtain  melting  points,  solidus 
lines,  and  liquidus  lines. 

For  this  work,  a  DuPont  1090  Thermal  Analyzer  (controller) 
with  a  Dupont  1600°C  DTA  cell  was  used.  This  uses  Pt-Pt  13  Rh  'j 

sample,  reference,  and  control  thermocouples.  The  samples  ] 

. —  .J 

(-.01  gm)  and  reference  material  were  sealed  in  evacuated 

(»10-6  torr )  quartz  DTA-style  ampoules,  which  were  cleaned 

prior  to  use  with  a  10:1  HN03/HF  solution  and  rinsed  with  D.I. 

water.  These  were  placed  in  55  mm.  Pt  crucibles,  which  fit  0 

over  the  thermocouples.  The  reference  sample  was  A1203,  which 

does  not  undergo  any  phase  transformation  over  the  temperature 

range  examined  (i.e.,  room  temp,  to  1400°C).  The  normal  , 

heating  rate  was  20°/min.;  10°/min.  was  used  occasionally,  with 
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identical  results.  The  cited  temperature  precision  for  the 


instrument  is  ±  2°C.  The  phase  transition  temperatures  and 
melting  points  determined  in  this  work  are  conservatively 
estimated  to  be  +  5°C.  The  system  was  calibrated  using  a 
silver  sample  (m.p.  of  960. 8°C).  Each  sample  was  heated  to  the 
maximum  temperature  (usually  1400°C)  and  cooled  to  300°C,  and 
the  cycle  was  repeated  twice  to  insure  complete  mixing  and 
reaction  of  the  components. 

To  ensure  that  the  samples  w^re  composed  of  the  correct 
compositions,  a  large  amount  (  =>2  gm)  of  each  sample  was  made 
from  appropriate  amounts  of  GaAs,  Ga,  Cr,  and  CrAs,  which  was 
annealed  in  a  large  quartz  ampoule  in  a  annealing  furnace  at 
=950°C  for  8  hours.  This  was  slow-cooled  to  room  temperature 
(it  was  found  that  if  it  was  fast-cooled ,  or  quenched,  liquid 
gallium  often  formed  which  made  the  sample  difficult  to  crush), 
then  crushed,  stirred,  and  =».01  gm  of  it  poured  into  the  DTA 
ampoule,  which  was  evacuated  and  sealed. 

3.2.3  Synthesis  and  Annealing  of  Samples 

In  this  section,  the  general  techniques  for  the  synthesis 
and  annealing  of  the  samples  are  given;  the  details  of  the 
synthesis  of  the  various  compounds  will  be  presented  in  the 
following  section.  For  both  synthesizing  the  various  binary 
compounds  and  annealing  the  different  mixtures  of  Cr,  Ga,  and 
As  to  determine  which  phases  are  stable  for  that  composition 
for  the  phase  diagram  determination,  the  following  procedure 
was  used:  the  material  was  weighed,  placed  in  quartz  ampoules, 
which  were  evacuated  to  =10"^  torr,  sealed,  and  placed  in  a 
furnace  for  the  appropriate  time  and  temperature.  The  furnace 


used  was  a  one  zone  tubular  Marshall  furnace,  with  a  maximum 
temperature  of  1200 °C.  A  sodium  heat  pipe  was  used  to  provide 
a  uniform  temperature  zone  (±1°C  along  the  12  inch  length),  and 
dampened  any  temperature  decrease  when  placing  a  cold  ampoule 
in  the  furnace.  The  temperature  was  controlled  by  a  Eurotherm 
919  controller,  with  a  45  amp  SCR.  A  chromel-alumel 
thermocouple  placed  right  next  to  the  ampoule  monitored  the 
temperature.  The  thermocouple  was  calibrated  by  performing 
cooling  curve  experiments  with  silver  (m.p.  =  960. 5°C)  and  zinc 
(m.p.  =  419. 5°C).  The  quartz  ampoules  used  were  constructed  by 
me  (made  from  Amersil  and  Spectrophysics  semiconductor  grade 
quartz  tubing),  and  were  cleaned  prior  to  use  with  a  10%  HF  in 
HNO 3  solution,  and  rinsed  in  deionized  water. 

3.2.4  The  Cr-Ga  System;  Experimental 

The  Cr-Ga  compounds,  Cr3Ga,  CrGa,  Cr4Ga3,  and  CrGa4,  were 
synthesized  and  analyzed  in  order  to  confirm  their  melting 
points  and  x-ray  patterns.  They  were  synthesized  by  annealing 
appropriate  amounts  of  Cr  and  Ga  (99.9999%  pure,  from 
Atomergic)  in  evacuated  quartz  ampoules  for  24  hours  at  950°C, 
crushing  and  mixing  the  contents,  and  resealing  and  reannealing 
for  24  more  hours  at  950 °C.  Each  was  then  cooled  to  40 °C  below 
its  reported  melting  point  (except  Cr3Ga,  since  its  melting 
point  is  too  high  for  the  annealing  system),  and  water  quenched 
at  that  temperature.  Each  was  x-rayed.  Three  intermediate 
compositions,  between  CrGa  and  Cr3Ga4,  between  Cr3Ga4  and 
CrGa4,  and  between  Cr3Ga  and  CrGa,  were  annealed,  x-rayed,  and 
analyzed  by  DTA  to  check  the  melting  points  and  transition 


temperatures 


The  Cr-As  compounds,  CrAs,  Cr4As3,  Cr2As,  and  Cr4As,  were 
synthesized  by  first  synthesizing  CrAs  and  then  adding 
appropriate  amounts  of  Cr  to  it  and  annealing  to  synthesize  the 
others.  This  minimizes  the  number  of  times  that  elemental 
arsenic  has  to  be  used.  The  synthesis  of  CrAs  was  no  easy 
task,  as  evidenced  by  the  fact  that  three  commercial  chemical 
companies  were  unsuccessful  at  synthesizing  it.  It  was  finally 
necessary  to  synthesize  pure  CrAs  myself.  This  procedure 
consisted  of  annealing  appropriate  amounts  of  powdered  chromium 
(99.99%  pure,  from  Atomergic)  and  arsenic  (99.9999%  pure,  from 
Atomergic)  in  a  thick-walled,  evacuated  and  sealed,  quartz 
ampoule.  (The  arsenic  was  first  heated  to  300°C  in  an 
evacuated  condenser  apparatus  which  removed  any  oxide  layer  on 
the  arsenic  by  depositing  the  oxide  at  the  cold  end  of  the 
tube.  The  pure  arsenic  was  then  crushed  and  added  to  the 
chromium,  with  no  excess  arsenic. )  The  temperature  in  the 
furnace  was  ramped  to  600 °C  at  a  rate  of  10 °C/  hour  and  left 
there  for  200  hours.  After  examining  the  contents  by  x-ray 
diffraction  methods,  and  crushing  and  resealing  the  sample,  it 
was  returned  to  the  furnace  at  600 °C  and  ramped  up  to  800°C  at 
10°/hr. ,  and  annealed  for  200  additional  hours.  No  evidence  of 
elemental  arsenic  was  observed,  as  evidenced  by  the  absence  of 
yellowish  or  greenish  vapor,  characteristic  of  arsenic.  The 
sample  was  slowly  cooled  to  room  temperat'  and  crushed  into  a 
fine  powder.  Subsequent  x-ray  analysis  produced  a  clean  x-ray 
pattern  of  CrAs.  The  melting  point  was  determined  by  DTA.  The 
other  Cr-As  compounds  were  then  synthesized  by  adding  Cr 


CrAs  and  annealing  in  a  similar  fashion.  Clean  x-ray  patterns 
were  obtained  and  are  consistent  with  previously  reported 
patterns.  Melting  points  were  determined  by  DTA.  In  addition, 
as  intermediate  composition  between  Cri,As3  and  Cr2As  was 
annealed,  x-ray'd,  and  analyzed  by  DTA. 

3.2.6  The  Ga-As-Cr  Ternary  System;  Experimental 

As  discussed  earlier,  once  the  respective  binary  systems 
are  known,  ternary  phase  diagrams  are  usually  determined  by 
annealing  various  three-component  compositions  at  various 
temperatures  and  quenched  to  establish  which  phases  are  stable 
for  that  composition  and  temperature.  In  this  manner, 
appropriate  tie  lines  may  be  drawn  connecting  stable  phases  and 
enclosing  tie-triangle  regions.  In  addition,  DTA  may  be  used 
to  determine  the  transition  temperatures  of  liquidus  and 
solidus  lines  at  various  compositions.  Using  this  information 
along  with  the  binary  information,  one  may  construct  ternary 
phase  diagrams.  The  major  features  of  a  ternary  phase  diagram 
are  established  by  the  tie  lines  between  compounds  and  may  be 
predicted  by  determining  which  pairs  of  compounds  have  a  lowest 
Gibbs  energy  at  cross-over  points.  However,  one  must  know  the 
AGf  values  of  all  the  compounds  present  in  the  system.  Since 
only  one  of  these  values  were  known  for  the  Ga-As-Cr  system 
(out  of  9  compounds),  and  because  experimentally  determining 
these  values  is  a  difficult  task,  requiring  specialized 
equipment,  this  method  was  not  used.  It  is  easier  to  anneal 
the  various  compositions  and  find  which  phases  are  present. 

Good  tutorials  and  discussions  of  ternary  phase  diagrams  and 
their  determination  are  found  in  Rhines'  Phase  Diagrams  in 
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Metallurgy  -  Their  Development  and  Application  (32 )  and  in 
Alper's  Phase  Diagrams  (33). 

Two  hundred  samples  of  various  Ga-As-Cr  compositions  (50 
different  compositions,  shown  in  figure  3-1),  lying  below  the 
50  a/o  arsenic  composition  line,  were  made  up  by  mixing 
appropriate  amounts  of  powdered  Cr,  Ga,  GaAs,  and  CrAs.  Some 
samples  of  the  same  composition  were  made  up  two  or  more  ways, 
i.e.,  using  GaAs,  Cr,  and  Ga ,  or  CrAs,  Cr,  and  Ga.  The  samples 
were  each  heated  in  evacuated  quartz  ampoules  to  1050 °C  and 
then  annealed  at  various  temperatures  (from  600  to  1050 °C)  for 
10  hours.  They  were  then  quenched  in  water.  (The  "room  temp" 
samples  were  slow-cooled  at  *100°C/hr.  to  room  temperature, 
instead  of  quenching.)  The  samples  were  then  crushed,  x-ray'd, 
reannealed,  and  re-x-ray'd.  By  comparing  the  x-ray  diffraction 
patterns  thus  obtained  to  the  known  patterns  of  the  elements 
and  compounds,  the  phases  present  for  each  composition  and 
temperature  were  determined.  (Mote:  no  ternary  compounds  were 
discovered.  All  peaks  on  the  x-ray  patterns  could  be 
attributed  to  binary  compounds  or  elements.  )  The  presence  of 
the  Ga-rich  liquid  (<*10  a/o  Cr,  20  a/o  As,  and  70  a/o  Ga  at 
700°C,  for  example)  could  be  identified  by  the  physical 
condition  of  the  annealed  sample  (i.e.,  very  malleable  and 
difficult  to  crush  into  a  powder  without  using  liquid  nitrogen 
to  cool  it  first)  and  by  the  presence  of  very  small  amounts  of 
Ga  and  CrGa4  as  seen  in  the  XRD  pattern;  I  experimentally 
determined  that  this  Ga-rich  liquid  becomes  CrGa4,  Ga,  and 
GaAs,  upon  quenching. 

Twenty-six  compositions  were  analyzed  by  DTA  (see  figure 
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3-2  for  the  compositions).  As  mentioned  earlier,  larger 
amounts  were  pre-annealed  and  x-ray'd  before  loading  them 
into  DTA  ampoules.  Each  went  through  3  heat/cool  DTA  cycles. 
These  were  again  x-ray'd  to  make  sure  that  the  same  phases  were 
present;  they  were  in  all  cases.  In  figures  3-3  and  3-4  are 
shown  the  DTA  curves  for  two  of  the  samples,  HI  and  B2 ,  whose 
compositions  are  shown  in  figure  3-2.  Table  3-1  lists  the 
transition  temperatures,  as  determined  from  the  DTA  work,  for 
the  17  ternary  compositions  -  the  remaining  9  were  binary 
compositions. 

3 . 3  Results  and  Discussion 

3.3.1  The  Cr-Ga  System 

The  phase  diagram,  as  proposed  by  Bornand  and  Feschette 
(14),  and  modified  slightly  by  Phillipe  et  al.  (16),  was 
confirmed  within  experimental  error.  The  x-ray  diffraction 
patterns  for  the  various  compounds,  as  reported  in  the  various 
references,  were  confirmed. 

3.3.2  The  Cr-As  System 

As  stated  earlier,  the  x-ray  diffraction  pattern  of  CrAs 
was  reported  previously  (20),  but  unambiguous  identification  of 
the  compound  was  not  possible  until  after  I  synthesized  the 
compound  and  produced  a  diffraction  pattern  showing  the 
indwidual  peaks.  The  diffraction  pattern  I  obtained,  shown  in 
figure  3-5,  with  the  corresponding  (hkl)  indices  listed  in 
table  3-2,  is  consistent  with  the  reported  (20)  lattice 
parameters  and  structure  (orthogonal)  of  CrAs.  The  lattice 
parameters  of  CrAs,  as  determined  from  the  angles  of  the 
various  reflections,  were  found  to  be: 
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Table  3-1.  Transition  Temperatures  for  the  17  Ternary 
Compositions. 


Sample  Designation 


Transition  Temps . ( °C ) 


1150,  1115,  670 
1160,  1060,  670 
1160,  970,  670 
1070,  910,  870,  770,  675 
1040,  960,  880,  770,  740 
990,  980,  880,  770,  740 
1010,  990,  880,  770,  740 
1320,  1270,  1015,  880,  740 
1070,  670 

1060,  870,  770,  675 
1140,  1100,  670 
1190,  1170 
1050,  975,  670 


a  =  3.48  A  (±0.01  A) 
b  =  6.22  A 
c  =  5.75  A 

These  compare  favorably  with  Nowotny  and  Arstad's  (20): 
a  =  3.486 
b  =  6.222 
c  =  5.741 

As  mentioned,  I  did  get  good  separation  between  peaks;  for 
instance,  between  the  (022)  and  (112)  plane  diffractions. 

The  melting  points  of  the  Cr-As  compounds  were  determined 
to  be  the  following  (all  ±  5°C): 

CrAs:  1360°C 

Cr4As3:  1185°C  ( incongruently  to  CrAs  +  j l) 

1300 °C  (to  all  liquid) 

Cr2As:  1175°C  (no  high  temperature  phase  was 

observed ) 

Cr4As:  >1400°C  (the  limit  for  quartz  ampoules) 

The  diffraction  patterns,  reported  by  others,  of  the  other 
Cr-As  compounds  were  confirmed. 

I  was  unable  to  determine  if  "Cr4As"  was  Cr4As  or  Cr3As, 
or  if  "Cr4As3"  was  Cr4As3  or  Cr3As2,  except  that  in  assuming 
the  former  stoichiometry  in  each  case  (based  on  the  most 
current  reports  in  the  literature),  and  adding  appropriate 
amounts  of  Cr  to  CrAs  on  that  basis,  I  observed  only  the  single 
phase  after  annealing  (by  x-ray  analysis),  leading  me  to 
suspect  that  these  assumptions  were  correct. 

With  regard  to  Yuzuri's  estimated  phase  diagram,  I  did 
agree  with  him  regarding  which  compounds  melt  congruently 
(Cr4As,  Cr2As,  and  CrAs)  and  which  incongruently  (Cr4As3),  but 
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the  actual  melting  points  differ  considerably.  The  relatively 
high  melting  point  of  CrAs  that  I  determined  would  lead  one  to 
predict  it  as  being  a  dominant  compound  in  the  Ga-As-Cr  s/'tem, 
an  unlikely  prediction  based  on  Yuzuri's  results.  Because  of 
this,  and  because  of  the  problems  with  the  reported  XRD  pattern 
of  CrAs,  it  was  only  after  I  concluded  the  work  on  the  Cr-As 
system  that  I  was  able  to  make  sense  of  my  ternary  phase 
diagram  results. 

3.3.3  The  Ga-As-Cr  Ternary  System 

From  the  200  anneals  and  x-ray  results,  and  the  26  DTA 
results,  isothermal  sections  from  room  temperature  to  1300 °C 
were  constructed  for  the  Ga-As-Cr  system  (for  As  a/o  £50). 
Several  of  these  are  shown  in  figures  3-6  to  3-15.  The 
liquidus  and  solidus  lines  which  are  dotted  are  only 
approximate,  as  are  the  sizes  of  the  single  solid  phase  regions 
and  the  two-phase  tie-lines;  the  size  of  the  GaAs  phase  is 
exaggerated  for  clarity. 

In  addition  to  the  isotherms,  a  vertical  section  (an 
’’isopleth")  was  constructed  along  the  GaAs-Cr  composition  line, 
from  0  to  70  a/o  Cr,  and  is  shown  in  figure  3-16.  The 
experimentally  determined  points  are  included  in  this  figure. 
This  is  a  vertical  slice  through  the  three  dimensional  ternary 
phase  diagram  in  contrast  to  hoizontal  slices  for  isotherms. 
This  paticular  vertical  section  shows  what  happens  as  different 
amounts  of  Cr  are  added  to  GaAs.  GaAs  is  present  only  in  the 
regions  on  the  right  side  of  the  diagram,  and  Cr  and  GaAs  are 
never  present  in  the  same  region.  This  diagram  also 
illustrates  the  complexity  of  the  system. 
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Figure  3-7. 


Isothermal  section  of  Ga-As-Cr  phase  diagram  for 
670°C.  Dashed  lines  indicate  estimated  features 


700  °C 


Figure  3-8.  Isothermal  section  of  Ga-As-Cr  phase  diagram  for 
700°C.  Dashed  lines  indicate  estimated  features 
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Figure  3-10.  Isothermal  section  of  Ga-As-Cr  phase  diagram  for 
880°C.  Dashed  lines  indicate  estimated  features 
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Figure  3-11.  Isothermal  section  of  Ga-As-Cr  phase  diagram  for 
1000 °C.  Dashed  lines  indicate  estimated  features 
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Figure  3-13.  Isothermal  section  of  Ga-As-Cr  phase  diagram  for 
1150 °C.  Dashed  lines  indicate  estimated  features 


Cr3Ga 


CrGa  Cr3Ga4 


A 

CrGa4 


Ga 


Figure  3-14.  Isothermal  section  of  Ga-As-Cr  phase  diagram  for 
1200°C.  Dashed  lines  indicate  estimated  features 
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Figure  3-16.  Vertical  section  of  Ga-As-Cr  phase  diagram  along 
GaAs-Cr  line.  Compositions  of  DTA  samples  are 
indicated  at  top. 
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The  major  significance  of  how  these  results  relate  to 
solubility  and  diffusion  studies  of  Cr  in  GaAs  are  the 
following : 

1.  CrAs  and  C^Ga,  and  not  GaAs,  are  the  dominant 
compounds  in  the  Ga-As-Cr  system.  This  means  that  GaAs  is  not 
in  equilibrium  with  nearly  as  many  phases  as  was  previously 
presumed . 

2.  Cr  and  GaAs  are  not  in  equilibrium  with  each  other 
under  any  conditions;  there  is  no  tie  line  connecting  those  two 
phases  at  any  temperature.  Therefore,  a  solubility  or 
diffusion  experiment  in  which  Cr  and  GaAs  are  placed  in  contact 
with  each  other,  or  in  the  same  ampoule,  would  result  in  either 
the  Cr  phase  or  GaAs  phase,  or  both,  depending  on  where  the 
total  composition  lies  on  the  ternary  phase  diagram,  ultimately 
disappearing  during  the  anneal.  (This  is  in  contrast  to  what 
was  done  in  previous  studies  (5,34)  and  which  will  be  discussed 
in  more  detail  in  chapter  6. )  This  also  means  that  some  other 
phase  (CrAs(g)  or  Ga^j),  rather  than  Cr,  would  precipitate 
from  GaAs  if  the  Cr  content  exceeds  the  solubility  level. 

3.  There  is  a  large  tie  triangle  region  that  includes 
GaAs,  plus  CrAs  ( s )  and  Ga^j;  this  is  a  logical  choice  as  a 
Cr  source  for  solubility  and  diffusion  experiments.  As 
discussed  in  chapter  2,  as  long  as  the  entire  system  (GaAs  +  Cr 
source)  is  anywhere  in  that  region,  the  system  is  invarient  and 
well-defined  solubility  and  diffusion  experiments  are  achieved 
for  isothermal  conditions.  Also,  when  putting  this  composition 
in  an  ampoule  with  undoped  GaAs  (for  solubility  and  diffusion 


experiments),  the  system  would  initially  be  removed  from  the 
final  equilibrium  state  only  by  the  amount  of  Cr  that  can  be 
dissolved  in  GaAs.  This  is  the  "quasi-equilibrium"  condition, 
discussed  in  the  preceding  chapter,  which  is  necessary  for 
valid  solubility  and  diffusion  studies. 

4.  A  Ga-rich  liquid  is  in  equilibrium  with  GaAs.  (This 
corresponds  to  the  liquid  solvent  used  in  Liquid  Phase 
Epitaxial  (LPE)  growth  of  GaAs.)  This  liquid  has  a  very  high 
Cr  solubility,  up  to  =10  a/o,  which  is  consistent  with  the 
results  of  Andre  and  LeDuc  (30),  as  mentioned  in  section  3.1. 
The  significance  of  this  in  my  solubility  and  diffusion  studies 
will  be  discussed  in  the  upcoming  chapters. 

3 . 4  Summary  of  Phase  Equilibria  Studies 

In  this  part  of  the  study  of  chromium  in  gallium  arsenide, 
the  following  was  accomplished: 

1.  The  various  Cr-As  compounds  were  synthesized,  melting 
points  were  determined,  and  x-ray  diffraction  patterns  were 
confirmed  or  corrected; 

2.  The  Cr-Ga  phase  diagram,  including  the  same  type  of 
information  as  in  the  Cr-As  system,  was  confirmed; 

3.  The  Ga-As-Cr  ternary  phase  diagram  was  experimentally 
determined  for  As  <_  50  a/b,  from  room  temperature  to  1300°C. 

4.  Important  implications  of  these  results  on  Cr  in  GaAs 
solubility  and  diffusion  studies  were  deduced,  namely: 

a)  GaAs  is  not  a  dominant  compound  in  the  Ga-As-Cr 
system; 

b)  GaAs  and  Cr  are  never  in  equilibrium  with  each  other; 

c)  A  suitable  tie  triangle  region  was  found  for  use  in 


solubility  and  diffusion  studies; 

d)  A  Ga-rich  liquid,  which  has  a  high  solubility  of  C 
is  in  equilibrium  with  GaAs  at  the  temperatures  of  interest 


4.  THE  SOLUBILITY  OF  CHROMIUM  IN  GALLIUM  ARSENIDE 


4 . 1  Introduction 

Valid  and  well-defined  solubility  studies  can  be 
performed  once  the  ternary  phase  equilibria  information  is 
known.  One  can  then  correctly  apply  these  results  to  other 
studies,  such  as  diffusion  studies.  Failure  in  the  past  to 
recognize  the  importance  of  phase  equilibria  has  led  to 
improper  determinations  and  uses  of  solubility  values. 

In  this  chapter,  previous  reports  and  uses  of  "solubility" 
values  of  Cr  in  GaAs  are  reviewed.  Next,  a  short  review  of 
solubility  in  relation  to  the  phase  diagram,  as  discussed  in 
chapters  2  and  3,  is  given.  Then,  the  experimental  techniques 
for  determining  solubility  are  described,  and  the  results  and 
discussion  of  the  solubility  experiments  are  presented. 

Finally,  the  influence  of  the  type  of  Cr  source,  and  other 
conditions  for  the  solubility  experiments,  on  the  surface 
morphology  of  GaAs  is  described,  which  provides  evidence  for 
the  relevance  and  correctness  of  the  phase  equilibrium  results. 

4 . 2  Literature  Review  of  the  Solubility  of  Cr  in  GaAs 

In  the  literature  there  are  numerous  invalid  and 

inconsistent  results,  references  to  results,  and  uses  of  the 
results  regarding  the  solubility  of  Cr  in  GaAs.  Most  agree 
(26,  35,  36)  that  the  maximum  solubility  of  chromium  in  undoped 
gallium  arsenide  at,  or  just  below,  its  melting  point  (1238°C) 
is  about  2  x  lO-'-7  cm-^ ,  but  that  is  where  the  agreement 
ends. 

The  most  obvious  of  the  errors  committed  in  measuring  the 
"solubility"  of  Cr  in  GaAs  are  the  use  of  noi. -equilibrium 
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techniques.  Morkoc  et  al.  (37)  propose  a  solubility  of  Cr  in 
GaAs  of  1  x  1()16  at  530 °C  from  molecular  beam  epitaxy 

( MBE )  experiments.  This  value  is  orders  of  magnitude  higher 
than  any  others  reported  for  that  temperature.  Bonnet  et  al. 
(38)  state  that  the  "solubility  limit  of  Cr  in  GaAs  is  around  5 
x  1016  cm"3  '•  (temperature  not  known)  from  their 
metal-organic  chemical  vapor  deposition  (MOCVD)  results.  Linh 
et  al.  (39)  point  out  that  their  own  MBE  and  MOCVD  results  for 
Cr  solubility  differ  from  each  other  as  well  as  from  other 
reported  results.  Both  these  methods  of  layer  growth,  however, 
are  generally  regarded  to  be  non-equilibrium  techniques,  and 
are  probably  not  suitable  for  equilibrium  solubility  studies. 

Even  solubility  results  taken  from  liquid  phase  epitaxial 
(LPE)  growth  have  problems  associated  with  them.  Woodard's 
(29)  work  on  LPE  growth  of  Cr-doped  GaAs  is  an  often  quoted 
source  for  Cr  solubility  values.  Even  if  those  results  are 
accurate,  they  are  only  for  the  particular  case  of  GaAs  is 
equilibrium  with  a  Ga-rich  liquid  of  a  particular  Cr  content 
(0.83  a/o  in  Woodard's  case).  This  leads  to  the  most  common 
misuse  of  Cr  solubility  values:  blindly  applying  solubility 
values  from  one  particular  experimental  system  to  another, 
without  any  understanding  and/or  concern  of  ternary  phase 
equilibria  considerations.  A  classic  example  of  this  is  in  two 
articles  by  Feng  et  al.  (40,41).  They  use  Cr  concentration 
levels  from  Tuck  and  Adegboyega ' s  (5)  diffusion  results  as  the 
Cr  solubilities  for  various  temperatures  in  order  to  interpret 
their  Cr  outdiffusion  experiments.  They  also  claim  that 
Woodard's  (29)  results  are  consistent  with  these  values. 


There  are  three  errors  in  this.  First,  Tuck's  Cr  bulk  levels 
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came  from  diffusing  Cr  from  a  pure  Cr  source  into  GaAs  in  a 
closed  quartz  ampoule.  However,  it  has  been  determined  in  this 
work  that  chromium  is  not  in  equilibrium  with  GaAs,  at  any 
temperature.  Therefore,  Tuck's  experiments  were  not  performed 
under  either  equilibrium  or  "quasi-equilibrium"  conditions,  and 
the  concentration  levels  should  not  be  used  as  "solubility" 
levels.  Secondly,  the  Cr  levels  they  interpret  as  being  the  Cr 
in  GaAs  solubilities  at  various  temperatures  are  from 
indiffusions  that  were  not  performed  long  enough  to  reach  the 
final  value  in  the  bulk  for  all  but  one  temperature;  this  point 
was  made  clear  in  Tuck's  paper.  Finally,  comparing  these 
results  to  Woodard's  is  in  error  for  the  reason  mentioned 
above,  namely,  one  must  consider  the  Cr  source  itself.  Another 
example  of  inadequate  consideration  of  the  Cr  source  is  given 
by  Woodard  (29)  himself.  Woodard  cites  Huber's  (42) 
outdiffusion  results  and  interpretations  as  being  consistent 
with  his  own  solubility  numbers.  In  Huber's  experiments  as  in 
Feng's,  however,  the  GaAs  is  not  in  equilibrium  with  any  liquid 
phase,  let  alone  one  with  0.83  a/o  Cr.  Thus,  since  the  Cr 
source  conditions  were  different  in  these  studies,  the 
solubility  results  should  not  be  compared. 

There  are  many  examples  in  which  the  results  are  not 
consistent  with  Woodard's  numbers,  and  the  researchers  can  not 
figure  out  why.  For  example,  Linh  et  al.  (39)  look  at  the  Cr 
levels  in  their  vapor  phase  epitaxy  (VPE)  layers  and  conclude 
that  they  are  about  four  times  higher  than  Woodard's  results 
would  suggest.  They  were  "unable  to  resolve  this 
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discrepancy".  Again,  they  assumed  that  Woodard's  solubility 
values  were  for  all  conditions. 


These  examples  show  the  general  lack  of  concern  for 
ternary  phase  equilibria  when  reporting  and  using  Cr-in-GaAs 
solubility  data.  It  is  shown  in  this  thesis  how  critical  it  is 
to  have  the  system  at  (or  near)  equilibrium  and  to  be 
well-defined  for  meaningful  and  valid  solubility  (and 
diffusion)  data. 

4 . 3  Solubility  in  Relation  to  the  Ternary  Phase  Diagram 
As  discussed  in  chapters  2  and  3,  there  are  several 
features  of  ternary  phase  equilibria  that  are  important  in 
solubility  studies  of  an  impurity  in  a  compound  semiconductor. 
In  summary,  they  are:  1)  the  solubility  of  a  component  in  a 
phase  (Cr  in  GaAs  in  this  work)  depends  on  the  location  of  the 
system  on  the  ternary  phase  diagram;  2)  the  solubility  is 
well-defined  if  the  temperature  is  specified  and  the  system 
lies  in  a  three  condensed  phase  region  (i.e.,  a  tie  triangle); 
and  3)  a  composition  within  a  three  phase  region  containing 
GaAs  is  a  suitable  source  for  solubility  experiments  of  Cr  in 
GaAs.  Specifically  a  composition  near  the  middle  of  the 
GaAs  ( s  )-CrAs  ( s  )-Ga  )  tie  triangle  region  will  be 
utilized  extensively  as  a  Cr  source  for  the  solubility  studies 
in  this  work.  Not  only  does  it  lead  to  well-defined 
conditions,  but  the  region  covers  a  very  large  composition 
region  around  the  GaAs  phase,  therefore  the  results  are 
applicable  to  many  experimental  situations.  In  addition,  the 
Cr  solubility  in  this  region  is  expected  to  be  greater  than 
anywhere  in  the  GaAs ( s )  -  Ga/^)  region  (which  is  where  the 


only  other  results  have  come  from);  therefore,  the  results  one 
obtains  (by  SIMS,  etc.)  will  be  statistically  much  more 
reliable.  Finally,  it  should  be  reiterated  that  elemental  Cr 
should  not  be  used  as  a  source  for  GaAs  studies  since  it  is 
never  in  equilibrium  with  GaAs.  (The  effect  of  doing  this  will 
be  demonstrated. ) 

The  above  statements  apply  not  only  to  solubility,  but  to 
all  the  other  intensive  properties  as  well,  such  as  partial 
pressures,  equilibrium  defect  concentrations,  and  related 
properties  such  as  diffusivity,  and  will  be  applied  to  the 
study  of  the  diffusion  of  Cr  in  GaAs  and  also  to  the  defect 
chemistry  of  GaAs  in  the  upcoming  chapters. 

4 . 4  Experimental  Techniques 
4.4.1  General  Techniques 

For  most  of  the  experiments,  the  GaAs  used  was  Cominco 
LEC  high  resistivity,  undoped  GaAs.  According  to  Cominco,  the 
resistivity  was  _>  1°7  ohm  cm.  and  the  dislocation  density  was  £ 
104  cm-2.  This  material  was  used  in  most  cases,  rather 
than  Bridgman  GaAs,  for  example,  because  of  its  very  low 
impurity  level,  which  is  the  reason  why  it  is  semi-insulating 
without  Cr  additions.  Other  GaAs  material  used  included: 
Crystal  Specialties  undoped  GaAs  (Bridgman);  Hewlett-Packard 
Te-doped  and  undoped  GaAs  (LEC);  and  Morgan  semiconductor 
Si-doped  GaAs.  The  GaAs  wafers,  350  to  400  microns  thick,  were 
cut  into  approximately  1  x  1.5  cm  squares.  Before  annealing, 
the  squares  were  cleaned  according  to  the  following  schedule 
which  was  adapted  from  Hewlett-Packard's  GaAs  cleaning 


schedule: 


1.  Gently  wipe  both  sides  of  the  wafers  with  Q-tips  in  a 
aquet  detergent  solution. 

2.  Rinse  substrate  with  DI  water  several  times. 

3.  Rinse  substrate  with  acetone  and  ultrasonic  clean  for 
2  min. 

4.  Replace  acetone  with  trichloroe thane  and  ultrasonic 
clean  for  2  min. 

5.  Rinse  twice  with  acetone. 

6.  Rinse  twice  with  methanol. 

7.  Ultrasonic  clean  in  acetone  for  2  min. 

8.  Boil  in  hot  trichloroe thane  for  ten  minutes. 

9.  Rinse  twice  with  acetone  and  agitate  in  ultrasonic 
cleaner. 

10.  Rinse  twice  with  methanol  and  agitate  in  ultrasonic 
cleaner . 

11.  Boil  in  HC1  for  10  minutes. 

12.  Pour  out  most  of  HC1  and  replace  with  cold  HC1. 

13.  Rinse  with  methanol  several  times  and  ultrasonic  clean 
for  two  minutes. 

14.  Rinse  in  isopropyl  alcohol  and  boil  for  two  minutes. 
Quartz  beakers  and  wafer  holders  were  used  throughout.  The 
first  10  steps  are  to  remove  large  particles  (GaAs  dust,  etc.), 
to  remove  organic  impurities,  and  to  degrease  the  sample. 

The  HC1  step  is  to  remove  the  native  oxide  that  usually  forms, 
and  the  last  two  steps  remove  the  HC1. 

The  Cr  sources  used  were  appropriate  Ga-As-Cr  phases  which 
were  made  prior  to  the  solubility  anneals  by  mixing  appropriate 
amounts  of  Cr,  Ga ,  GaAs,  and  CrAs,  and  annealing  these 


mixtures  -  this  was  so  that  the  sources  were  already  in  their 
final,  equilibrium  state  before  the  actual  solubility 
experiments  began.  The  annealing  procedure  was  outlined  in 


chapter  3.  The  source  material  was  then  crushed,  and  small 
samples  of  it  were  x-rayed  to  be  sure  the  appropriate  phases 
were  present.  The  ampoules  used  were  made  from  high  purity 
quartz,  and  cleaned  in  the  same  manner  as  described  in  chapter 
3.  For  the  solubility  (and  diffusion)  experiments,  a 
constriction  was  made  in  the  ampoule  to  keep  the  powder  source 
material  and  the  the  GaAs  substrate  physically  separated.  The 
total  interior  volume  of  the  ampoule  was  *>  2  ml.  A  diagram  is 
shown  in  figure  4-1. 

Immediately  after  cleaning  the  GaAs  wafer,  the  wafer  and 
source  were  placed  in  the  ampoule,  which  was  evacuated  to 
»10-6  torr,  and  sealed.  The  vacuum  inside  the  ampoule  was 
checked  by  a  Tesla  coil  before  and  after  each  anneal.  The 
annealing  was  done  in  the  furnace  system  described  in  chapter 
3.  After  annealing  the  samples  for  the  appropriate  time  and 
temperature,  the  ampoules  were  quickly  removed  from  the  furnace 
and  quenched  in  water.  The  GaAs  wafers  were  then 
ultrasonically  cleaned  in  methanol  and  boiled  in  HC1  to  remove 
any  loose  particles  and/or  dissolve  any  precipiates  from  the 
surface,  and  then  rinsed  in  isopropyl  alcohol.  The  GaAs  wafers 
were  then  analyzed  by  SIMS.  The  following  section  will  discuss 
SIMS  in  detail,  first  discussing  the  basic  concepts  of  SIMS, 
and  then  presenting  the  specific  details  of  the  present  work. 
4.4.1  Secondary  Ion  Mass  Spectrometry 


Secondary  Ion  Mass  Spectrometry  (SIMS)  has  become  a 
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4-1.  Schematic  diagram  of  the  annealing  system  used 
the  solubility  and  diffusion  experiments. 


powerful  and  popular  technique  for  impurity  analysis  in 
semiconductors.  Its  high  sensitivity,  allowing  one  to  measure 
concentrations  on  the  order  of  1  x  10-^  cm"^  (ppb),  and  its 
ease  of  obtaining  concentration  vs.  depth  profiles  are  the  main 
reasons  for  this.  Good  reviews  of  SIMS  are  found  in  an  article 
by  V.R.  Deline  (43)  and  a  book  by  A.W.  Czanderna  (44).  A  good 
review  of  SIMS  applications  is  an  article  by  E.  Zinner  (45). 

The  method  uses  low  energy  ionized  atoms  to  bombard  a 
sample  surface,  and  by  chemical  and  kinetic  processes,  the 
surface  atoms  are  ionized  and  ejected  (or  "sputtered")  from  the 
sample.  These  "secondary  ions"  (as  opposed  to  the  "primary" 
ions  of  the  bombarding  beam)  can  be  analyzed  by  mass 
spectrometry.  The  resulting  ion  intensity  (in  counts  per 
sampling  time  and  as  measured  by  the  mass  spectrometer)  versus 
sputtering  time  data  can  then  be  easily  converted  to 
concentration  versus  depth  data  by  appropriate  calibration 
techniques,  as  will  be  discussed  l^ter. 

Because  a  mass  spectrometer  is  used,  only  atoms  (or  ions) 
of  a  specific  mass  are  counted,  and  because  of  the  great 
sensitivity  of  mass  spectrometers,  one  can  count  very  low 
numbers;  therefore,  SIMS  is  a  very  sensitive  and  specific 
analysis.  Even  if  two  species  have  the  same  mass,  such  as  P 
and  SiH,  one  can  distinguish  between  them  using  the  "voltage 
offset"  technique,  which  considers  that  ions  with  the  same  mass 
usually  have  different  kinetic  energy  distributions. 

The  primary  beam  is  usually  Ar+,  02+,  or  Cs+, 
depending  on  the  electronegativity  of  the  species  to  be 
analyzed,  and  the  energy  of  this  beam  is  relatively  low  (1-15 


KeV)  so  that  the  damage  to  the  matrix  material  is  minimized. 

Because  of  the  complexity  of  the  sputtering  process,  there 
are  several  limitations  on  the  resolution  of  both  the 
concentration  and  depth  determinations,  leading  to 
uncertainties  regarding  the  accuracy  of  the  values  and  profile 
shapes  obtained  by  this  method.  However,  there  are  several 
steps  that  can  be  taken  to  minimize  these  problems.  These  will 
now  be  discussed. 

An  easy  and  effective  way  to  convert  ion  intensity  to 
elemental  concentration  -  very  difficult  to  do  theoretically  - 
is  to  first  analyze  standard  samples  in  which  the  concentration 
of  the  desired  element  in  the  matrix  material  is  known  and 
obtain  a  conversion  factor  for  that  particular  system.  This 
first  requires  that  the  concentration  in  the  standard  be 
known.  This  is  easily  done  using  ion  implantation.  Here  the 

standard  is  made  up  by  implanting  a  known  total  "dose"  of  the 
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desired  element  into  the  matrix  material.  Then  by  integrating 
the  SIMS  ion  intensity  -  time  profile  and  comparing  to  the 
total  dose  (and  knowing  the  size  of  the  sampling  area),  one  may 
obtain  a  counts-to-atom  conversion  factor.  This  number  is 
multiplied  by  the  atomic  density  in  the  solid  of  one  of  the 
matrix  elements  (e.g..  As  in  GaAs )  and  then  divided  by  the  ion 
intensity  of  that  element  for  that  profile  to  get  a  modified 
conversion  factor.  This  takes  into  account  the  fact  that  the 
ion  intensity  value  -  which  is  in  "counts  per  sampling  time"  - 
depends  on  the  sputtering  rate  of  that  particular  analysis. 
Therefore,  each  time  an  analysis  is  made,  a  different 
conversion  factor  does  not  need  to  be  used;  one  simply 


measures  the  matrix  ion  intensity  value  at  the  end  of  the 
sputtering  and  uses  this  and  the  modified  conversion  factor  (or 
"sensitivity  factor")  to  make  the  conversion  from  ion  intensity 
to  concentration. 

The  sputtering  rate  affects  both  the  concentration  results 
and  the  depth  resolution,  and  there  are  several  problems 
associated  with  it;  yet,  again,  steps  may  be  taken  to 
alleviate  them.  One  problem  is  a  variable  sputtering  rate 
which  may  be  caused  by  such  features  as  inclusions, 
precipitates,  and  thin  films  (e.g.,  oxides)  in  the  matrix.  If, 
for  example,  there  is  a  region  of  enhanced  sputtering,  then  the 
measured  depth  in  that  region  will  be  smaller  and  the  measured 
concentration  will  be  higher  than  the  actual  numbers.  There  is 
no  reasonable  way  to  correct  for  this;  however,  one  can  tell  if 
this  is  occurring  by  simply  profiling  a  matrix  element  and 
seeing  if  it  remains  constant  or  not.  Another  problem  is  that 
as  the  SIMS  crater  (the  hole  created  by  the  sputtering  process) 
gets  deeper,  the  primary  beam  becomes  defocused,  and  the  crater 
area,  and  hence  sampling  area,  may  change.  In  addition,  atoms 
from  the  sides  of  the  crater  may  be  sputtered  and  be  counted  as 
coming  from  much  a  deeper  region.  These  two  effects  may  be 
minimized  by  using  a  laterally  controlled  beam  which  rasters 
out  a  large  crater,  and  by  using  an  aperature  in  front  of  the 
detector  to  sample  from  a  small  and  uniform  sized  area. 

Other  problems  include:  atoms  in  the  matrix  being  knocked 
deeper  into  the  matrix  by  the  primary  ions,  leading  to 
erroneous  profiles;  damage  created  by  the  primary  ions 
(usually  this  damage  is  uniform  and  hence  not  important);  and 
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oxide  layers  at  the  surface  which  usually  enhance  the  secondary 
ion  yield  by  making  the  atoms  easier  to  ionize.  Also,  a  common 
observation  is  a  small  dip  in  ion  intensity  near  the  surface  in 
the  profile,  which  is  believed  to  be  due  to  the  depth  being 
shallower  than  the  projected  range  of  the  primary  beam.  Under 
normal  conditions,  this  is  usually  -200A  -  so  for  deep  profiles 
this  is  not  a  problem. 

For  this  work,  a  Cameca  "IMS  3F"  SIMS  system,  at  Charles 
A.  Evans  and  Associates,  was  used.  For  Cr  in  GaAs,  an  02+ 
primary  beam  was  used,  with  an  energy  of  12.5  keV.  The  beam 
size  was  either  75  or  ,150  um  in  diameter  (it  varied  from  day  to 
day),  which  rastered  out  an  area  of  250x250  um  (or  500x500  um 
when  the  150  um  was  used)),  and  the  sampling  area,  controlled 
by  an  aperature  in  front  of  the  detector  was  75  um2  (or  150 
um2 ) .  The  usual  primary  beam  current,  which  determined  the 
sputtering  rate,  was  -6  uamps ,  resulting  in  a  sputtering  rate 
of  -50  A/sec. 

At  the  start  of  each  set  of  analyses,  a  "low  chromium" 
sample  was  analyzed  by  SIMS  to  determine  the  background  level 
for  that  day.  This  was  usually  *1  x  10^5  cm~3  and 
sometimes  less. 

The  time-to-depth  conversion  was  determined  by  measuring 
the  depth  of  the  SIMS  crater  after  each  analysis  with  a  Dektak, 
which  uses  a  very  sensitive  stylus  to  measure  step  distances, 
and  then  assuming  that  the  sputter  rate  is  constant  throughout 
the  profile  (confirmed  by  periodically  profiling  of  Ga  and/or 
As  to  establish  that  they  give  constant  concentration-depth 
profiles).  There  is  an  uncertainly  factor  in  the  absolute 


depth  measurement  of  15%.  In  addition,  the  depth  measurement 
was  taken  across  the  entire  SIMS  crater,  giving  a  depth  vs. 
distance  profile.  This  showed  the  crater  to  be  very  uniform. 
Micrographs  were  also  taken  (see  figure  4-2)  which  confirm  the 
uniformity  of  the  SIMS  sputtering  of  GaAs,  as  opposed  to  some 
matrix  materials,  which  do  not  sputter  uniformly. 

To  determine  the  ion  intensity-to-concentration  conversion 
factor,  the  ion  implantation/integration  method,  as  discussed 
previously,  was  used.  An  example  of  this  is  shown  in  figure 
4-3.  This  is  a  SIMS  profile  of  Cr  implanted  into  GaAs  with  a 
dose  of  1  x  1014  cm-2.  Using  the  integral  of  this  profile 
(calculated  by  the  SIMS  computer),  and  knowing  the  dose,  the 
size  of  the  sampling  area,  and  the  arsenic  ion  intensity,  the 
conversion  or  "sensitivity"  factor  was  determined.  This  was 
done  occasionally  to  ensure  that  it  did  not  change  appreciably 
over  the  span  of  the  research  (it  varied  by  less  than  30%).  In 
addition,  one  sample  (a  12  hr,  800°C  Cr  indiffusion  sample)  was 
occasionally  analyzed.  The  standard  deviation  (after  five  Cr 
analyses  over  2  years)  was  50%.  This  was  therefore  taken  to  be 
the  maximum  expected  error  in  the  concentration  values  obtained 
in  SIMS  profiles.  (See  the  appendix  for  a  full  discussion  on 
the  error  analyses. ) 

4 . 5  Solubility  Studies  -  Results  and  Discussion 

Anneals  of  GaAs  with  a  "D  source",  as  shown  in  figure  4-4, 
were  done  at  800°C  for  6,  12,  and  24  hours  in  order  to 
determine  the  length  of  time  needed  to  reach  the  Cr  solubility 
level  in  the  bulk  of  the  GaAs  (i.e.,  >  4  urns  from  the  surface); 
SIMS  Cr  profiles  indicated  that  12  hours  was  sufficient  time  to 
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Figure  4-3. 


Ion  intensity  versus  sputtering  time,  as  determined 
by  SIMS,  of  GaAs  implanted  with  Cr  at  the  dose  and 
energy  indicated  above. 


do  so  at  800°C.  The  Cr  solubility  was  thus  determined  to  be 
3  x  10-^  cm"3  (+  1.5  x  10^6  cm"3  )  at  800°C.  Eight 
anneals  were  then  done  at  that  temperature  using  sources  D,  H, 
I,  and  J  (see  figure  4-4)  for  6,  12,  or  24  hours,  as  indicated 
in  figure  4-5.  As  mentioned  earlier,  according  to  the  phase 
rule  and  the  Ga-As-Cr  ternary  phase  diagram,  GaAs  samples  in 
equilibrium  with  sources  D,  H,  and  I  should  have  the  same  Cr 
solubility,  while  a  sample  in  equilibrium  with  source  J  should 
be  different,  and  probably  less.  The  SIMS  Cr  profiles,  shown 
in  figure  4-5,  show  this  to  be  so.  The  bulk  Cr  levels  for  the 
D,  H,  and  I  samples  are  all  around  3  x  10-*-®  cm-^  ,  while 
that  for  J  is  »4  x  10-*-^  cm-^ .  (Note:  the  strange  behavior 
at  the  surface  will  be  discussed  in  later  chapters.)  These 
results  confirm  the  ternary  phase  diagram  findings  of  chapter  3 
and  also  confirm  the  applicability  of  the  phase  diagram  and  the 
phase  rule  for  this  type  of  study. 

These  results  explain  why  Woodard's  (29)  solubility  values 
are  low  compared  to  the  values  obtained  for  sources  D,  H,  and 
I.  As  discussed  earlier,  his  results  are  for  GaAs  in 
equilibrium  with  Ga ( ^ )  whose  Cr  content  is  0.86  a/o.  This 
would  correspond  to  approximately  point  W  in  figure  4-3.  This 
is  to  the  right  of  even  J  in  the  GaAs ( s )-Ga ( 4 )  region,  and 
hence  one  would  expect  a  lower  Cr  solubility  in  GaAs  for  that 
case. 

Solubility  anneals  were  then  done  with  H  source  (which 
tended  to  result  in  less  surface  build-up  of  Cr)  at  700,  750, 
900,  and  1000°C  in  order  to  obtain  Cr-in-GaAs  solubility 
values,  again,  for  the  GaAs { s )-CrAs ( s )-Ga ( ^  )  tie 
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triangle.  Thef"  were  done  for  18,  18,  8,  and  7  hours 
respectively  -  long  enough  to  reach  the  maximum  level  in  the 
bulk.  The  resulting  SIMS  Cr  profiles,  including  that  for 
800°C,  are  shown  in  figure  4-6.  The  Cr  solubility  vs.  1/T  is 

plotted  in  figure  4-7.  As  expected,  the  Cr  solubility 

extrapolates  to  the  maximum  Cr  solubility  of  2  x  lO^  cm"^ 
just  below  the  melting  point  of  GaAs. 

As  is  the  case  with  most  impurities  in  semiconductors,  Cr 
exhibits  retrograde  solubility  in  GaAs;  i.e.,  the  solubility 
increases  with  temperature  above  the  eutectic  point,  670°C  in 
this  case.  Swalin  (46),  using  the  results  of  Thurmond  and 
Struthers  (47),  gives  a  good  review  of  retrograde  solubility, 
including  a  good  explanation  for  why  it  is  common  in 

semiconductors.  Basically  it  is  due  to  the  large  positive 

heats  of  solution  for  the  impurities  in  semiconductors 
(because  of  the  strong  bonding  between  the  lattice  atoms,  and 
hence  the  relatively  small  solubilities  of  impurities,  compared 
to  solutes  in  metals).  It  is  shown  by  Thurmond  and  Struthers 
how  large  heats  of  solution  then  lead  to  retrograde  solubility 
based  on  the  equations  for  the  solidus  and  liquidus  lines. 

To  confirm  that  the  bulk  levels  reached  in  the  solubility 
experiments  correspond  to  equilibrium  conditions,  an  experiment 
was  done  in  which  a  GaAs  wafer  was  annealed  at  900°C  (with  H 
source)  for  a  long  time  (12  hours),  and  then  annealed  at  800 °C 
(again  with  H  source).  SIMS  analysis  was  done  after  each 
anneal.  The  Cr  profiles  are  shown  in  figure  4-8.  As  one  can 
see,  the  bulk  Cr  level  after  the  900°C  anneal  (curve  A)  was  at 
8  x  1016  cm"3 ,  the  determined  solubility  level  for  that 
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temperature.  The  bulk  Cr  level  after  the  subsequent  800°C 
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anneal  (curve  C)  is  lower,  =>4  x  10-*-^  cm"^ .  Also  shown 
(curve  B)  is  the  profile  of  the  original  800°C  solubility 
determination,  with  a  bulk  level  of  3.5  x  lO-^  cm-^ .  This 
shows  that  the  same  Cr  concentration  (within  experimental 
error)  was  reached  both  from  above  and  from  below,  confirming 
that  it  is  the  equilibrium  solubility  concentration. 

From  the  solid  solubility  results,  and  from  the  Ga-As-Cr 
phase  diagram  results  (specifically,  the  Cr  concentration  in 
the  Ga-rich  liquid  in  equilibrium  with  CrAs(S)  and 
GaAs(S)),  the  distribution  coefficient,  k,  of  Cr  in  GaAs  was 
determined  for  various  temperatures;  the  distribution 
coefficient  equals  [Cr ] s/ (Cr I ^ ,  the  solid  being  GaAS(s) 
and  the  liquid  being  Ga ( ^  )  in  this  case.  These  numbers,  as 
well  as  the  liquid  and  solid  solubility  numbers  from  which  they 
were  calculated,  are  tabulated  in  table  4-1.  The  distribution 
coefficients  vs.  1/T  are  plotted  in  figure  4-9.  The 
distribution  coefficient  at  800°C  was  also  determined  from 
sample  J  (GaAs  in  equilibrium  with  J  source),  and  is  indicated 
by  "J"  in  the  graph.  The  results  show  that  k  does  not  vary 
much  with  Cr  concentration;  even  though  the  solubilities  vary 
by  an  order  of  magnitude  between  J  and  D,  k  varies  by  only  a 
factor  of  two.  Also  plotted  in  figure  4-9  are  the  results  of 
Woodard  (29),  from  his  LPE  studies  (which  are  for  GaAs(S)  in 
equilibrium  with  a  Ga-rich  liquid  of  0.86  a/o  Cr 
concentration).  These  are  very  consistent  with  the  results  for 
higher  Cr  concentrations.  Therefore,  while  the  solubilities 
vary  as  the  total  composition  moves  across  the  phase  diagram, 
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the  ratio  of  solid  solubility  to  liquid  solubility  remains 
approximately  constant,  as  one  would  expect  for  dilute 
solutions.  In  addition,  the  distribution  coefficients 
determined  from  Cr-doped  GaAs  crystal  growth,  which  are  for 
temperatures  just  below  the  melting  point  of  GaAs,  are  shown 
(from  references  24-29).  As  one  can  see,  the  results  of  this 
work  fill  the  gap  nicely  between  the  lower  temperature  results 
of  Woodard  for  LPE  conditions  and  the  higher  temperature 
results  from  the  bulk  crystal  growers.  (Note:  upon  further 
examination  of  the  "crystal  growth"  values,  one  observes  the 
trend:  the  more  Cr  in  the  melt,  the  smaller  the  k  value.  This 
makes  sense  since  as  the  Cr  content  in  the  melt  increases,  the 
melting  point  of  GaAs  decreases,  and  it  is  at  that  temperature 
at  which  the  distribution  of  Cr  takes  place.  Therefore,  the 
temperatures  for  which  the  k  values  correspond  to  for  those  six 
cases  would  be  different.  Taking  this  into  account,  the  six 
points  should  not  be  plotted  at  the  same  temperature,  but  for 
decreasing  temperatures,  and  would  follow  the  trend  of  k 
decreasing  with  decreasing  temperature.) 

To  see  if  there  is  a  difference  in  Cr  solubility  between 
one  type  of  undoped  GaAs  and  another,  solubility  anneals  were 
done  using  GaAs  substrates  from  different  suppliers  (Cominco, 
Hewlett  Packard,  and  Crystal  Specialties).  The  differences 
between  the  resulting  solubilities  were  within  the  expected 
error  (t  50%)  of  the  experiments.  Therefore,  the  way  the  GaAs 
is  produced  apparently  has  no  appreciable  effect  on  the  final 
equilibrium  solubility. 

To  see  the  effect  of  doping  on  the  solubility  of  Cr  in 


GaAs,  solubility  anneals  (at  800°C  using  H  source)  were 
performed  in  Te-doped  (n-type)  GaAs.  After  an  18  hour  anneal, 
10  yms  of  the  GaAs  was  etched  off  each  surface  (using  a  16:1:1 
H20:H202:H2S04  solution).  This  was  done  to  ensure  that  any 
surface  artifacts  are  removed,  in  addition  to  removing  any  Te 
depletion  region  caused  by  Te  outdif fusion.  (However,  since 
the  diffusivity  of  Te  at  800°C  is  very  small,  this  was  not  a 
problem. )  SIMS  analysis  indicated  the  Cr  solubility  to  be 
8  x  10^6  cm~3 ,  almost  three  times  that  in  undoped  GaAs  at 
that  temperature.  This  is  consistent  with  the  model  of  Cr 
being  a  deep  acceptor  in  GaAs  (on  a  Ga  site ) (48-51 ) *  since  the 
presence  of  a  donor  usually  increases  the  solubility  of  an 
acceptor  (46,52),  based  on  charge  equilibrium  considerations. 
(Likewise,  the  presence  of  another  acceptor  usually  decreases 
the  solubility  of  an  acceptor. )  These  results  are  consistent 
with  those  of  Brozel  et  al.  (36)  who  found  in  their  crystal 
growth  experiments  that  the  Cr  solubility  increases  with 
increasing  Si  concentraion ,  the  Si  acting  as  a  donor.  C-V 
carrier  concentration  measurements,  using  gold  contacts  to  the 
GaAs  surface  and  an  H-P  automated  C-V  measuring  system,  were 
made  on  the  Te-doped  sample  both  before  and  after  the  Cr 
indiffusion.  Before  Cr  diffusion,  the  n-type  carrier 
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10^7cm~^  carriers  were  compensated  by  the  chromium.  If  the 
Cr  acts  as  a  double  acceptor,  as  claimed  by  Brozel  et  al.  (36) 
then  that  means  that  1/2  x  (1.8  x  10^7)  cm-^  ,  or 
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9  x  10^-®  cm"^  Cr  atoms  are  electrically  active.  This 
compares  very  well  to  the  measured  3  x  10^®  cm-^  chromium 
(the  difference  is  well  within  the  experimental  error).  This 
means  that  the  bulk  chromium  seen  in  the  SIMS  profiles  is  in 
solution,  and  are  not  precipitates  or  clusters.  It  is  also 
good  evidence  of  the  accuracy  of  the  SIMS  analysis  in  this 
work. 

Chromium  in  GaAs  is  interesting  in  regards  to  its 
electrical  properties;  Cr  is  said  to  be  a  deep  acceptor  with  an 
energy  level  (for  the  first  ionization  state)  just  above  midgap 
and,  therefore,  above  the  Fermi  level  in  otherwise  undoped 
GaAs.  Therefore  Cr-doped  GaAs  is  semi- insulating , and  not 
p-type,  even  at  high  temperatures.  (Occasionally  it  is  found 
to  be  p-type  (53);  however,  that  is  probably  due  to  other 
acceptors  being  present. )  But  if  free  electrons  are  present, 
due  to  the  presence  of  n-type  impurities  or  donors,  the 
chromium  atoms  will  accept  the  electrons  (one  or  two  per  Cr 
atom,  depending  on  the  electron  concentration  (36)),  and 
compensate  the  donors.  This  is  the  reason  why  Cr  is  put  into 
GaAs  in  the  first  place:  to  compensate  the  n-type  impurities 
and  render  the  GaAs  semi-insulating.  Even  if  too  much  Cr  is 
added,  the  GaAs  will  still  remain  semi-insulating  and  not 
become  p-type.  It  is  for  these  reasons  that  Cr  is  sometimes 
called  an  "electron  trap"  in  GaAs,  rather  than  an  acceptor. 

And  as  White  (2)  points  out,  the  normal  models  and  equations 
for  acceptors,  such  as  those  regarding  occupancy  and  the 
position  of  the  Fermi  level,  may  not  apply  for  a  species  such 
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However,  in  heavily  doped  GaAs  (n-type),  all  the  Cr  is 
probably  ionized  (36),  and  acts  more  like  a  normal  acceptor 
(e.g.r  its  solubility  increases  with  an  increase  in  the  n-type 
doping  level),  and  the  equations  that  are  used  to 
quantitatively  predict  the  increase  in  solubility  may  be  used 
as  a  first  order  approximation.  The  increase  in  the  solubility 
of  an  acceptor  in  the  presence  of  a  donor  is  given  by 
(54,46): 

Cs"  »  Cs'  (1  +  gn/nj_ )/( 1  +  g)  (4-1) 

where  Cs"  is  the  new  solubility  level,  Cs '  is  the 
solubility  in  the  intrinsic,  or  undoped  case,  n  is  the  electron 
concentration  (due  to  the  doping),  n^  is  the  intrinsic 
carrier  concentration,  and  g  is  the  degeneracy  factor,  (which 
according  to  Sze  (55),  is  equal  to  1/4  for  acceptors  in  Ge,  Si, 
and  GaAs).  Testing  this  for  the  case  discussed  above,  at  800°C 
(n^  =  7  x  10-1-6  cm-^  (13)),  a  Te  doping  of  3.8  x  lO^7, 
and  an  intrinsic  Cr  solubility  of  3  x  10^  cm~^ ,  one  would 
predict  from  this  equation  a  new  solubility  of  6  x  10^-® 
cm“3 .  The  actual  solubility,  from  SIMS,  is  8  x  10^6 
cm-3 .  a  similar  calculation  at  850 °C  (with  the  same  Te 
doping),  would  predict  a  solubility  of  8  x  10^6  cm-^  for 
the  intrinsic  case.  A  solubility  experiment  was  done  at  850°C 
in  the  same  way  as  in  the  previous  case,  and  the  SIMS  Cr 
profile  indicated  a  new  solubility  of  9  x  10^-®  cm"^.  As 
predicted,  at  high  temperatures  the  %  increase  is  not  as  great, 
due  to  the  fact  that  n^  increases.  As  a  first  order 


approximation,  equation  4-1  is  adequate. 


A  similar  experiment  was  done  in  n-type  Si-doped  GaAs  at 
300°C,  for  48  hours.  The  exact  Si  concentration  was  not  known 
(it  was  believed  to  be  in  the  low-to-mid  10^7’s  as  with  the 
Te  case),  but  the  qualitative  result  was  the  same.  The  Cr 
solubility  was  higher  than  in  the  undoped  case:  7  x  101® 
cm“3  vs.  3  x  10^6  cm“3 .  So,  even  though  the  donor  was  on 
a  different  site  (on  a  Ga  site,  rather  than  an  As  site  as  is 
the  case  of  Te),  the  effect  was  the  same,  indicating  that  it  is 
an  electrical  effect,  rather  than  a  steric  and/or  lattice 
stress  effect. 

4 . 6  Surface  Morphology  and  Different  Sources 

Visual  observations  of  the  GaAs  surface  after  solubility 
experiments  confirm  and  nicely  illustrate  the  applicability  of 
the  phase  diagram  studies  to  these  solubility  studies.  As 
mentioned  in  chapter  3,  a  proper  choice  of  the  Cr  source,  as 
determined  by  the  ternary  phase  diagram,  is  essential  for  valid 
solubility  and  diffusion  studies,  and  the  condition  of  the  GaAs 
surface  after  such  experiments  clearly  demonstrates  this. 

All  the  solubility  experiments  reported  so  far  in  this 
chapter  were  done  by  using  a  well-defined  Cr  source  which  was 
in  equilibrium  with  GaAs  (e.g.  H  source  or  J  source).  An 
optical  micrograph  of  a  typical  surface  after  a  solubility 
anneal  is  shown  in  figure  4-10.  This  was  after  a  2  hour,  800°C 
anneal  with  H  source  (shown  at  500X).  In  contrast,  an  anneal 
was  done  under  the  same  conditions,  except  using  pure  Cr  as  the 
source  (see  figure  4-11).  As  is  obvious,  the  GaAs  surface 
greatly  degraded  in  the  latter  case.  As  predicted  by  the  phase 


Figure  4-10.  Optical  micrograph  (500X)  of  GaAs  surface  after 
800  °C ,  2  hr.  Cr  indiffusion  anneal  with  H  source 


Figure  4-11.  Optical  micrograph  (200X)  of  GaAs  surface  after 

800°C,  2  hr.  Cr  indiffusion  anneal  with  elemental 
Cr  source. 


diagram,  Cr  and  GaAs  are  not  in  equilibrium  with  each  other, 
and  new  phases,  probably  CrAs(S)  and  Ga^j  (which  contains 
Cr  and  As),  will  form  from  the  starting  materials.  (In 

addition,  a  SIMS  analysis  showed  a  lower  Cr  level  in  the  bulk 

of  the  GaAs  than  in  the  H  source  case,  »5  x  10 5  Vs.  »3  x 

10^5  cm~3 .  this  is  perhaps  due  to  the  GaAs  being  in 

contact  with  Ga^j.)  These  results  confirm  the  findings  of 
the  phase  diagram  determinations,  and  show  why  pure  Cr  should 
not  be  used  as  the  Cr  source  in  GaAs  solubility-diffusion 
studies,  as  has  been  done  in  the  past  (which  will  be  discussed 
in  chaper  6).  More  concerning  surface  phenomena  will  be 
presented  in  chapter  7. 

4 . 7  Summary 

In  this  chapter,  the  work  on  the  solubility  of  Cr  in  GaAs 
is  presented.  It  was  shown  that  previous  studies  paid  little 
attention  to  phase  equilibria  considerations  in  regards  to 
solubility  values,  both  in  determining  them  and  in  using  them. 
Using  the  information  acquired  from  the  phase  equilibria 
experiments,  valid  and  well-defined  solubility  studies  were 
performed  for  Cr  in  GaAs  for  a  range  of  temperatures  and  for 
different  regions  of  the  Ga-As-Cr  phase  diagram.  The  Cr 
solubility  in  GaAs  was  determined  from  700 °C  to  1000°C  for  the 
GaAS(S)  -  CrAS(S)  -  Ga ( 4  )  tie  triangle  region.  In 
addition,  the  solubility  was  determined  at  800 °C  for  GaAs  in 
equilibrium  with  a  Ga-rich  liquid  which  contained  5  a/o  Cr,  for 
comparison;  as  expected,  the  Cr  solubility  was  less  for  this 
case.  The  distribution  coefficient  was  determined  for  these 
same  temperatures,  and  the  results  were  consistent  with 


literature  values  for  lower  and  higher  temperatures  and  for 
different  Cr  concentrations  in  the  Ga-rich  liquid.  The  Cr 
solubility  in  Te-doped  (n-type)  GaAs  was  measured,  which  was 
found  to  be  higher  than  in  undoped  GaAs,  consistent  with  the 
model  of  Cr  acting  as  an  acceptor,  or  "electron  trap".  C-V 
carrier  concentration  measurements  of  the  Te-doped  GaAs,  both 
before  and  after  Cr  solubility  anneals,  confirmed  that  the  Cr 
acts  as  a  double  acceptor  for  high  donor  concentration  and  that 
the  Cr  measured  by  SIMS  is  in  solution  and  at  the  levels  that 
the  SIMS  analyses  indicate.  Finally,  observations  of  the 
surface  morphology  after  solubility  experiments,  using  a 
well-defined,  "quasi-equilibrium"  source  in  one  case,  and  pure 
Cr  in  another,  confirm  the  findings  of  chapter  3  concerning 
phase  equilibria  and  show  why  such  considerations  are  important 
in  solubility  studies. 


HE  SUBST ITUT 1 ONAL/ I NTERST IT  I AL  DISSOCIATIVE  DIFFUSION  MECHANISM 


5 . 1  Introduction 

The  substitutional/interstitial  dissociative  diffusion 
mechanism  has  been  effective  in  describing  diffusion  of 
impurities  in  semiconductors.  In  the  past,  it  was  used  mostly 
for  impurities  in  compound  semiconductors,  but  it  is  now  being 
applied  to  impurities  in  elemental  semiconductors,  including 
silicon.  This  mechanism  has  been  effective  in  explaining  a 
variety  of  complex  diffusion  profiles  which  are  observed  for 
different  impurity /matrix  material  systems  and  also  for 
different  profiles  in  the  same  system  that  arise  for  different 
diffusion  conditions.  This  chapter  is  a  brief  tutorial  on  this 
mechanism.  The  motivation  for  the  model  and  a  description  of 
its  development  is  presented  in  the  "History  of  the  Model" 
section,  along  with  a  listing  of  the  systems  that  are 
successfully  described  by  the  model.  The  governing  reactions 
and  equations,  and  the  different  cases  -  which  are  used  to 
simplify  the  rather  complex  equations  -  are  given,  followed  by 
a  discussion  of  how  several  real  systems  fit  these  different 
cases.  In  addition,  one  step  of  this  diffusion  process, 
namely,  the  substitutional-to-interstitial  conversion,  is 
examined,  and  two  mechanisms  for  this  step  are  presented  and 
discussed.  Finally,  activation  energies  and  D0  values  are 
examined  in  order  to  help  one  identify  the  mechanism  in  a 
specific  system  and  in  determining  which  is  the  appropriate 
case  for  a  particular  condition  in  a  system.  In  chapter  6, 
this  mechanism  is  used  to  model  the  complex  behavior  of 


chromium  in  gallium  arsenide. 

5 . 2  History  of  the  Model 

The  substitutional/interstitial  dissociative  diffusion 
mechanism  was  conceived  and  developed  in  response  to  unusual 
diffusion  behavior  observed  in  several  systems.  It  not  only 
explains  uniform  rises  or  drops  in  the  bulk  impurity 
concentration  level,  and  even  “uphill  diffusion"  at  the  back 
side,  but  also  explains  what  appears  to  be  anomalously  high 
diffusion  values  for  impurities  that,  according  to  their 
electrical  activity,  are  predominantly  located  on 
substitutional  sites.  It  also  explains  different  diffusion 
profiles  for  different  conditions,  and  even  two  or  more 
profiles  superimposed  on  each  other. 

Frank  and  Turnbull  (56),  in  1956,  first  proposed  the  basic 
mechanism  while  analyzing  other  workers'  results  of  diffusion 
of  Cu  in  Ge.  In  the  diffusion  profiles  (concentration  vs. 
depth)  of  Tweet  and  Gallagher  (57)  for  Cu  diffusing  into  Ge 
from  one  side,  it  was  noted  that  the  distribution  corresponded 
to  Fick's  Law  (normal  error  function  diffusion)  only  at  small 
depths.  At  greater  depths  there  was  a  concentration  of  copper 
which  was  virtually  independent  of  distance  and  which  rose 
uniformly  with  time.  In  addition,  the  copper  that  was  measured 
was  all  electrically  active  and  corresponded  to  substitutional 
copper  (i.e.,  acceptors).  Frank  and  Turnbull  showed  that 
dislocation  pipes  could  not  explain  the  behavior  (which  had 
previously  been  guessed  to  be  the  cause).  They  proposed  the 
following  mechanism.  Copper  dissolves  in  germanium  in  two 
states,  interstitial  and  substitutional.  In  interstitial 


solution,  it  has  a  very  high  diffusivity,  but  a  low 
solubility.  Conversely,  as  a  substitutional,  it  has  a  low 
diffusivity,  but  a  high  solubility.  However,  they  do  not 
diffuse  independently  as  in  the  normal  two-stream  model.  There 
is  a  dissociative  mechanism  corresponding  to  the  inter¬ 
conversion  of  interstititial  and  substitutional  species.  This 
probably  requires  a  supply  of  vacancies,  which  is  available 
only  at  the  free  surface  or  near  dislocations.  The  basic 
diffusion  mechanism  begins  with  copper  substitutional  atoms 
dissociating  near  the  surface  to  vacancies  and  interstitials, 
which  then  diffuse  independently.  Then,  the  interstitital 
atoms  convert  back  to  substutitional  (to  maintain  local 
equilibrium  between  the  two  species)  by  combining  with  a 
vacancy  which  could  come  either  by  diffusing  from  the  surface, 
or  by  generation  in  the  bulk  (from  dislocation  climD).  The 
first  case  accounts  for  the  normal  error  function  profile  at 
smaller  depths,  and  the  second  case  accounts  for  the  nearly 
uniform  level  in  the  bulk  which  rises  with  time. 

Since  the  mechanism  was  first  proposed,  it  has  been 
applied  effectively  to  the  following  impurity/matrix  systems: 
indium  in  GaAs  (58),  Ag  in  InP  (59),  Zn  in  InSb  (60),  Cd  in  InP 
(61),  Au  in  InAs  (62),  Au  in  Si  (63-65),  P  in  Si  (66),  and  most 
extensively,  Zn  in  GaAs  (10,  67-76).  In  his  review  of 
diffusion  in  III-V  semiconductors  (77),  Kendall  concludes,  upon 
a  study  of  the  literature,  that  the  following  systems  also 
follow  this  mechanism:  Cu,  Ag ,  Au,  Mn,  Sn,  Mg,  Cd,  and  perhaps 
S,  in  GaAs,  Cu  in  InSb,  and  Au  in  InP.  The  mechanism  has  also 
been  proposed  for  Cr  diffusion  in  GaAs  (5,78).  Several  of 


these  cases  will  be  discussed  in  more  detail  later. 

Almost  all  these  diffusing  elements  are  either  shallow  or 
deep  acceptors  in  III-V  compound  semiconductors,  and  are 
believed  to  occupy  the  group  III  element  site  in  the  compound. 
There  is  a  structure  motivated  rationale  for  a 
substitutional/interstitial  diffusion  mechanism  impurities  in 
III-V  compounds  since  an  impurity  on  a  group  III  site  must 
travel  through  a  group  V  plane  in  reaching  another  group  III 
site.  This  is  probably  done  interstitially  (79).  This 
mechanism  is  likely  whenever  a  slower  diffusing  species  has  a 
higher  solubility  than  a  faster  diffusing  species,  and  there  is 
a  finite  interconversion  rate  for  the  two  species.  Impurities 
in  semiconductors  often  have  these  properties. 

5 . 3  Governing  Equations 

The  substitional/interstitial  dissociative  diffusion  model 
is  described  in  detail  by  Tuck  in  his  book,  Diffusion  in 
Semiconductors  (7),  who  followed  the  treatment  of  Sturge  (80). 
The  basic  reactions  and  equations  for  the  model,  from  these 
references,  are  as  follows. 

As  described  earlier  for  this  mechanism,  the  impurity  atom 
enters  the  diffusion  host  at  the  surface  where 
interstitial/substitutional  equilibrium  is  attained 
and  equilibrium  solubility  values  are  achieved.  These  species, 
along  with  vacancies,  which  arise  either  from  the  surface  or 
from  the  vacancy  generation  attending  the  formation  of 
the  interstitials  from  substitut: onals,  diffuse  in  the 
material.  It  is  assumed  that  substitutional  diffusion  is 
negligibly  slow  compared  to  the  interstitial  and  vacancy 
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diffusion.  The  concentrations  of  the  substitutional  impurity 
species,  the  interstitial  impurity  species,  and  the  vacancies, 
will  be  denoted  by  Cs,  Cj_,  and  Cv  respectively,  and  are 
functions  of  x,  the  distance  from  the  surface.  The 
equilibrium,  or  solubility,  values  of  these  are  denoted  by 
Cs ' ,  Ci',  and  Cv'.  The  basic  conversion  reaction  is: 

Ii  +  V  =■  is  (5-1) 

where  I  denotes  impurity  atom,  V  denotes  vacancy,  and  the 
subscripts  i  and  s  denote  interstitial  and  substitutional, 
respectively.  Inside  the  material,  the  interstitials,  which 
have  a  low  solubility,  and  vacancies  can  combine  by  this 
reaction,  resulting  in  the  impurity  becoming  substitutional. 
These  vacancies  can  either  be  the  ones  diffusing  from  the 
surface  -  or  wherever  they  were  created  during  the  formation 
of  interstitials  from  substitutionals ,  to  be  discussed  in 
chapter  6  regarding  "junction  diffusion"  -  or  can  be  generated 
in  the  bulk.  More  interstitials  can  then  diffuse  in  and 
convert  to  substitutional  species  until  the  high  substitutional 
solubility  level  is  reached  throughout  the  material.  Figure 
5-1  shows  this  schematically.  (Note:  here  we  are  considering 
uncharged  species.  For  charged  species,  equation  5-1  would 
include  either  holes  or  electrons.  This  will  be  considered 
later.  ) 

The  assumption  is  made  that  the  equilibrium  between 
substitutional  and  interstitial  atoms  is  always  maintained,  so 
that: 
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where  K  is  the  equilibrium  constant  for  equation  5-1.  This,  of 
course,  is  also  true  when  the  concentrations  are  at  their 
equilibrium,  or  solubility,  values,  so: 

K  =  Cs  '  /C  '  Cv  '  (5-3) 

Fick's  Law  for  two  species  diffuson  is: 


(  5-4  ) 


where  Dj_  and  Ds  are  the  dif f us ivities  for  the  interstitial 
and  substitutional  species  respectively. 

We  assume  Ds  =0,  so 


(  5-5  ) 


The  continuity  equation  for  vacancies  is: 
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where  the  first  term  on  the  right  is  the  indiffusion  of 


vacancies  from  the  surface,  the  second  term  is  the  loss  of 


V'  '  -  V-  r  1  J  - 


vacancies  due  to  the  interstitials  converting  to 
substitutionals,  and  the  third  term  is  the  bulk,  production  of 
vacancies  by  dislocation  climb  or  some  other  process,  assuming 
that  the  production  of  vacancies  is  proportional  to  the 
deviation  from  the  equilibrium  value. 

The  solution  to  these  equations,  in  order  to  determine  the 
concentration-depth  profiles  as  a  function  of  time,  cannot  be 
solved  explicitly.  However,  several  cases  can  be  considered  in 
which  appropriate  approximations  can  be  made  resulting  in 
simple  solutions.  In  each  case,  it  is  assumed  the  total 
impurity  concentration  equals  the  substitutional  atom 
concentration  (since  CS>>C^),  and  therefore  we  wish  to 
solve  for  Cs(x,t)  in  each  case. 

Case  I.  Vacancy  equilibrium  is  maintained  throughout  the 
crystal.  In  this  case,  whenever  a  vacancy  is  needed  for  the 
interstitial-to-substitutional  conversion,  it  is  availabla. 
Therefore  Cv  =  Cv'  and  equation  5-2  becomes:  = 

CS/KCV'.  This  is  substituted  in  equation  5-5,  leading  to: 


3C 

— i  (1  +  KC  '  ) 
at 


(5-7) 


Rearranging,  and  substituting  for  K  from  equation  5-3,  leads 
to: 
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This  is  just  normal  error  function-type  diffusion,  with  an 
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effective  diffusion  coefficient  of: 


C  . 


Def f =  Di 


(C  .  '  +  C  '  ) 
i  s 


=  "D 


fast 


( 5-8a ) 


(Designated  as  such  because  it  is  usually  much  higher  than  the 
other  effective  D  to  be  discussed  next). 

For  indiffusion  with  constant  surface  conditions  (with 
C(x,0)  =  0,  C  ( 0  ,  t )  =  Cs  * ,  and  a  semi-infinite  medium), 
equation  5-8  can  be  solved,  resulting  in: 


C  =  Cs  =  Cs'  erf c  (x/2 ( Dfastt ) 1/2 )  (5-9) 


Case  I  would  correspond  to  a  high  dislocation  density  (rapid 
vacancy  generation),  a  high  initial  vacancy  concentration,  or 
long  diffusion  times  (to  be  discussed  in  conjunction  with  Case 
IV). 

Case  II.  There  is  no  production  of  vacancies  in  the  bulk 
(k  =0),  so  that  all  vacancies  for  the  interstitial-to- 
substitutional  conversion  must  come  by  vacancy  diffusion. 

Since  this  is  slow  compared  to  interstitial  diffusion,  Cj_ 
will  be  assumed  to  be  Cj.'  for  all  x.  Therefore  equation  5-6 
becomes : 

—  (1  +  KC^  =  Dv  (5-10) 
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which,  upon  substitution  for  K,  results  in: 


9 


M 


c 

s 


C  ' 
s 


erf  c 


2/U 


A 


slow 


for  indiffusion,  where 


D 


slow 
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(5-12) 


This  gives  the  same  kind  of  concentration  profiles  as  in 
Case  I,  except  Dslow<<Dfast  (vacancy  diffusion  is  usually 
much  slower  than  interstitial  diffusion)  and  therefore  the 
penetration  of  the  impurity  is  much  shallower  for  any  given 
time.  Case  II  would  dominate  when  the  dislocation  density  is 
low  (little  or  no  vacancy  generation)  and  for  short  times/low 
temperatures  (before  vacancy  generation  can  appreciably  replace 
the  vacancies .  ) 

Case  III.  This  is  when  ki>0,  but  small.  Again,  the 
interstitial  diffusion  is  fast  enough  so  that  C^  =  C^’  for 
all  x.  But  in  contrast  to  Case  II,  there  is  some  vacancy 
generation  in  the  bulk.  Equations  5-3  and  5-6  result  in: 


3C  ■ 


at 


slow 


(5-13) 


where  9  =  ( C v *  +  Cs')/kiCv'.  This,  and  the  same 

boundary  and  initial  conditions  for  indiffusion,  can  be  solved 

by  Laplace  transforms,  giving: 


C  =  Cs '  ( 1-exp (-t/Q ) )  erfc ( x/2 (Dsiowt ) l/2 ) 


(5-14  ) 


This  is  plotted  for  different  t/3  values  in  figure  5-2  (from 
ref.  7).  cor  small  x,  this  is  approximately  Case  II  (vacancy 
diffusion  limited),  and  for  large  x,  C  *  Cs '( 1-exp (-t/a )  and 
the  diffusion  is  vacancy  generation  limited  and  independent  of 
x.  In  finite  crystals,  vacancies  can  diffuse  in  from  the 
opposite  side  of  the  crystal  as  well.  According  to  this  model, 
the  interstitials  could  convert  to  substitutionals  there,  and 
this  results  in  a  U-shaped  profile,  even  when  the  impurity  is 
diffused  in  from  one  side  only.  This  "uphill  diffusion"  has 
become  a  test  for  the  substitutional/interstitial  dissociative 
diffusion  mechanism,  as  disscussed  by  Kendall  (81).  A  classic 
example  of  this  is  Au  in  Si  (7).  (Note:  this  is  not  uphill 
diffusion  in  the  sense  of  diffusing  up  a  concentation  gradient 
while  diffusing  down  an  activity  gradient,  rather,  it  is  uphill 
diffusion  in  the  sense  that  the  net  profile  appears  to  be 
diffusing  uphill;  the  profile  is  developed  by  interstitials 
actually  diffusing  down  their  concentration  gradient, 
andconverting  to  substitutionals  more  extensively  at  the 
surfaces  of  high  vacancy  concentrations.) 

Case  IV.  This  is  for  k]_>0,  and  large  O  small).  Here, 
the  generation  of  vacancies  is  fast  enough  so  that  the 
interstitial-to-substitutional  conversion  in  the  bulk  starts  to 
take  place  before  the  interstitial  concentration  is  uniform 
throughout  the  bulk.  If  we  assume  that  Dv  is  so  slow  that  it 
can  be  ignored,  and  also  assume  a  quasi-equilibrium  in  which 
3Cv/3t  =  0  in  the  bulk,  equation  5-6  becomes: 


3C  C 

3  =  Jl_  (C  1  -  _ !_) 

KCl  5  erf c  u 


where  u  =  x/2 ( Dfastt ) .  Using  equation  5-3  and 
approximating  9  by  Cg'/kCv',  gives: 
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(5-17) 


Simplifying  and  solving  yields: 


C  =  Cs '  erfc  u  (1  -  exp(-t/(e  erfcu)) 


(5-18) 


for  indiffusion.  This  is  plotted  in  figure  5-3  (from  ref.  7). 
As  t  gets  large,  this  approaches  Case  I.  If  there  is 
appreciable  vacancy  diffusion  at  the  surface,  then  the 
expression  C  =  Cs  erf ( x/2 ( Ds^owt ) )  would  apply  near 
the  surface,  and  would  be  superimposed  on  equation  5-18.  Or, 
as  Tuck  suggests,  there  may  be  a  region  near  the  surface  where 
vacancy  equilibrium  is  achieved  (due  to  a  high  dislocation 
density  from  surface  damage  or  from  the  impurity  atoms 
themselves)  and  Case  I  would  apply  for  that  region,  with  Case 
III  or  IV  obeyed  in  the  deeper  regions. 

Another  case,  suggested  by  Asbeck  (82),  is  similar  to  Case 
I  (where  there  is  no  vacancy  generation  -  or  vacancy 
annihilation  in  the  case  of  outdif fusion ) ,  except  that  both 
vacancy  and  interstitial  diffusion  are  both  limiting  factors. 
Asbeck  derives  an  effective  D  of: 


D 


(5-19) 
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which  is  smaller  than  Dsj_ow  above,  and  is  concentration 
dependent. 

The  different  cases  above  can  all  be  easily  applied  to 
outdif fusion.  The  governing  differential  equation  is  the  same 
in  each  case,  but  with  different  boundary  and  initial 
conditions.  The  solutions  are  error  functions  rather  than 
complementary  error  functions,  and  a  vacancy  annihilation,  or 
sink,  term  is  used  rather  than  a  vacancy  generation  term. 

For  the  four  cases  described  by  Tuck,  the  diffusion 
coefficients  are  concentration  independent;  they  are 
proportional  to  equilibrium  concentration  values,  which  are 
constants.  Implicit  in  these  analyses  is  the  assumption  that 
the  activity  coefficient  of  the  impurity  is  constant  with 
concentration  and  equal  to  unity  (Henry's  Law).  This  is 
reasonable  for  most  impurities  in  semiconductors,  especially  Cr 
in  GaAs,  since  the  concentrations  are  so  low.  If  the  activity 
coefficient  was  not  constant,  concentration  dependant  D’s  could 
result.  In  the  next  section,  it  will  be  seen  how  concentration 
dependent  D's  can  result  when  the  diffusing  species  are 
ionized . 


5 . 4  Concentration  and  Pressure  Dependences 

The  above  treatment  assumed  that  the  diffusing  species  are 
all  uncharged  and,  as  a  result,  the  effective  diffusion 
coefficients  were  independent  of  concentration.  However,  if 
one  or  more  of  the  species  is  charged,  then  concentration 
dependent  diffusion  can  result.  For  the  case  in  which  the 
interstitial  atom  is  a  singly  ionized  donor  and  the 
substitutional  atom  is  a  singly  ionized  acceptor,  for  example. 


equation  5-1  becomes: 


I . +  +  V  =  I  +  2h+  (5-20  ) 

1  s 

and  the  mass  action  law  is: 

KC . +  C  =  C  "  ( yP)2  (5-21) 

where  p  is  the  hole  concentration  and  y  is  the  hole  activity 
coefficient.  For  high  doping  (p>n^)  and  all  species  ionized, 
ci  =  ci'»  Cs  =  Cs ' ,  and  p  =  Cs  (y=l),  which  leads  to: 

K  CjCy  =  Cg3  (5-22  ) 

and 

K  C i '  Cv 1  =  Cg '3  (5-23  ) 


which  replace  equations  5-2  and  5-3.  Equations  5-5  and  5-6 
would  remain  the  same.  For  Case  I  (vacancy  equilibrium),  this 
leads  to 


5Cs  .  i_  ,D  !!• 

-  —  1  eff  - 

at  ax  ax 


where 


Deff  *  3  Ci-DiCs2/cs>3 


(5-24  ) 


(5-25) 


Therefore,  the  diffusion  coefficient  is  concentration  dependent 
(due  to  the  Cs2  term). 

It  is  important  to  note  that  it  is  possible  to  have 
charged  species  and  not  have  concentration  dependent 
diffusion.  This  would  be  the  case  if  in  the  conversion 
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reaction  equation,  I  j_  +  V  =  Is,  the  charges  on  each  side  of 
the  reaction  exactly  cancelled  each  other,  and  electrons  or 
holes  were  not  necessary  to  balance  the  equation.  One 
possibility  would  be:  I  j_  +  V“  =  is--  Chang  and 
Pearson  (83)  do  suggest  that  in  GaAs ,  the  gallium  vacancy  is  a 
singly  ionized  acceptor.  However  they  may  have  been  looking  at 
an  impurity-vacancy  complex,  thus  making  it  appear  that  the 
vacancy  was  charged.  Furthermore,  when  modelling  Zn  diffusion 
in  GaAs,  the  gallium  vacancy  is  assumed  to  be  neutral  and  the 
theoretical  profiles  fit  the  data  very  well;  if  the  vacancies 
were  ionized,  the  model  would  not  fit.  In  regards  to 
substitutional  chromium,  there  is  considerable  data  confirming 
that  substitutional  chromium  is  un-ionized  in  undoped  GaAs. 
Another  possible  reaction  is:  Ij_+  +  V  =  Is.  The 
interstitial  Zn  species  is  assumed  to  be  a  donor  (due  to  the 
d-electrons  that  are  available),  so  it  is  reasonable  to  think 
that  the  interstitial  Cr  atom  might  also  be.  Since  no 
electrons  or  holes  appear  in  this  reaction,  this  too  would 
result  in  diffusion  that  is  independent  of  concentration.  One 
more  possible  reaction  is:  e”  +  Ij_+  +  V  =  Is.  Here,  an 
electron  _is  needed  to  balance  the  reaction,  but  concentration- 
independent  diffusion  would  probably  still  occur.  This  is 
because  the  source  of  electrons  would  be  the  interstitials  and 
thus  the  electron  concentration  profile  would  be  virtually 
uniform  and  the  electron  concentration  level  would  be  very  low, 
probably  much  less  than  the  intrinsic  carrier  concentration 
level,  and  therefore  would  not  produce  concentration-dependent 
diffusion.  Considering  all  this,  impurity  diffusion  in  GaAs 
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that  fits  the  "uncharged"  model  rather  than  the  "charged"  model 
probably  means  that  the  substitutional  species  is  uncharged. 

Another  important  point  is  that  if  the  impurity 
concentration  is  much  below  the  intrinsic  carrier  con¬ 
centration,  then  concentration-independent  diffusion  will 
result  no  matter  if  the  species  are  charged  or  not.  The  reason 
for  this  is  because  the  hole  and  electron  concentrations  are 
dictated  by  ni  and  not  by  the  substitutional  impurity 
concentration . 

The  presence  of  charged  species  not  only  affects  the 
diffusion  by  changing  the  relative  concentrations  of  the 
different  species,  but  also  by  creating  an  electric  field. 
Charged  species  -  both  atomic  and  electronic  -  diffusing  at 
different  speeds,  set  up  an  internal  electric  field  which 
speeds  up  the  slower  moving  species  (usually  the  atomic 
species)  and  slows  down  the  faster  moving  species  (usually  the 
electrons  or  holes).  However,  in  the  following  treatments,  the 
electric  field  term  will  be  assumed  negligible  (based  on  the 
results  of  Muller  et  al.  (78)). 

Weisberg  and  Blanc  (84)  numerically  calculated  the 
solution  for  Case  I  with  charged  species  (from  equations  5-24 
and  5-25).  They  did  this  for  three  sub-cases,  in  which  the 
charge  difference  between  the  interstitial  and  substitutional 
species  is  1,  2,  and  3,  respectively  (equation  5-25  is  for  the 
second  of  those  three  sub-cases;  similar  equations  are  easily 
derived  for  the  other  two).  These  were  compared  to  the 
uncharged  case.  The  results  are  shown  in  figure  5-4.  The 
profiles  for  the  charged  species  drop  off  more  quickly  than  the 


uncharged  (error  function)  case  due  to  the  concentration 
dependence  of  Deff.  Another  way  to  look  at  the  effect  is  to 
consider  the  conversion  reaction:  Ij_+  +  V  =  Is“  +  2h+. 

Near  the  surface,  where  the  substitutional  level  is  high,  the 
hole  concentration  is  high  (since  for  doping  >  nj_,  the  source 
of  holes  is  the  acceptor  atoms).  Therefore,  the  equilibrium  is 
shifted  right  to  left  (relative  to  low  concentration).  Since 
the  effective  diffusion  coefficient  is  approximately 
proportional  to  the  Cj_/Cs  ratio,  Deff  is  higher 
here.  Conversely,  for  lower  substitutional  concentrations, 
(i.e.,  at  the  tail  of  the  profile),  the  equilibrium  is  shifted 
to  the  right  relative  to  high  concentrations,  and  Deff  is 
less  there.  This  same  reasoning  applies  whenever  the  hole  (or 
electron)  concentration  is  changed,  either  locally  or  overall, 
resulting  in  a  higher  or  lower  Deff*  Case  II  diffusion  with 
charged  species  would  probably  also  result  in  concentration 
dependent  diffusion  -  shifting  equation  5-20  to  the  right  or 
left  due  to  a  change  in  the  electron  or  hole  concentration  - 
probably  also  changes  the  ratio  CV'/(CV'  +  Cs '  )  ,  and 
hence  changes  DS]_ow  for  different  concentrations. 

It  has  been  suggested  by  Ghandhi  (79)  that  a 
concentration-dependent  diffusion  coefficient  might  also  result 
for  uncharged  species  if  there  is  a  problem  in  achieving 
equilibrium  between  C^,  Cs,  and  Cv  locally.  It  was 
assumed  in  the  present  work  that  achieving  this  equilibrium 
is  no  problem. 

Tuck  utilized  Weisberg  and  Blanc's  solutions  to  solve  Case 
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IV  (k]_>0  and  large)  for  charged  species.  His  theoretical 
results  for  Zn  in  GaAs ,  from  (73),  are  shown  in  figure  5-5.  H 
claims  that  they  fit  the  qualitative  shape  of  the  actual 
profiles  in  the  bulk. 

To  estabablish  the  dependence  of  the  diffusion  on  the 
component  vapor  pressures,  we  first  write  equation  5-25  as: 


(5-26) 


Then  for  the  case  of  Cr  diffusion  in  GaAs,  in  which  the 
impurity  is  on  a  Ga  site,  the  value  of  Cv'  can  be  related  to 
the  arsenic  vapor  pressure  by  the  defect  reaction:  GaAs  +  As2 
=  GaAs  +  2Vq3;  leading  to: 
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where  X6  is  the  equilibrium  constant  between  Vq3  and 
pAs2 *  (According  to  Arthur  (85),  Asz  is  the  dominant 
arsenic  species  in  the  vapor  phase  on  the  gallium-rich  side; 
As4  is  dominant  on  the  arsenic  side).  Equation  5-27  shows  the 
dependence  of  Deff  on  the  arsenic  vapor  pressure. 

The  above  treatment  regarding  vapor  pressure  dependence 
was  for  the  charged  diffusing  species  case;  similar  treatment 
of  the  uncharged  case  (from  equation  5-8a,  rather  than  from 
5-26),  and  assuming  Cj_'<<Cs',  leads  to: 
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Figure  5-5. 


Theoretical  Zn  in  GaAs  diffusion  profiles,  f 
ref .  73 . 
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Specific  systems  can  now  be  classified  according  to  their 
respective  categories.  For  example,  the  diffusion  of  Cu  in  Ge, 
as  described  before,  would  be  Case  III,  with  the  erfc  diffusion 
(Case  II,  or  vacancy  diffusion  controlled)  at  the  surface  and 
the  uniform  levels  deep  in  the  bulk  building  up  with  time.  For 
shorter  times,  however,  it  was  observed  (56)  that  the  bulk 
profile  did  drop  off  slightly,  presumably  because  of 
insufficient  time  for  the  interstitials  to  diffuse  all  the  way 
in  prior  to  the  initiation  of  interstitial-to-substitutional 
conversion.  This  is,  therefore.  Case  IV.  Van  der  Maesen  and 
Brenkman's  (86)  profiles  of  Cu  diffusion  in  Ge  (as  shown  in 
Tuck  (7))  clearly  show  this  behavior.  Tweet  and  Gallagher  (57) 
found  that  the  rate  of  increase  in  the  bulk  level  depends  on 
the  number  of  dislocations  in  the  crystal  -  the  more 
dislocations,  the  faster  the  increase.  This  is  consistent  with 
the  model,  since  k;  would  increase  with  increasing  vacancy 
sources,  leading  to  a  faster  rise. 

The  diffusion  of  Ag  and  Cu  in  InP  would  be  classified  the 
same  way  as  Cu  in  Ge,  and  Sn,  Mg,  Cd,  and  Zn,  at  low 
concentrations,  would  be  vacancy  diffusion  controlled,  or  Case 
II  (77).  At  high  concentrations,  these  would  be  classified  in 
Cases  I,  III,  or  IV;  the  higher  concentrations  cause  a  higher 
vacancy  generation  rate.  The  diffusion  of  Zn  in  GaAs  has  been 
studied  extensively,  as  previously  mentioned.  Many  have 
assumed  Case  I  (68,  71,  72,  for  example).  Weisberg  and  Blanc's 


(84)  theoretical  fitting  of  Case  I  (vacancy  equilibrium), 
charged  species,  to  the  data  of  Cunnell  and  Gooch  (72)  for 
various  Zn  surface  concentrations  is  very  good,  except  for  a 
small  inflection  or  dip  in  from  the  surface.  Tuck  and  Kadhim 
(73),  who  observe  a  more  pronounced  inflection  in  their 
profiles  (which  are  for  diffusions  done  for  various  times), 
propose  that  it  really  follows  Case  IV,  and  associate  the 
inflection  with  the  switch  of  the  profile  from  the  vacancy 
equilibrium  region  to  the  vacancy  generation  region.  They  show 
that  the  build-up  with  time  is  consistent  with  this.  Shaw  and 
Showan  (69)  claim  that  the  idea  of  "uniform  vacancy  generation 
in  the  bulk"  is  inappropriate  for  Zn  in  GaAs  since  the  levels 
are  too  high  for  normal  vacancy  generation.  They  say  that 
diffusion-induced  dislocations  cause  the  regeneration  of 
vacancies  in  the  bulk.  (This  would  explain  why  at  a  low 
concentration  of  Zn,  vacancy  diffusion  would  control  Zn 
diffusion.)  Chiaretti  and  Cognetti  (75)  claim  that  a  certain 
amount  of  time  (between  1  and  2  hours  at  800 °C,  for  example) 
must  elapse  before  "vacancy  equilibrium  is  maintained"  and  deep 
diffusion  at  relatively  high  concentrations  is  achieved  (i.e., 
Case  III  or  IV).  Others  (87,  88)  have  made  similar  statements. 

It  is  clear  from  this  discussion  and  from  the  model  itself 
that  the  initial  state  of  the  host  material,  especially  the 
dislocation  density  and  the  initial  vacancy  concentration,  is 
very  important  in  establishing  which  case  is  dominant  and 
therefore  the  corresponding  diffusion  behavior.  If  there  are 
many  vacancy  sources  present,  or  if  the  initial  vacancy 
concentration  is  high  (which  depends  on  the  history  of  the 
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material),  the  diffusion  might  be  Case  I  at  the  outset. 
Otherwise,  it  may  take  some  time  before  that  occurs  (either 
Case  III  or  IV).  Luckily,  however,  both  these  cases  ultimately 
become  Case  I;  therefore  the  Dfast  values  can  be  measured 
accurately  no  matter  what  the  starting  conditions  are,  as  long 
as  enough  time  is  allowed  to  reach  vacancy  equilibrium. 

5 . 6  Interstitial  to  Substitutional  Conversion 

In  the  preceding  treatment  of  the  diffusion  mechanism,  the 
interstitial  to  substitutional  conversion  was  assumed  to  take 
place  via  the  reaction: 

Ii  +  V  =  Is  (5-29) 

This  is  commonly  called  the  Frank-Turnbull  (F-T)  mechanism,  or 
sometimes  the  vacancy  mechanism.  (This  is  not  to  be  confused 
with  substitutional  dopant  diffusion,  which  is  sometimes 
referred  to. as  "diffusion  via  vacancies"  or  "a  vacancy  model 
for  diffusion".  Because  of  this  confusion,  the  above 
conversion  mechanism  will  be  referred  to  as  the  "F-T 
mechanism",  which,  like  the  conversion  mechanism  to  be 
discussed  below,  is  part  of  an  interstitial  diffuson 
mechanism. )  It  was  assumed  that  a  shortage  of  vacancies  during 
the  diffusion  suppresses  the  bulk  concentration  levels 
(resulting  in  Case  III  or  IV  rather  than  Case  I).  However, 
some  researchers  have  recently  proposed  that  in  some 
impurity /diffusion  host  systems,  the  interstitial  atom  does  not 
wait  for  a  vacancy  to  be  formed,  but  instead  creates  its  own 
space  by  "kicking  out"  a  lattice  atom,  which  then  becomes  an 
interstitial  (this  is  actually  an  adaptation  of  the 
interstitialcy  model  of  diffusion,  this  time  applied  to  an 
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impurity  in  a  host  material).  This  is  commonly  called  the 
"kickout  mechanism."  The  conversion  reaction  in  this  case  is: 

Ii  +  =  Is  +  Li  (5-30) 

where  I  and  L  denote  impurity  atom  and  lattice  atom, 
respectively,  i  denotes  interstitial  site,  l  denotes  lattice 
site,  and  s  denotes  substitutional  site.  (Note:  a  lattice 
site  for  a  lattice  atom  is  the  same  position  as  a 
substitutional  site  for  an  impurity  atom).  This  has  recently 
been  applied  to  the  case  of  gold  in  silicon  by  GOsele  et  al. 
(89)  and  Hill  et  al.  (65).  Schaake  et  al.  (66,  90)  have 
applied  this  mechanism  to  explain  anomalous  diffusion  of 
phosphorus  in  silicon,  and  impurity  getter ing  in  mercury 
cadmium  telluride. 

For  long  times,  high  temperatures,  and/or  high  dislocation 
densities,  the  results  are  exactly  the  same  in  both  cases:  the 
impurity  diffuses  with  an  effective  diffusion  coefficient  equal 
to  Dj_Cj_'/(Ci'  +  Cs '  ) .  However,  for  short  times,  low 
temperatures,  and/or  low  dislocation  densities,  the  diffusion 
is  not  limited  (the  bulk  level  is  not  suppressed)  by  vacancy 
diffusion  or  generation,  but  by  the  kickout  reaction  (usually 
considered  to  be  insignificant)  and  by  removing  of  the  excess 
interstitial  lattice  atoms,  either  by  outdif fusion,  or  by 
internal  interstitial  sinks  (analogous  to  internal  vacancy 
sources ) . 

The  governing  equations  of  the  kickout  mechanism  are  very 
similar  to  those  of  the  F-T  mechanism.  The  surface  vacancy 
indiffusion  equations  are  replaced  by  interstitial  lattice  atom 
outdiffusion  equations,  and  the  bulk  vacancy  generation  term  is 
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replaced  by  a  combination  of  the  kickout  reaction  term,  the 
interstitial  lattice  atom  outdiffusion  term  (sometimes 
including  a  surface  evaporation  term),  and  the  interstitial 
lattice  atom  annihilation  term  (due  to  sinks  -  perhaps 
dislocations).  In  Hill's  (65)  analysis  of  Au  in  Si,  the  bulk 
rise  in  the  impurity  level  is  explained  by  just  the 
interstitial  lattice  atom  outdiffusion  (from  a  finite  crystal 
so  that  the  bulk  level  of  excess  interstitial  lattice  atoms 
decreases  with  time)  and  the  surface  evaporation  term.  In  the 
F-T  model,  the  bulk  rise  of  impurity  atoms  cannot  be  reasonably 
explained  by  vacancies  diffusing  into  a  finite  crystal  because 
of  the  smaller  number  and  slower  diffusivity  of  vacancies 
(compared  to  interstitial  lattice  atoms).  However,  Hill  must 
still  use  an  interstitial  sink  to  explain  Au  behavior  in 
previously  unannealed  Si.  And  Gosele  (39)  uses  the 
interstitial  sink  term  in  his  Au  in  Si  analysis. 

Because  the  governing  equations  and  reactions  are  so 
similar,  both  models  predict  similar  behavior  for  the  different 
conditions  (i.e.,  high  or  low  dislocation  densities, 
temperatures,  etc.),  and  it  is  very  difficult  to  distinguish 
the  two  in  many  cases.  However,  some  (65,  89)  propose  the 
following  distinguishing  differences  between  the  two  models: 

(1)  the  bulk  impurity  level  build-up  initially  follows  a  /t 
dependence  (whether  the  bulk  build-up  is  controlled  by 
interstitial  outdiffusion  or  annihilation),  as  opposed  to  a  (1 
-  exp(-t/0))  dependence  for  the  F-T  mechanism.  However,  as  the 
impurity  level  approaches  the  solubility  level,  the  kickout 
model  also  follows  an  exponential  approach;  (2)  the  surface 


diffusion,  which  dominates  for  short  times/low  dislocation 
densities,  is  theoretically  governed  by  a  concentration 
dependent  diffusion  coefficient  proportional  to  Cs“2  in  the 
kickout  case,  as  opposed  to  Deg£  =  DVCV’/(CV'  +  Cs ' ) , 
which  is  independent  of  concentration,  in  the  F-T  mechanism. 

The  first  difference,  which  was  used  to  distinguish  the  two 
models  in  the  Au  in  Si  case,  is  not  always  easily 
distinguishable.  An  initial  /t  dependence  which  later  becomes 
an  expontential  dependence  is  very  similar  to  a  straight 
exponential,  especially  when  the  concentrations  are  initially 
low  and  there  is  some  error  in  the  experimental  measuring 
technique.  Even  then,  Huntley  et  al.  (91)  propose  that  the  /t 
dependence  can  be  explained  by  the  F-T  (vacancy  generation) 
model.  The  second  difference  has  not  been  proven.  For  the 
cases  that  the  kickout  model  has  been  proposed,  the  impurity 
profiles  at  the  surface  were  not  accurately  measured,  or  other 
surface  effects  obscured  them.  Or  in  the  case  of  P  in  Si, 
neither  model's  fit  to  the  experimental  points  was  good. 

Another  difficulty  in  distinguishing  the  two  mechanisms 
is  that  virtually  all  phenomena  explained  by  one  model  can  just 
as  easily  be  explained  by  the  other.  These  phenomena  include 
gettering,  enhanced  diffusion,  and  effects  of  changing  the 
vapor  pressure,  dislocatiom  density,  and/or  vacancy 
concentration.  The  reason  for  this  is  that  an  equilibrium 
exists  between  vacancies  and  interstitial  lattice  atoms  (just 
as  one  exists  between  vacancies  and  interstitial  impurities), 
with  the  following  reaction: 


Lj_  +  V  = 


(5-31) 


In  fact,  the  reaction  for  the  kickout  mechanisms  can  be 
converted  to  the  vacancy  mechanism  by  just  adding  this  reaction 
to  it.  An  example  of  how  either  model  can  explain  a  phenomenon 
is  as  follows:  increasing  the  initial  vacancy  concentration 
(by  changing  the  vapor  pressure,  for  instance)  would  mean  that 
less  vacancies  would  have  to  be  generated  in  the  F-T  model, 
leading  to  a  faster  build-up  of  impurities.  Or  using  the 
kickout  model,  increasing  the  vacancy  concentration  shifts 
reaction  5-31  to  the  right,  resulting  in  fewer  interstitial 
lattice  atoms.  Therefore,  it  is  easier  to  remove  the  excess 
interstitial  lattice  atoms,  allowing  a  faster  increase  in  the 
impurity  atom  concentration  level.  Similar  arguments  can  be 
made  for  the  other  phenomena. 

Because  of  these  considerations,  it  is  very  difficult  to 
distinguish  the  two  models  in  many  cases.  As  Gosele  et  al. 

(39)  point  out,  while  some  systems,  such  as  Au  in  Si,  can  be 
explained  better  by  the  kickout  mechanism,  and  some,  such  as  Cu 
in  Ge ,  by  the  F-T  mechanism,  there  is  a  large  body  of 
experimental  results  which  are  compatible  with  both  models  and 
are  impossible  to  distinguish.  Furthermore,  some  systems  might 
operate  in  one  manner  under  some  conditions  and  in  the  other 
manner  under  other  conditions. 

There  are  still  some  problems  with  the  kickout  model.  The 
initial  reason  for  proposing  the  kickout  mechanism  is  the 
failure  af  the  F-T  model  to  explain  such  features  as 
anomalously  fast:  (1)  bulk  vacancy  generation;  and  2)  vacancy 
diffusion.  (Some  investigators  (69)  have  explained  these 
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anomalies  with  the  proposal  that  high  levels  of  impurities 
introduce  dislocations,  which  in  turn  generate  vacancies,  and 
also  enhance  vacancy  diffusion.)  However,  certain  phenomena 
and  anomalies  in  the  kickout  mechanism  are  just  as  difficult  to 
explain.  For  instance,  the  interstitial  annihilation  process, 
which  is  needed  to  explain  some  phenomena,  is  not  understood, 
and  the  equations  regarding  the  process  are  not  known.  Also, 
an  extremely  wide  range  of  Si  interstitial  diffusion 
coefficients  are  needed  to  model  impurity  diffusion  in  silicon 
(10--*-7  and  10-7  cm^/sec  for  P  and  Au  respectively). 

Also,  as  mentioned  earlier,  the  fit  to  experimental  data  is  not 
extremely  good.  As  one  can  see,  there  are  still  questions 
regarding  both  models. 

Finally,  if  one  assumes  that  the  vacancy  diffuses  in  an 
interstitial  manner  (suggested  by  Schaake  (92))  then  the  two 
models  are  really  one  in  the  same.  Near  the  surface,  vacancies 
diffuse  in  by  lattice  interstitials  diffusing  out  (which  are 
created  by  forming  a  vacancy).  The  bulk  level  is  limited  by 
vacancy  generation  which  is  caused  by  lattice  atoms 
transforming  to  interstitial  atoms  and  being  annihilated  and/or 
diffusing  out.  This  is,  in  essence,  how  Zahre  et  al.  (93) 
explain  the  anomalous  diffusion  behavior  of  S  in  GaAs  (the  same 
type  of  behavior  seen  with  P  in  Si). 

To  emphasize  how  similar  the  two  models  are  in  regards  to 
their  governing  equations  and  reactions,  and  since  they  are,  in 
fact,  indistinguishable  in  many  cases,  the  vacancy  regeneration 
parameter,  0,  will  from  this  point  on  be  called  the 
"equilibrium  time  constant".  This  parameter,  with  the  same 


values  which  were  determined  in  this  research,  can  then  be  used 
with  either  model,  or  any  other,  if  necessary.  It  will  be  used 
in  describing  the  final  attainment  of  either  vacancy 
equilibrium  or  interstitial  lattice  atom  equilibrium,  or  any 
other  phenomenon  that  limits  the  interstitial-to-substitutional 
conversion  in  an  exponential  manner.  For  instance,  just 
"finding "  the  vacancies,  even  if  there  are  enough  or  more  than 
enough,  might  limit  the  conversion,  and  would  probably  have  an 
exponential  time  dependence  since  the  rate  of  finding  the 
vacancies  is  probably  proportional  to  how  many  are  left.  The 
equilibrium  time  constant  parameter,  and  the  same  values  for 
it,  could  be  used  in  this  case,  and  equations  virtually 
identical  to  the  ones  described  earlier  would  be  used  to  model 
the  diffusion  process. 

In  conclusion,  the  differences  between  the  two 
interstitial-to-substitutional  conversion  models  (and  perhaps 
even  other  conversion  models)  are  extremely  subtle,  and  in 
fact,  may  be  the  same  anyway.  The  important  thing  is  in 
identifying  which  general  diffusion  mechanism,  for  example, 
interstitial,  substitutional,  or  substitutional/interstitial 
dissociative,  is  operative  for  a  particular  system  and/or 
situation.  Whether  the  substitutional/interstitial 
dissociative  mechanism  operates  via  the  kickout  mechanism  or 
the  F-T  mechanism,  or  some  other,  is  not  terribly  important 
since  both  predict  virtually  the  same  behavior. 

5 . 7  Dn  and  Ea  Values 

The  diffusion  coefficient  can  be  written  in  the  form: 
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D  =  D0  exp  (-Ea/kT)  (94),  and  DQ  and  Ea  values  can  be 
easily  obtained  from  In  D  vs.  1/T  plots.  There  are 
correlations  between  the  particular  diffusion  mechanism  and 
reported  D0  and  Ea  values,  thus,  these  quantities  may  be 
used  to  identify  the  mechanism  for  a  given  system. 

For  purely  interstitial  diffusion  in  compound 
semiconductors,  the  Ea  values  are  usually  low  (e.g.,  0.52  eV 
for  Cu  in  InAs  (95))  with  D0  values  in  the  IQ-3  cm2 /sec. 
range  (3.6  x  10“3  cm2/sec  for  Cu  in  InAs).  This  results  in 
high  D  values  (1.3  x  10-3  cm2/sec  for  Cu  in  InAs  at 
800 °C).  At  the  other  extreme,  pure  substitutional  diffusion 
has  high  Ea  values  (4.16  eV  for  Se  in  GaAs  (96)),  and  fairly 
high  D0  values  (3  x  103  cm2/sec  for  Se  in  GaAs), 
resulting  in  very  low  D  values  (1  x  10-^6  cm2/sec  for  Se  in 
GaAs  at  800 °C).  The  value  of  Ea  is  high  for  substitutional 
diffusion  because  there  are  two  parts:  the  energy  for  atom 
motion  and  the  energy  for  vacancy  formation,  the  latter 
depending  on  bond  strength  and  therefore  quite  large.  On  the 
other  hand,  Ea  for  interstitial  diffusion  depends  only  on 
atom  motion  (squeezing  through  the  other  atoms),  and  is  lower. 

For  atoms  which  are  presumed  to  diffuse  by  the 
substitutional/interstitial  dissociative  process,  the  D  values 
usually  lie  between  those  of  purely  substitutional  diffusers 
and  purely  interstitial  diffusers  and  most  have  Ea  values  in 
the  0.5  to  2.0  eV  range.  The  D0  values  are  lower  than  the 
interstitial  diffusers  (in  the  10-4  to  10-7  cm2/sec. 
range)  resulting  in  D  values  usually  less  than  interstitial, 
but  higher  than  substitutional  diffusers.  The  D  values  are 
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mechanism.  For  vacancy  diffusion  controled  impurity  diffusion 
(Case  II),  the  Ea  values  are  relatively  higher  (-1  to  2  eV), 
with  corresponding  D0  values  giving  net  D  values  at  the  lower 
end  of  the  sub./int.  D  range.  On  the  other  hand,  in  those  case 
controlled  by  vacancy  generation  and  interstitial  diffusion 
(Cases  I,  III,  and  IV),  the  Ea  values  are  usually  less  than  1 
eV  and  with  corresponding  0o  values  that  give  net  D  values  at 
the  higher  end  of  the  range. 

Of  course,  all  of  this  is  reasonable  when  one  looks  at  the 


equations  for  the  effective  diffusion  coefficients.  Dfast  = 

Di  CiMCi'  Cs');  therefore,  the  activation  energy 
should  be  close  to  the  activation  energy  for  interstitial 
diffusion,  modified  slightly  by  the  temperature  dependences  of 
the  equilibrium  concentrations  of  the  interstitials  and 
substitutionals .  The  D0  value  would  be  less  because 
C i * / ( C i *  +  Cs ' )  is  usually  much  less  than  unity. 

Similarly,  the  Ea  value  for  DS]_ow  (which  equals  Dv 
CV'/(CV'  +  Cs '  )  )  would  be  close  to  the  vacancy  diffusion 
value,  which  according  to  (83)  is  2.1  eV  for  VQa  in  GaAs. 

An  example  of  this  second  correlation  is  the  case  of  Cu  in 


InSb.  For  substitutional/interstitial  diffusion  believed  to  be 


vacancy  diffusion  controlled,  Stocker  (97)  reported  DQ  and 
Ea  values  of  9  x  lO'^cm^/sec.  and  1.08  eV  respectively, 
leading  to  a  D  value  of  5  x  10-^  cm^/sec.  at  480°C.  In 
addition,  Stocker,  and  also  Boltaks  and  Sokolov  (98),  observed 
deeper  diffusion  branches,  with  DQ  and  Ea  values  of  3  x 
10~5  cm^/sec.  and  0.37  eV  respectively,  resulting  in  a  D 
value  of  1  x  10""7  cm^/sec.  at  480°C.  This  was  attributed 
to  the  interstitial  diffusion  controlled  case. 

5 . 8  Summary 

This  chapter  reviewed  the  substitutional/interstitial 
dissociative  diffusion  mechanism,  which  is  being  used  with 
increasing  frequency  to  describe  impurity  diffusion  in 
semiconductors.  The  governing  reactions  and  equations,  for 
different  rate  limiting  cases  depending  on  which  features  of 
the  process  dominate,  have  been  given.  In  this  treatment  of 
the  model,  the  Frank-Turnbull ,  or  "vacancy",  model  for  the 
interstitial-to-substitutional  conversion  has  been  assumed. 
Specific  systems  which  are  presumed  to  follow  the 
substitutional/interstitial  model  have  been  considered  to 
establish  which  of  the  different  rate  limiting  cases  they  fit. 
A  new  mechanism  for  the  conversion  process,  the  so-called 
"kickout"  mechanism,  has  been  discussed.  It  has  been  found 
that  in  some  cases,  this  model  fits  the  data  better,  and  in 
other  cases,  the  Frank-Turnbull  model  fits  the  data  better; 
however,  in  most  cases,  the  two  are  indistinguishable.  It  was 
shown  how  D0  and  Ea  values  for  the  effective  diffusion 
coefficients  could  be  used  to  establish  if  the 
substitutional/interstitial  diffusion  mechanism  is  operative 


for  a  specific  system  and  if  so,  establish  which  is  the  rate 


limiting  case, 


6.  CHROMIUM  DIFFUSION  IN  GALLIUM  ARSENIDE 


6 . i  Introduction 

Using  the  same  principles  concerning  ternary  phase 
equilibria  that  were  used  in  the  solubility  studies, 
e] ’.-defined  chromium  diffusion  experiments  were  performed; 
this  chapter  reports  on  those  studies.  First,  a  literature 
review  of  previous  and  current  work  on  chromium  redistribution 
in  gallium  arsenide  is  presented,  emphasizing  the  following 
major  points:  1)  there  have  been  several  recent  studies  on 
this  subject;  however,  these  studies  were  made  without 
knowledge  of,  and  usually  without  concern  for,  phase 
equilibria;  and  2)  apparently,  there  is  a  complex  diffusion 
mechanism,  more  complicated  than  simple  interstitial  or 
substitutional  diffusion,  or  even  a  simple  combination  of  the 
two.  Next,  the  approach  for  this  work  concerning  Cr  diffusion 
is  outlined  and  discussed,  followed  by  the  experimental 
procedures.  The  results  are  then  given,  which  are  shown  to  fit 
the  substitutional/interstitial  dissociative  diffusion 
mechanism,  and  values  for  key  parameters  pertaining  to  that 
model  are  given.  In  addition  it  is  shown  how  the  defect 
chemistry  (especially  concerning  vacancies  and  impurities)  of 
gallium  arsenide  can  be  studied  using  chromium  as  a  "probe". 
Finally,  in  this  chaper  and  the  next  (which  deals  with  several 
surface  phenomena),  some  of  the  results  of  others  are  explained 
in  terms  of  the  results  of  this  work. 


6.2  Literature  Review  of  Chromium  Redistribution 


depositing  each  element  on  a  GaAs  surface  and  annealing  the 
sample  under  an  arsenic  over-pressure.  He  calculated  the 
diffusion  coefficients  by  measuring  the  junction  depths  (by  the 
stain-etching  method)  and  assuming  an  error  function  diffusion 
profile.  cor  Cr  at  900°C,  he  reported  a  D  of  3  x  lO-^ 
cm^/sec.  However,  his  results  may  be  in  doubt  since  the  Cr 
concentrations  he  reports  -  derived  from  the  junction  depths  - 
are  much  higher  than  the  maximum  Cr  solubility  ((27)  and  from 
the  present  work),  in  addition  to  the  fact  that  elemental  Cr  is 
not  a  proper  diffusion  choice  since  Cr  and  GaAs  are  not  in 
equilibrium  with  each  other. 

Sato  (34)  in  1973  first  speculated  that  surface  conversion 
in  n-type  GaAs  arose  from  chromium  diffusion  out  of  the  GaAs 
(where  the  chromium  presumably  compensated  the  n-type 
impurities)  and  into  the  Si02  cap  (which  is  used  in  GaAs 
processing  to  prevent  arsenic  vaporization  during  high 
temperature  anneals).  Sato  calculated  a  diffusion  coefficient 
at  850°C  from  an  experiment  virtually  identical  to  that  of 
Khludkov,  with  a  reported  D  value  of  1.5  x  10-^  cm^/sec. 
However,  the  Cr  concentrations  were  again  much  higher  than  the 
maximum  solubility  level,  and  elemental  Cr  was  used  as  the 
source.  Sato  also  studied  the  diffusion  of  Cr  from  Cr-doped 
GaAs  into  the  oxide  layer.  He  calculated  diffusion 
coefficients  by  performing  electrical  measurements  and  fitting 
a  composite  medium  diffusion  solution  to  the  results,  again 
assuming  a  simple  error  function  profile  and  that  Cr  acts  as  a 
singly  ionized  acceptor.  The  resulting  D  values  were  in  the 
10**14  to  10”12  range  for  T  =  700°C  to  800°C,  with  an 


activation  energy  of  3.6  eV. 


Little  work  was  done  on  this  subject  until  1977  when 
several  researchers  (100,  101)  reported  problems  with  GaAs 
devices  that  they  attributed  to  a  depletion  of  chromium  near 
the  surface  after  heat  treatments.  Deveaud  (48),  that  same 
year,  observed  "abnormal"  Cr  diffusion  toward  the  surface  - 
depletion  near  the  surface  and  buildup  right  at  the  surface  - 
after  it  was  implanted  in  GaAs  and  annealed.  In  1978,  Partin 
(102)  estimated  diffusion  coefficients  from  hole  lifetime 
measurements  that  were  two  to  three  orders  of  magnitude  higher 
than  the  values  previously  reported  for  the  same  temperatures. 
This  was  the  first  indication  that  a  deep  diffusion  process 
might  be  occurring. 

From  1979  to  the  present,  much  work  has  been  done  studying 
the  "unusual"  behavior  of  chromium  m  gallium  arsenide.  Many 
have  reported  on  Cr  gettering  to,  or  building  up  at,  regions  of 
lattice  damage  or  stress,  either  regions  damaged  by  ion 
implantation,  the  damaged  and/or  stress  surface  region,  or  just 
the  free  surface  (5,  39,  40,  42,  103-120).  Many  also  report  on 
the  depletion  of  chromium  near  the  surface,  which  may  be  a 
result  of  the  gettering  to  the  surface,  or  simple  outdif fusion 
(and  vaporization),  or  both  (41,  42,  50,  78,  105,  107-112,  118, 
121,  122).  However,  there  is  disagreement  regarding  the  nature 
and  exact  cause  or  causes  of  the  redistribution,  which  will  now 
be  discussed. 

The  most  often  reported  type  of  concentration  profile  of 
Cr  following  an  anneal  of  Cr-doped  GaAs,  is  one  with  a  surface 
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build-up  of  Cr,  as  high  as  three  orders  of  magnitude  above  the 
bulk  level,  over  the  first  few  tenths  of  microns  from  the 
surface  of  the  GaAs,  and  then  a  Cr  depletion  over  the  next  two 
or  three  microns.  (In  most  of  these  cases,  the  Cr  profile  was 
obtained  using  SIMS.)  Favennec  et  al.  (105),  for  example,  show 
a  clear  example  of  this  behavior.  Many  investigators  believe 
this  is  due  to  gettering  of  chromium  to  the  encapsulent,  which 
is  usually  Si02  or  Si3N,*  (40-42,  103  ,  105,  112  ,  123  ).  Huber  et 
al.  (42)  and  Evans  et  al .  (103)  proposed  that  the  Cr 
preferentially  segregates  to  the  strained  mismatched  region 
near  the  cap/GaAs-crystal  interface,  implying  a  lowering  of  the 
Cr  activity  there.  This  is  consistent  with  the  observation  of 
Magee,  et.  al .  (123-125)  that  Cr  is  preferentially  segregated 
to  regions  of  high  dislocation  density  that  could  arise  from 
residual  stress.  However,  Magee  (125)  suggests  that  the  actual 
gettering  may  be  due  to  stress-induced  diffusion  (a  field 
effect).  The  "cap-gettering "  theory,  for  whatever  reason,  is 
supported  by  the  results  of  Kasahara  et  al.  (112),  Feng  et  al. 
(41),  and  Eu  et  al.  (40),  who  report  that  when  caps  are  not 
used,  but  with  an  AsH3/H2  flow  in  the  ambient  to  prevent  loss 
of  As,  the  build-up  is  reduced.  In  fact,  Kasahara  (112) 
reported  that  the  build-up  was  completely  eliminated.  Eu  (40) 
suggested  that  the  Cr  build-up  and  depletion  associated  with 
the  use  of  caps  may  occur  because  the  bulk  chromium  levels  are 
above  the  solubility  limit,  and  Cr  precipitates  nucleate  at  the 
surface  dislocations.  He  proposed  an  alternate  explanation, 
that  the  redistribution  may  be  due  to  arsenic  diffusing  through 
the  caps,  leaving  a  Ga-rich  phase  which  getters  Cr. 
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Udagawa  (108),  who  observed  a  high  Cr  surface  build-up  at 
high  temperatures  when  a  cap  was  not  used,  suggested  that  the 
Cr  diffuses  into  the  Ga  pools  or  droplets,  which  form  when  the 
GaAs  is  annealed  (at  T  >  800°C  in  a  hydrogen  atmosphere  in  his 
case)  and  "which  have  a  high  Cr  solubility."  At  T  <  8Q0°C,  no 
Ga  droplets  are  observed,  and  no  build-up  is  reported. 

Vasudev  et  al.  (109)  proposed  that  the  chromium  was  not 
gettered  to  the  encapsulent;  the  build-up  was  worse  when  a  cap 
was  not  used.  They  proposed  that  the  build-up  and  depletion  is 
controlled  by  a  balance  between  the  surface  solubility  limit  (a 
fixed  quantity)  and  the  original  Cr  density  in  the  substrate. 
Tuck  et  al.  (120)  reported  even  stranger  behavior  after 
annealing  Cr-doped  GaAs  in  an  evacuated  closed  quartz  ampoule 
(with  and  without  elemental  arsenic  in  the  ampoule);  they 
observed  a  large  chromium  build-up  in  the  first  40  microns  from 
the  surface,  and  a  depletion  throughout  the  rest  of  the  water. 
They  attributed  this  to  a  sudden  shift  in  the  vacancy/chromium 
equilibrium  after  an  initial  vacancy  undersaturation. 

Several  workers  have  attemped  to  model  the  diffusion 
behavior.  Most  have  assumed  simple  error  function  diffusion 
profiles,  with  or  without  surface  evaporation  terms,  and  fit 
these  to  measured  profiles  after  annealing  uniformly  Cr-doped 
GaAs,  or  Cr  implanted  GaAs,  or  Cr-doped/undoped  GaAs  junctions 
(39,  41,  50,  82,  105,  108,  110,  112,  122,  126).  Most  of  the 
reported  D  values  -  determined  from  the  fitting  -  are  in  the 
lO--^  to  10"^0  cm^/sec  range  from  700  to  1000°C.  Both 
Udagawa  et  al.  (108)  and  Kasahara  et  al.  (112)  find  that  while 
the  diffusion  coefficient  was  constant  with  time  for  short 
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annealing  times  (0  to  45  minutes  ?  900°C,  for  example),  it 
became  much  larger  at  longer  times. 

Just  recently,  Morkoc  et  al.  (37)  and  Palmateer  et  al. 
(127),  who  studied  Cr  redistribution  during  MBE  GaAs  growth  and 
processing,  claimed  that  the  fast  redistribution  of  Cr  observed 
in  their  work  cannot  be  explained  by  a  normal  diffusion  model, 
and  that  interstitial  Cr  and/or  Cr  segregation  may  be  the 
driving  force. 

Several  have  also  studied  the  effects  of  changing  the 
arsenic  overpressure  in  the  chromium  diffusion  (5,  40,  108, 

112,  122).  Kasahara  et  al.  (112)  and  Udagawa  et  al.  (108) 
reported  that  increasing  the  arsenic  pressure  apparently 
increased  the  Cr  diffusion;  more  Cr  depletion  was  observed 
after  outdiffusion  when  the  Pas2  was  increased.  But 
Kasahara  (112)  and  Asbeck  (82)  both  suggest  that  it  is  the 
surface  concentration  that  is  changing,  and  not  the  diffusion 
coefficient.  (If  one  fits  erf  diffusion  profiles  to  Kasahara's 
results,  but  with  a  constant  D  and  a  varying  surface 
concentration,  one  sees  that  they  fit  virtually  perfectly.) 

Tuck  et  al.  (5,  120)  have  been  the  only  ones  to  study  the 
indiffusion  of  Cr  into  GaAs  since  Khludkov  (99)  and  Sato  (34). 
They  did  this  by  annealing  GaAs  and  a  piece  of  radioactive 
Cr-^  in  an  evacuated  closed  quartz  ampoule.  (They  estimated 
that  the  GaAs  initially  contained  -10^-®  cm”3  Cr^2 , 
perhaps  even  more.)  Using  an  etch  and  count  method,  they 
obtained  Cr  profiles  for  various  temperatures  and  times.  The 
profiles  were  much  like  Case  III  and  Case  IV  diffusion  as 
discussed  in  chapter  5,  and  in  their  articles  they  suggested 
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that  the  substitutional/interstitial  dissociative  diffusion 
mechanism  was  operative.  In  addition  they  found  that 
increasing  the  arsenic  pressure,  by  adding  elemental  arsenic, 
decreased  the  amount  of  Cr  that  diffused  in. 

The  fact  that  for  short  annealing  times  and  low 
temperatures  (<1Q00°C)  -  the  conditions  at  which  most  of  the 
other  workers  did  their  experiments  -  this  complex  mechanism 
reduces  to  a  simple  error  function  diffusion  process  (Case  II) 
explains  why  the  others  got  relatively  good  fits  to  their 
results  in  assuming  the  simple  profile.  It  also  explains  why 
for  longer  times,  the  simple  model  no  longer  could  explain  the 
results. 

Others  have  recently  used  this  model,  or  models  similar  to 
it,  in  analyzing  their  results.  One  example  is  Asbeck's  (128) 
analysis  of  outdiffusion  as  discussed  in  section  5.3  of  chapter 
5).  Muller  et  al.  (78)  assumed  a  two  species  diffusion  model 
for  Cr  when  they  theoretically  predicted  the  redistribution  of 
Cr  and  other  dopants  during  LPE  growth,  in  which  induced 
electric  field  effects  were  considered.  They  found  that  at  the 
impurity  levels  considered,  the  electric  field  had  virtually  no 
effect  on  Cr  movement. 

In  an  attempt  to  explain  the  results  of  Tuck's 
indiffusion  experiments  and  Kasahara's  outdiffusion 
experiments,  especially  the  apparently  contradictory  results 
when  arsenic  pressure  is  increased  (the  Cr  diffusion  apparently 
increased  in  one  case  and  decreased  in  the  other),  Muller  (129) 
proposed  a  model  similar  to  the  substitutional/interstitial 
dissociative  model,  except  assuming  that  the  substitutional 


chromium  species  is  sometimes  complexea  with  a  "D^s"  species, 
a  residual  donor  due  to  an  arsenic  atom  on  some  unusual  site; 
it  is  this  complexing  which  supposedly  renders  the  GaAs 
semi-insulating.  According  to  Muller,  the  Cr  is  not  complexed 
in  the  conditions  of  Tuck's  experiments  (indiffusions  into 
GaAs),  but  is  in  the  conditions  of  Kasahara's  experiments 
(outdif fusions  from  semi-insulating  Cr-doped  GaAs);  this 
difference  leads  to  the  different  behavior  of  Cr. 

It  should  be  pointed  out  that  Tuck  et  al.  (5,  120) 
concluded  that  before  their  results  could  be  interpreted,  the 
Ga-As-Cr  ternary  phase  diagram  must  be  known.  The  fact  that 
they  found  that  the  Cr  source  picked  up  large  amounts  of  Ga  and 
As  during  the  course  of  the  diffusion,  and  also  that  Cr 
concentrations,  upon  addition  of  arsenic  to  the  ampoule, 
decreased  instead  of  the  expected  increase  that  would  be 
predicted  based  on  a  Cr-on-Ga-si te  model,  suggested  to  them 
that  perhaps  the  Cr  was  reacting  with  Ga  and  As  and  forming  a 
new  phase.  The  need  for  the  Ga-As-Cr  phase  diagram  to  study  Cr 
in  GaAs  is  clear. 

In  summary,  much  work  has  been  done  studying  the 
redistribution  and/or  diffusion  of  Cr  in  GaAs.  Many  have 
calculated  diffusion  coefficients  (almost  all  for  short  times 
and  low  temperatures)  based  on  a  simple  error  function  model; 
yet  these  cannot  explain  the  deep  penetration  of  chromium  seen 
in  some  cases.  The  anomalous  Cr  build-up  at  the  GaAs  surface 
has  been  studied,  and  a  few  reasons  for  this  have  been 
suggested,  including  gettering  to  a  damaged  or  stressed  region, 
or  to  a  Ga-rich  region.  No  one,  except  for  Tuck  et  al.,  has 
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been  concerned  with  ternary  phase  equilbria  considerations. 

Most  experiments  were  done  in  open  furnace  tubes,  and  hence 
were  not  even  closed  systems.  And  Tuck's  experiments,  even 
though  the  importance  of  phase  equilibria  was  known,  were  done 
and  analyzed  without  the  actual  knowledge  of  it.  Finally,  a 
complex  diffusion  mechanism  was  suggested  for  Cr  in  GaAs,  but 
no  detailed  quantitative  study  had  been  done. 

6 . 3  Introduction  to  Diffusion  Studies 

The  purpose  of  this  study  of  Cr  diffusion  in  GaAs  was  to 
perform  diffusion  experiments  under  well-defined  conditions  and 
to  fit  the  results  to  a  diffusion  mechanism  model.  Three  types 
of  diffusion  experiments  were  done:  1)  Cr  indiffusions  in 
evacuated  quartz  ampoules,  with  the  selection  of  Cr  sources 
based  on  the  phase  equilibria  as  discussed  in  chapter  3;  2)  Cr 
outdif fusions ,  by  annealing  Cr-doped  GaAs  in  quartz  ampoules, 
with  and  without  well-defined  Cr  sinks;  and  3)  Cr  "junction" 
diffusions,  which  will  now  be  discussed.  Because  the 
substitutional/interstitial  dissociative  diffusion  mechanism 
was  strongly  suspected  -  based  on  Tuck’s  results  (5,  120)  as 
well  as  the  initial  results  of  the  present  work  -  a  slow,  and 
high  concentration,  diffusion  was  expected  in  the  surface 
region  for  indiffusion  (vacancy  diffusion  controlled). 

Likewise,  a  slow  diffusion  case  was  expected  at  the  surface  for 
outdif fusion.  However,  because  of  the  problems  in  the  GaAs 
surface  region  (anomalous  Cr  build-up  and  depletion,  to  be 
discussed  in  chapter  7),  these  diffusion  features  were 
obscured.  Therefore  interdiffusion,  or  "junction  diffusion," 
experiments  were  performed  in  which  the  diffusion  of  Cr  could 


be  studied  without  the  surface  phenomena  interfering.  This  was 
done  by  annealing  samples  of  undoped  GaAs  LPE  layers  on 
Cr-doped  GaAs  substrates.  The  "junction  sample"  (so-called 
because  of  the  sharp  junction  with  respect  to  the  Cr 
concentration  between  the  substrate  and  the  epilayer)  is  ideal 
for  Boltzmann-Matano  type  diffusion  experiments  and  analysis, 
and  both  outdif fusion  (out  of  the  substrate)  and  indiffusion 
(into  the  epilayer)  are  studied,  all  in  one  sample.  And  all 
this  occurs  away  from  the  GaAs  free  surface.  These  samples 
were  annealed  just  like  the  outdiffusion  samples,  with  or 
without  well-defined  Cr  sources  or  sinks. 

Considerations  of  ternary  phase  equilibria  are  relevant  to 
diffusion  studies  for  the  following  reasons:  1)  only  phases 
which  are  in  equilibrium  (or  "quasi-equilibrium,  as  previously 
discussed)  with  GaAs  should  be  used  as  Cr  sources  or  sinks;  2) 
the  Cr  solubility  (the  ultimate  level  of  Cr  in  the  GaAs  at  the 
end  of  a  diffusion,  and  the  level  of  Cr  at  or  near  the  surface 
for  all  times  for  both  indiffusions  and  outdif fusions )  depends 
on  the  location  of  the  overall  system  composition  on  the 
ternary  phase  diagram,  which  depends  on  the  amounts  of  the 
source  and  sink,  and  their  compositions;  and  3)  the  other 
intensive  variables,  such  as  C^',  Cs  ’ ,  and  Dcr,  also 
depend  on  the  location  in  the  phase  diagram.  However,  it  will 
be  shown  that  the  values  of  the  diffusion  coefficients  (both 
Dfast  an<3  Dslow)  were  found  not  to  vary  (within 
experimental  error)  over  the  range  of  composition  studied. 

This  will  be  discussed  later. 


6 . 4  Experimental  Techniques 

The  techniques  and  procedures  used  for  the  diffusion 
sutdies  were  the  same  as  those  used  in  the  solubility  studies 
and  discussed  in  chapter  4,  i.e.,  annealing  in  evacuated, 
closed  quartz  ampoules  with  pre-annealed  source  material,  and 
profiling  the  Cr  level  with  SIMS. 

The  undoped  GaAs  was  obtained  from  Cominco  (LEC,  high 
resistivity),  Hewlett-Packard  (LEC,  high  resistivity),  and 
Crystal  Specialties  (Bridgman,  low  resistivity).  The  Cr-doped 
GaAs  was  obtained  from  Crystal  Specialties,  the  Te-doped  from 
Hewlett-Packard  and  the  Si-doped  from  Hewlett-Packard  and 
Morgan  Semiconductor. 

The  LPE  "junction"  samples  were  obtained  from 
James  S.C.  Chang  at  Stanford  University  (see  his  Ph.D.  Thesis 
(130)  for  details  concerning  the  growth  of  these  -  some 
discussion  of  the  growth  is  given  in  the  next  chapter). 

All  samples,  including  the  junction  samples,  were  cleaned 
prior  to  each  anneal  according  to  the  schedule  given  in  chapter 
4.  Also,  each  was  cleaned  with  HC1  and  methanol  prior  to  SIMS 
analysis  to  clean  away  any  loose  particles  or  precipitates, 
except  when  specifically  studying  the  precipitates. 

6 . 5  Chromium  Diffusion  Experiments 
6.5.1  Introduction 

In  this  section,  the  results  of  the  the  three  types  of 
diffusion  experiments  -  indiffusion,  outdiffusion  and  junction 
diffusion  -  are  reported  and  discussed.  In  addition,  the 
implications  of  the  results  are  discussed,  including  how  the 
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reports  of  other  investigators  can  be  explained  in  terms  of 
these  results. 

6.5.2  Indiffusion  Experiments 

A  series  of  indiffusion  experiments  were  performed  for  a 
variety  of  times  and  temperatures,  using  H  source  and  Cominco 
undoped  GaAs.  The  temperatures  were  from  750  to  1000 °C  at  50 
degree  intervals,  and  the  times  ranged  from  1/2  to  12  hours  for 
the  lower  temperatures  and  1/4  to  2  hours  for  the  higher 
temperatures.  The  SIMS  Cr  profiles  for  the  800 °C  series,  shown 
in  figure  6-1,  are  the  results  of  sputtering  to  a  relatively 
shallow  depth  to  show  more  detail  of  the  surface  region.  In 
figure  6-2,  several  of  the  profiles  of  the  800  °C  samples  are 
shown  again,  this  time  profiled  to  a  much  greater  depth,  and 
figures  6-3  and  6-4  show  the  deep  Cr  profiles  of  the  850  and 
900°C  indiffusions. 

As  seen  in  the  shallow  profiles  of  the  800  °C  series 
(figure  6-1),  the  Cr  concentration  in  the  bulk  is  virtually 
level,  and  rises  uniformly  with  time.  As  discussed  earlier, 
this  suggests  that  the  substitutional/interstitial  dissociative 
diffusion  mechanism  is  operative.  As  seen  in  the  deep  profiles 
(figures  6-2,  6-3,  and  6-4),  the  "bulk"  Cr  levels  rise  with 
time,  as  in  figure  6-1,  but  here  one  can  see  that  the  Cr  levels 
tail  off  deep  in  the  interior  of  the  GaAs.  This  strongly 
resembles  Case  IV  of  the  diffusion  mechanism,  as  outlined  in 
chapter  5.  One  sees  that  the  bulk  rise  (to  the  final 
solubility  level)  is  faster  as  the  temperature  increases. 

Also,  the  chromium  penetration  is  much  deeper  than  one  would 
calculate  using  the  diffusion  coefficients  reported  in  most  of 
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Cr  CONCENTRATION  in  GaAs  (cm-3) 


Figure  6-4. 


Cr  concentration  profiles  from  900°C  indiffusion 
experiments  -  H  source. 


the  previous  work  (e.g.,  D  =  1  x  10“^  cm^/sec.  at  900°C 
according  to  (112)). 

In  most  of  the  indiffusion  profiles,  anomalous  features 
are  seen  in  the  surface  region  (0-4  pms  from  the  surface),  as 
discussed  earlier;  in  particular,  a  large  and  erratic  build-up 
of  Cr  right  at  the  surface  and  a  region  in  which  the  Cr  level 
is  less  than  the  bulk  level  just  to  the  right  of  it  in  the 
profiles.  These  features,  which  will  be  examined  as  to  their 
causes  in  chapter  7,  obscure  what  is  going  on  at  the  surface; 
however,  they  do  not  seem  to  affect  what  is  happening  in  the 
bulk.  (From  this  point  on,  the  "bulk"  shall  be  the  region  that 
is  inside  the  surface  build-up  and  depletion  region. ) 

Using  the  equation  for  Case  IV  diffusion,  uncharged 
species  (equation  5-18)  and  the  Cs '  values  from  the 
solubility  results,  the  profiles  were  fit  with  theoretical 
curves.  For  each  temperature,  Dfast  and  0  were  adjusted  to 
fit  the  profiles,  Dfast  being  the  diffusion  coefficient  that 
is  proportional  to  the  interstitial  diffusivity,  and  0  being 
the  vacancy  generation  parameter,  or  in  general,  the 
equilibrium  time  constant  which  describes  the  bulk  rise  to  the 
vacancy  equilibrium  case.  The  values  for  these  two  parameters 
for  each  temperature  were  used  for  all  the  different 
profiles,  i.e.,  for  the  different  times.  The  theoretical 
curves,  plotted  on  the  experimental  profiles  for  those  same 
three  temperatures  in  figures  6-5  to  6-7,  show  a  good  fit  to 
the  bulk  profiles.  The  Dfast  values  for  these  and  the  other 
temperatures  are  tabulated  in  table  6-1  and  are  plotted  in 
figure  6-8  (the  top  curve  in  the  figure),  showing  nice 
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Figure  6-7. 


900°C  Cr  profiles  with  theoretical  profiles 
(dotted);  D  =  5.7xlQ-^  cm^/sec,  q  *  3,000  sec 


Table  6-1  Chromium  indiffusion  results  (LEC  GaAs,  H  source) 


Temp 

Diast 

©_ 

750  °C 

1.5 

X 

10-10  cm^ /sec 

1.3 

x  10 5 

800 

2.8 

X 

10-10 

3.2 

x  10  4 

850 

4.0 

X 

10-10 

6.6 

x  103 

900 

5.7 

X 

10-10 

3.0 

x  103 
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Figure  6-8.  Dfast  and  Dslow  vs.  1/T  for  Cr  in  GaAs 


Arrhenius  behavior.  From  this,  an  activation  energy  was 
determined  to  be  0.83  (±  .07)  eV  and  a  DQ  value  of  2.04  x 
10~6  cm^/sec.  The  0  values  are  also  tabulated  in  table  6-1 
and  are  plotted  in  figure  6-9  (the  circled  points).  Again, 
Arrhenius  behavior  is  observed.  In  addition,  CQ,  the 
extrapolated  surface  concentration,  was  plotted  vs.  time  in 
figure  6-10,  along  with  theoretical  Cs '  ( 1-exp ( - t/ 9 ) )  curves, 
which  show  a  good  fit.  This  exponential  dependency  rather  than 
a  /t  one  is  good  evidence  that  the  vacancy  conversion  mechanism 
is  operating  rather  than  the  kickout  mechanism,  as  discussed  in 
section  5.6  of  chapter  5.  (However,  the  kickout  mechanism 
cannot  be  ruled  out  solely  by  this  evidence  because  it  also 
predicts  an  exponential  time  dependency  as  the  concentration 
approaches  the  solubility  level;  however, it  is  unclear  when 
this  transition  from  /t  to  exponential  behavior  occurs.) 

As  discussed  earlier,  the  models  for  charged  species  give 
much  different  results  than  for  uncharged  species;  in  this  case 
the  uncharged  model  clearly  fits  best;  the  profiles  for  the 
charged  case  would  fall  off  abruptly  at  about  C  =0.5  CQ,  as 
opposed  to  the  gradual,  error  function- type  drop  in 
concentration  as  the  depth  increases.  This  is  consistent  with 
substitutional  Cr  being  uncharged  in  undoped  (semi-insulating) 
GaAs,  as  discussed  in  chapter  4. 

As  mentioned  in  chapter  5,  one  must  consider  the  intrinsic 
carrier  concentration,  because  if  it  is  much  higher  than  the 
impurity  concentration  (at  the  diffusion  temperature),  then  the 
diffusion  process  will  not  be  concentration  dependent  even  for 
charged  species.  In  GaAs,  n^  =  4  x  10^6  cm-^  at  750°C, 
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Figure  6-10.  CQ  vs.  time.  CQ  is  the  bulk  Cr  concentration, 

extrapolated  to  the  surface.  Theoretical  curves  are 
dotted. 
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and  the  solubility  of  Cr  is  -2  x  10^-6  cm-3  at  that 
temperature.  Therefore,  substitutional  Cr  could  be  charged  and 
still  not  exhibit  concentration  dependent  diffusion.  This  will 
be  considered  again  in  the  junction  diffusion  samples.  In  any 
case,  the  uncharged  model  fits  very  well. 

In  addition  to  a  good  fit  of  the  experimental  profiles  to 
the  theory,  the  Ea  and  D0  values  of  .83  eV  and  2.04  x 
10“6  cm^/sec,  respectively,  provide  additional  confirmation 
of  the  model  since  they  are  typical  values  for  the 
substitutional/interstitial  dissociative  diffusion  mechanism, 
as  discussed  in  chapter  5.  Other  results,  to  be  discussed  in 
this  chapter  and  the  next,  also  lend  support  to  this  diffusion 
model  for  Cr  in  GaAs. 

Since  vacancies  can  diffuse  in  from  the  surface,  one  would 
expect  to  see  Case  II  diffusion  at  the  surface,  with  the  Cr 
concentration  right  at  the  surface  equal  to  the  Cr  solubility 
value  for  the  temperature  in  question  -  6  x  10^6  cm~3  at 
850°C,  for  example  -  and  a  shallow  error  function  profile  down 
to  the  position  where  the  "bulk  Cr"  is  reached  (due  to  the  Case 
IV  diffusion).  However,  as  mentioned,  the  unusual  surface 
phenomena  obscure  this.  In  a  couple  of  profiles  (the  25  min., 
850°C  and  the  25  min.,  800°C  profiles,  for  instance),  this 
diffusion  case  is  evident.  From  these,  diffusion  coefficients 
for  800  and  850°C  were  determined,  1  x  lO”-^  and  7  x  10~^2 
cm^/sec.  respectively. 

A  diffusion  was  performed  in  doped  GaAs  (Te  doped, 
n  -  4  x  10^7  cm_2 ) ,  to  see  if  doping  has  an  effect  on  the 
diffusion  of  chromium.  The  SIMS  Cr  profile  for  the  2  hr.. 
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850°C  diffusion  (with  H  source)  is  shown  in  figure  6-11. 
Although  the  extrapolated  surface  concentration  of  9  x  10^-® 
cm~3  j_s  higher  than  the  Cr  solubility  in  undoped  GaAs,  as 
discussed  in  chapter  4,  the  effective  diffusion  coefficient 
(DfaSt)  of  *2  x  10“10  cm2/sec.  is  less  than  in  the 
undoped  case  (where  D  =  4  x  10“^^  cm2/sec).  This  can  be 
explained  by  the  model.  By  equation  5-9,  Dfast  =  Di 
C i  '/ ( C i ' +CS 1 ) .  In  highly  doped  n-type  GaAs,  Cr  is 
charged,  therefore,  the  conversion  equation  is  not: 

Ii  +  V  =  Is  (5-1) 

but  either: 

Ii+v=Is~+  h+  (6-1) 

or : 

Ii  +  V  +  e"  =  Is-  (6-2) 

Therefore,  increasing  the  electron  concentration  (by  increasing 
the  doping  level)  shifts  equation  6-2  to  the  right.  This 
decreases  the  ratio  Ci/Cs,  and  therefore  decreases 
Dfast*  (This  is  also  the  reason  why  Cs,  which  is  the  Cr 
solubility,  increases  with  n-type  doping.) 

Next,  Cr  diffusions  were  done  into  Crystal  Specialties 
undoped  (Bridgman)  GaAs,  for  various  times  at  750°C.  There 
were  two  distinct  differences  between  the  results  for  these  and 
for  the  Cominco  LEC  GaAs:  1 )  0  is  much  smaller  for  the  Crystal 
Specialties  (2.4  x  104  sec.  vs.  1.3  x  105  sec.)  meaning 
that  the  solubility  level  was  reached  much  more  quickly;  and  2) 
in  some  cases,  the  tail  end  of  the  Cr  profile  dropped  off  more 
quickly.  Figure  6-12  shows  four  SI MS  Cr  profiles  for  the 
7 50 °C ,  4  hr.  Cr  indiffusion,  all  from  the  same  sample  (also 
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Figure  6-11.  Cr  concentration  profile  from  850aC, 
indiffusion  experiment  in  Te-doped  ( 
cm“3 )  GaAs  -  H  source. 
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Figure  6-12.  Cr  concentration  profiles  from  75G°C,  4  hr. 

indiffusion  experiment  in  Crystal  Specialities 
(Bridgman)  GaAs  -  H  source.  The  three  profiles 
from  different  locations  on  same  GaAs  wafer. 


shown  is  the  theoretical  curve  for  the  uncharged  case).,  A  hint 
to  the  cause  for  this  is  the  fact  that  usually  3r idgman-type 
GaAs  is  not  semi- insulating  due  to  impurities  (usually  oxygen 
and  silicon)  and  defects,  and  it  is  usually  n-type  (the  reason 
why  Cr  must  be  added  to  it  to  make  it  semi-insulating),  whereas 
LEC  GaAs  is  intrinsically  semi-insulating.  Unannealed  Bridgman 
GaAs  was  checked  by  C-V  analysis,  and  it  was  found  that  the 
carrier  concentration  (n-type)  varied  from  point  to  point  on 
the  wafer,  from  a  value  of  -1  x  10^4  cm_3  to  the  mid 
1016's.  These  non-uniform  carrier  concentrations  can  explain 
the  diffusion  profiles  in  the  Bridgman  samples.  In  the  regions 
where  the  electron  density  is  low,  the  normal  uncharged 
diffusion  occurs.  However,  in  regions  of  moderate  electron 
density,  the  diffusion  is  influenced  by  the  presence  of  the 
electrons.  Near  the  surface,  where  the  Cr  level  is 
approximately  the  same  as  the  electron  concentration,  the  Cr 
compensates  the  free  electrons  and  the  material  becomes 
semi-insulating.  However,  deeper  in  the  sample,  the  Cr  level 
is  lower  and  the  net  electron  density  is  greater  than  zero. 
Thus,  by  the  arguments  used  earlier  regarding  a  shift  in  the 
conversion  reaction  due  to  a  change  in  the  electron  or  hole 
concentration,  the  effective  diffusion  coefficient  is 
decreased.  The  extent  of  the  decrease  in  the  diffusivity  is 
directly  related  to  the  electron  density;  hence,  for  different 
electron  densities,  the  extent  of  the  diffusion  reduction  will 
be  different,  resulting  in  different  shaped  curves.  To  test 
this,  a  Cr  diffusion  was  performed  with  a  Crystal  Specialty 
(Bridgman)  wafer  at  1000°C.  At  that  temperature,  the  intrinsic 
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carrier  concentration  is  =  5  x  lO-'-7  cm-^  (as  compared  to  =4 
x  10^6  at  750°C)  and  therefore,  any  extrinsic  electrical 
affect  should  be  eliminated.  Figure  6-13  shows  the  resulting 
SIMS  Cr  profile,  along  with  the  profile  for  a  similar  run  with 
Cominco  LEC  GaAs.  As  one  can  see,  the  profiles  are  virtually 
identical;  the  Bridgman  sample  profile  does  not  exhibit  the 
abrupt  drop-off  in  Cr  concentration  observed  in  the  lower 
temperature  samples.  Therefore,  there  is  good  evidence  that 
the  difference  in  Cr  profiles  observed  is  related  to 
electrically  active  impurities  or  defects  in  the  Bridgman  GaAs. 

The  reason  for  the  smaller  value  for  0  is  not  clear.  At 
first  one  would  think  that  since  LEC  GaAs  is  usually  grown  with 
a  slight  excess  of  arsenic,  whereas  Bridgman  GaAs  is  usually 
grown  exactly  stoichiometric  (to  avoid  high  arsenic  pressures), 
that  LEC  GaAs  would  have  a  higher  initial  concentration  of 
gallium  vacancies.  Therefore,  the  build-up  to  the  equilibrium 
value  for  that  reason  should  be  faster  for  the  LEC  material 
rather  than  the  Bridgman  material.  This,  of  course,  assumes 
that  Cr  is  incorporated  onto  a  gallium  site  rather  than  an 
arsenic  site.  To  test  this  assumption,  a  piece  of  LEC  GaAs 
was  pre-annealed  in  an  ampoule  with  a  piece  of  elemental 
arsenic  at  800°C  for  24  hours.  This  should  increase  the 
gallium  vacancy  concentration  and  decrease  the  arsenic  vacancy 
concentration.  Then  a  800°C  Cr  indiffusion  was  performed  for  4 
hours.  The  bulk  level,  as  determined  by  SIMS,  was  indeed 
hijher  than  in  the  normal  LEC  case,  with  a  0  value  of  1/2  the 
original.  Therefore,  the  original  assumptions  regarding  the 
diffusion  mechanism  seem  to  be  good.  Another  possibility  is 
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that  the  dislocation  densities  might  be  different;  more 
dislocations  would  generate  more  vacancies,  and  hence  0  would 
be  smaller.  But  the  dislocation  densities  (as  measured  by  the 
etch  pits),  are  about  the  same  (*5  x  10^  cm~2  for  the  LEC 
and  =>2  x  10^  for  the  Bridgman),  if  anything,  lower  for  the 
Bridgman.  However,  these  dislocation  densities  refer  to  all 
types  of  dislocations,  but  there  are  different  types  of  dis¬ 
locations,  such  as  screw,  edge,  edge  that  ends  with  a  Ga  atom, 
edge  that  ends  with  an  As  atom,  and  not  all  generate  gallium 
vacancies.  Therefore,  more  of  the  type  that  creates  gallium 
vacancies  may  be  present  in  the  Bridgman  GaAs  than  in  the  LEC 
GaAs.  Another  possibility  is  that  whatever  is  the  source  of 
the  net  electron  concentration  in  the  Bridgman  material  (oxygen 
impurities,  for  example)  might  also  affect  the  0  value. 

At  this  point,  the  exact  cause  of  this  difference  in  Cr 
behavior  is  not  clear.  For  whatever  reason,  this  difference  in 
0  for  different  types  of  GaAs  is  good  evidence  that  the  rise  in 
the  bulk  Cr  concentration  level  is  associated  with  processes 
occuring  inside  the  GaAs,  as  the  model  describes,  rather  than 
on  changing  surface  conditions  due  to  changing  source  and/or 
vapor  phase  conditions  (which  could  result  in  similar 
profiles).  In  addition,  this  difference  in  Cr  behavior,  as 
well  as  the  changes  in  diffusivity,  shows  how  Cr  can  be  used  as 
a  probe  to  study  such  things  as  vacancies  and  other  defects  and 
impurities,  which  are  normally  difficult,  if  not  impossible,  to 
study  directly. 

Finally,  a  Cr  indiffusion  was  done  into  GaAs  (LEC)  at 
800°C  for  12  hours,  this  time  using  J  source.  As  discussed 
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earlier,  this  results  m  a  lower  Cr  solubility,  i.e., 
substitutional  equilibrium  concentration,  than  using  H  source; 
here  any  difference  in  diffusion  was  studied.  From  the  SIMS  Cr 
profile  (figure  6-14),  a  diffusion  coefficient  of  4  x  10“^^ 
cm^/sec.  was  determined.  A  repeat  of  this  experiment 
resulted  in  the  same  Dfast  value.  This  value  is  just  barely 
within  the  (large)  expected  error  of  the  300°C,  H  source  value 
of  2.8  x  10“i0  (see  figure  6-8),  and  therefore  one  cannot 
really  say  that  the  value  of  Dfagt  changed  from  source  H  to 
J.  However,  a  higher  value  for  J  _is  reasonable  because  as  the 
source  changes  from  H  to  J,  the  value  for  Cs ' ,  the 
substitutional  solubility,  decreases  almost  an  order  of 
magnitude,  and  C^',  the  interstitial  solubility,  would 
probably  not  change  as  much.  Therefore,  Dfast,  which  equals 
Dici  * / ( C  ' +CS  '  ) ,  would  be  expected  to  slightly  increase 
as  the  source  composition  moves  towards  the  Ga  corner. 


6.5.3  Qutdiffusion  Experiments 

Well-defined  Cr  outdiffusion  experiments  were  performed  by 
annealing  Cr-doped  GaAs  (Crystal  Specialties)  in  an  evacuated 
closed  quartz  ampoule,  which  also  contained  a  large  amount  of 
"sink"  material.  This  is  the  same  idea  as  using  the  "source" 
material  in  the  indiffusions,  except  that  the  composition  of  it 
is  such  that  the  Cr  solubility  in  GaAs ( s )  for  that  system  is 
less  than  the  initial  Cr  level  in  the  GaAs.  Therefore,  Cr  will 
diffuse  out  of  the  GaAs  and  into  the  sink.  An  example  would  be 
annealing  GaAs,  with  1  x  lO^-7  cm-3  Cr,  with  J  source  at 
800 °C.  If  there  is  a  lot  more  sink  material  than  GaAs,  then 
the  total  system  composition  will  lie  near  the  composition  of 
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J.  The  GaAs  with  =>5  x  10^  cra'^  Cr  is  in  equilibrium  with 
the  sink,  and  therefore,  Cr  will  diffuse  out  of  the  GaAs  until 
that  concentration  level  is  reached  throughout.  During  the 
diffusion,  the  surface  concentration  of  Cr  should  be  at  the 
solubility  limit  (5  x  10^5  cm-^  in  this  case),  as  is  the 
case  with  indiffusions  as  well.  The  problem  with  doing 
outdif fusions  without  a  well-defined  sink  is  discussed  a  little 
later. 

Outdif fusions  were  performed  using  J  source  at  700,  800 
and  850 °C  (for  2,  1,  and  1  hr.,  respectively),  with 
the  corresponding  SIMS  Cr  profiles  shown  in  figure  6-15a.  The 
extent  of  outdiffusion  is  similar  to  that  of  the  assumed  Case 
II  for  indiffusion  (vacancy  diffusion  controlled);  therefore, 
error  function  diffusion  curves  were  fitted  to  the  profiles  to 
determine  Dsiow  values.  As  discussed  in  chapter  5,  solving 
equation  5-10  for  outdiffusion  conditions  leads  to  the  usual 
outdiffusion  equation,  (C-Cj  )/(C2-C  -L  )  =  erf  ( x/ ( 4Dt )  )  , 

where  D  in  this  case  is  "Dsiow"  and  equals  Dv  CV'/(DV' 

+  Cs '  ) ;  C 2  the  initial  bulk  Cr  level,  and  the  surface  is 
maintained  at  a  concentration  of  Clf  in  this  case  the 
solubility  of  GaAs  for  the  J  sink.  By  making,  the  best  fit  to 
the  profiles,  shown  in  figure  6-15b,  Ds]_ow  values  were 
determined;  they  are  tabulated  in  table  6-2  and  are  plotted  in 
figure  6-8.  As  one  can  see,  they  are  very  close  to  the  Dsiow 
values  determined  from  the  indiffusion  experiments.  It  should 
be  noted  that  the  fits  to  the  profiles  are  not  perfect, 
especially  for  the  850 °C  profile  over  the  first  micron  from  the 
surface;  this  will  be  discussed  later.  Also,  the  build-up  of 
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Table  6-2.  Chromium  outdif fusion  results. 


Temp 
700  °C 


Time 


2  hr. 


S  ink 


■2s  low 

2  x  10"13  cm2/sec 
1.6  x  10~12 
1.4  x  10“10 

1.4  x  10“2 

1.5  x  10"12 

1.6  x  10~12 
2.0  x  10"12 

2.2  x  10"12 
1.8  x  10'12 
1.0  x  lO'11 

1.3  x  10-11 
2.0  x  10'11 
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Cr  at  the  surface  was  ignored  in  these  analyses;  this  topic  is 
discussed  in  chapter  7. 

Additional  outdiffusion  experiments  were  performed  for 
different  annealing  times,  for  example,  1/2  and  2  hours  at 
800 °C,  and  4  hours  at  850 °C.  The  Ds]_ow  values  determined 
from  these  were  each  within  20%  of  the  other  values  for  those 
temperatures  -  within  the  expected  error. 

SIMS  profiles  were  obtained  from  several  of  the  samples 
before  they  were  annealed,  and  their  bulk  levels  were  compared 
to  the  bulk  levels  of  the  profiles  after  the  outdif fusions.  In 
each  case,  the  level  was  within  the  expected  error  of  the 
concentration  determinations  (±  50%);  however,  in  virtually  all 
the  cases  (90%),  the  bulk  level  after  the  anneal  was  at  or 
below  the  initial  bulk  level.  If  it  was  just  a  random 
fluctuation,  then  in  more  cases  than  this  the  bulk  level  after 
the  anneal  would  have  been  greater  than  the  initial  level. 
Therefore,  I  suspect  that  the  bulk  level  may  drop  a  little 
during  the  anneal. 

The  mechanism  of  Cr  outdiffusion  can  be  explained  in  terms 
of  the  proposed  model.  During  outdiffusion,  as  in  indiffusion, 
the  atoms  diffuse  by  converting  from  substitutional  to 
interstitial.  A  certain  number  make  this  conversion  initially 
and  can  quickly  diffuse  out.  (This  could  explain  the  drop  in 
the  bulk  level.)  However,  now  the  GaAs  is  saturated  with 
vacancies,  and  no  more  atoms  can  convert  to  inters’  _ial  and 
diffuse  out  until  the  vacancies  can  be  disposed  of.  This  can 
happen  either  by  the  vacancies  diffusing  to  the  surface,  which 
is  a  natural  vacancy  sink,  or  they  can  be  annihilated  in  the 
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bulk,  such  as  by  appropriate  dislocation  movement.  This  is 
completely  analogous  to  the  indiffusion  case,  and  it  is  why  the 
same  governing  equations  hold  for  each  case,  with  only  initial 
and  boundary  conditions  different.  The  fact  that  the 
switchover  from  Case  II  dominating  to  Case  I  dominating  is  not 
observed  -  even  after  4  hours  at  850  °C  -  is  probably  because 
the  bulk  vacancy  annihilation  mechanism  is  not  as  effective  as 
the  bulk  vacancy  generation  mechanism.  (As  mentioned,  other 
researchers  (112,  108)  did  observe  an  abrupt  increase  in  the 
effective  diffusion  coefficient  after  -  1  hour  at  800°C, 
indicating  a  transition  from  Case  II  to  Case  I,  probably  due  to 
a  different  defect  state  of  their  starting  material.) 

Outdif fusion  runs  were  also  done  at  800 °C,  for  1  hour 
each,  using  different  sinks:  Jl,  J2 ,  and  D,  whose  compositions 
are  shown  in  figure  4-4.  (Note:  Jl  is  the  same  as  J,  used 
earlier. )  This  results  in  different  surface  concentrations  for 
the  same  temperature.  In  the  case  of  "sink"  D,  the  solubility 
determined  by  D  is  above  the  initial  Cr  bulk  level,  thus  an 
indif f us  ion  would  result,  and  D  is  a  source  in  this  case  rather 
than  a  sink.  The  SIMS  profiles  are  shown  in  figure  6-16.  As 
expected,  the  surface  concentration  increases  as  the 
sink/source  composition  goes  from  J2  to  Jl  to  D.  The 
appropriate  theoretical  equation  (indiffusion  or  outdif fusion ) 
was  fit  to  each  curve,  and  only  Cx  was  changed  from  one  curve 
to  another,  and  not  D.  The  fits  are  quite  good.  Therefore, 
while  the  solubility  of  Cr  in  GaAs  changes  over  this  range  of 
composition  in  the  ternary  phase  diagram,  Ds]_ow  does  not 
appear  to.  Since  DS]_ow  =  Dv  CV'/(CV'  +  Cs ' ) ,  and 
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Figure  6-16.  Cr  concentration  profiles  from  800°C,  45  min. 

outdiffusion  experiments  with  various  sinks,  as 
indicated . 


increases  when  goin<_  from  J2  to  J1  to  D,  it  appears 


I 

since  Cs ' 

I  that  Cv '  increases  accordingly  so  that  the  net  change  in  Dv 

I 

is  small.  This  is  reasonable  since  one  would  expect  that  the 
equilibrium  gallium  vacancy  concentration  would  increase  as  the 
j  composition  moves  away  from  the  Ga  corner  of  the  diagram. 

For  comparison,  some  outdiffusion  experiments  were  done 
without  a  sink  in  the  ampoule  (just  annealing  Cr-doped  GaAs  in 
I  the  evacuated,  closed  ampoule).  The  SIMS  profile  for  a  750°C, 

i 

1  hr.  outdiffusion  is  shown  in  figure  6-17.  This  is  very 
similar  to  the  profile  obtained  using  J1  or  J2  sinks.  An  error 
function  curve  was  fitted  to  this  (ignoring  the  surface 
build-up  and  assuming  a  C1  value  of  1  x  10 cm~^ ) , 
resulting  with  a  diffusion  coefficient  consistent  with  the 
other  outdif fusions .  This  result  is  also  consistent,  both  in 
surface  concentration  and  in  the  diffusion  coefficient,  with 
most  of  the  results  of  other  researchers,  who  also  did 
short-time  outdif fusions  without  a  well-defined  sink.  This 
will  be  discussed  further  in  section  6.6. 

The  reason  why  the  surface  concentration  is  similar  to 
that  for  outdiffusing  when  the  J  sink  is  used  can  be  explained 
in  terms  of  the  phase  diagram  and  phase  equilibria.  At  first , 
one  would  think  that  because  the  system  is  made  up  of  only 
Cr-doped  GaAs,  in  which  the  Cr  level  is  above  or  at  the 
solubility  level  for  the  annealing  temperature,  then  it  would 
lie  in  the  GaAs ( 3 )-CrAs ( s )-Ga ( ^  )  tie  triangle  region  (see 
figure  6-13).  Therefore,  no  outdiffusion  would  occur. 

However,  as  mentioned  in  chapter  2,  one  must  take  into  account 
the  vapor  pressure.  Because  there  is  no  sink  or  source,  all 
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after  each  anneal.  A  typical  result  is  shown  in  figure  6-19, 
before  and  after  a  1  hour,  800°C  anneal.  Three  observations 
are  made:  1)  there  is  "normal"  Cr  diffusion  evident  from  the 
substrate  into  the  LPE  layer  (assumed  to  be  the  vacancy 
diffusion  controlled  case);  2)  the  bulk  level  on  the  substrate 
uniformly  decreases  a  little,  as  in  the  outdiffusion  case;  and 
3)  Cr  builds  up  at  the  surface,  which  is  evidence  of  fast 
diffusion  of  Cr  from  the  substrate,  through  the  LPE  layer  to 
the  surface.  These  features  are  characteristic  of  all  the 
junction  diffusion  samples. 

The  decrease  in  the  bulk  level  is  believed  to  arise  from 
the  same  mechanism  as  seen  in  the  outdiffusion  samples:  the 
initial  diffusion  of  interstitial  chromium  occurs  until  the 
vacancy  concentration  is  saturated.  The  "normal"  diffusion  at 
the  interface  is  assumed  to  be  the  vacancy  controlled  Cr 
diffusion.  The  proposed  mechanism  is  very  similar  to  the 
outdiffusion  case,  except  that  it  includes  both  outdiffusion 
and  indiffusion  -  outdiffusion  from  the  substrate  and 
indiffusion  into  the  LPE  layer.  The  fast  interstitials  quickly 
diffuse  out  of  the  substrate  and  into  the  LPE  layer  (and  to  the 
surfuoj ,  and  perhaps  even  out  of  the  sample).  Once  the 
substrate  is  saturated  with  vacancies,  no  more  interstitial  Cr 
can  be  produced  until  the  vacancies  are  removed.  Since  there 
is  no  evidence  for  efficient  bulk  annihilation,  this  can  only 
occur  by  vacancy  diffusion  from  the  substrate  and  into  the  LPE 
layer.  Therefore,  the  vacancy  diffusion  controlled  Cr 
diffusion  is  seen  both  out  of  the  substrate  and  into  the  LPE 
layer,  where  the  interstitials  can  recombine  with  the  vacancies 


and  become  substitutionals.  This  results  in  one  continuous 


diffusion  profile  across  the  junction. 

This  type  of  experiment  is  ideal  to  measure  the  effective 
diffusion  coefficient  as  a  function  of  concentration  because, 
if  one  assumes  a  semi-infinite  medium  on  each  side,  one  can  use 
the  Boltzmann-Matano  analysis  to  determine  D(C)  (131,  132). 
Shewmon  (94)  gives  a  detailed  review  of  this  graphical 
technique.  The  resulting  equation  for  D(C’),  for  the 
conditions:  C(x<0,t=0)  =  Cj[  and  C(x>0,t=0)  =  C2  and 

semi- inf inite  media,  is: 

D(C '  )  =  —  (—)  jC‘  x  dC  (6-3  ) 

2 1  dC  C'  C1 

where  x  =  0  is  defined  by 

C2 

/  x  dC  =  0  (6-4) 

C1 

Therefore,  by  measuring  the  appropriate  area  and  slope  one  may 
obtain  D  values  as  functions  of  concentration.  This  analysis 
was  done  on  all  the  samples.  It  was  found  that,  within  error 
of  the  technique,  D  does  not  vary  with  concentration.  The 
Dsiow  values  determined  are  also  tabulated  in  table  6-3. 

In  regards  to  the  question  of  the  Cr  being  charged  or  not, 
the  diffusion  coefficients  for  the  70C|OC  sample  were  determined 
to  be  within  15%  of  each  other  from  concentrations  levels  of 
4  x  1Q15  to  3  x  1C)16  cm-3  ,  while  nj_=  2  x  10^-6  cm-^  , 
thus  indicating  that  the  substitutional  Cr  is  not  charged. 
Therefore,  the  concen tra t  ion- independent  diffusion  can  not  be 
attributed  to  just  the  nL  effect,  since  n^  is  not  >>  Cs. 


.  0  ^ 
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Since  the  diffusion  coefficient  was  determined  to  be 


independent  of  concentration,  then  Case  II,  uncharged,  of  the 
diffusion  model  should  apply.  This  is  simple  "error  function” 
diffusion  applied  to  the  conditions  mentioned  previously, 
resulting  in: 

C  (  x ,  t )  =  (Ca-C2)/2  •  erf  (x/2  (D^t)1/2  )  (6-5) 

This  was  fitted  to  each  profile  with  the  appropriate  Ds]_ow 
value,  as  determined  by  the  Boltzmann-Matano  analysis.  Two 
examples  of  the  fit  are  shown  in  figures  6-20  and  6-21.  As  one 
can  see,  the  fits  are  extremely  good,  indicating  that  D  is 
independent  of  concentration,  and  that  the  model  quite 
adequately  describes  the  behavior.  Furthermore,  the  fits  are 
much  better  in  these  cases  than  for  the  outdiffusion 
experiments,  especially  at  the  higher  temperatures,  indicating 
that  the  diffusion  at  the  surface  is  probably  being  affected  by 
other  phenomena;  this  confirms  the  importance  of  doing  these 
junction  diffusion  experiments  to  study  the  diffusion  of  Cr  in 
GaAs . 

The  above  analysis  assumes  an  exact  step  formation  of 
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the  Dt  of  the  final  profile,  and  dividing  by  the  time  (and  Dj 
does  not  even  have  to  have  been  constant,  or  t  L  need  to  be 
known,  to  do  this). 

This  was  done  for  each  sample.  It  is  these  new  corrected 
D's  that  are  tabulated  (the  D's  for  the  short  time/low 
temperature  anneals  are  about  35%  smaller  than  the  uncorrected 
values  -  the  D's  for  the  long  time /high  temperatures  are  only 
about  10%  smaller). 

The  diffusion  was  performed  for  various  times  at  800  and 
900 °C,  and  the  D  values  determined  did  not  change  with  time. 
This  indicates  that  the  onset  of  Case  I  diffusion  did  not  occur 
for  the  times  used,  and  that  the  diffusion  coefficients 
determined  are  valid.  (The  D  values  would  not  remain  constant 
with  time  if,  for  example,  the  bulk  level  changed  after  the 
initial  drop. ) 

The  average  Dsiow  value  for  each  temperature  is  plotted 
in  figure  6-8.  These  are  quite  similar  to  those  determined 
from  the  outdif fusions  and  the  indiffusions  at  the  surface,  as 
expected . 

From  the  junction  diffusion  Ds]_ow  values  for  different 
temperatures,  values  for  Ea  and  D0  were  determined:  Ea  = 

1.7  eV  and  D0  =  1.7  x  10“^  cm^/sec.  These  are  very 
reasonable  numbers  for  the  vacancy  diffusion  controlled 
diffusion  case,  as  discussed  in  chapter  5.  (The  Ds^ow  values 
from  the  junction  diffusions  are  assumed  to  be  more  accurate 
than  the  outdiffusion  or  indiffusion  Ds^ow  value,  which  were 
determined  from  the  profiles  at  the  surface;  therefore  it  was 
with  these  values  from  which  E=  and  D0  were  determined.) 


Using  equation  5-12  (Dsi0w  3  Dv  CV’/(CV'  +  Cs ’ ) ) ,  one 
can  calculate  Dv  values  for  the  different  temperatures  from 
the  following  selected  values:  the  present  values  for  Dsiow' 
the  present  solubility  values  for  Cs ' ,  and  the  Cv'  values 
from  Chang  and  Pearson  (83).  These  results,  tabulated  in  table 
6-4,  are  within  an  order  of  magnitude  of  the  vacancy 
diffusivity  values  determined  by  Chang  and  Pearson  (83).  The 
activation  energy  for  Dv  was  determined  to  be  1.9  eV,  very 
close  to  the  value  of  2.1  eV  also  determined  by  Chang  and 
Pearson.  Therefore,  it  is  reasonable  to  conclude  that  this  Cr 
diffusion  is  limited  by  gallium  vacancy  diffusion. 

As  mentioned  previously,  evidence  for  the  fast 
interstitials  is  seen  in  the  build-up  of  Cr  at  the  surface, 
many  microns  away  from  the  source  of  Cr  in  some  cases.  If  the 
Cr  diffusion  were  just  via  the  normal  one-species  diffusion 
(which  also  usually  follows  "error  function  diffusion”),  then 
one  would  only  observe  the  diffusion  profile  at  the  junction. 
But  because  there  are  the  fast  interstitials,  which  can 
concentrate  by  converting  to  substitutional  or  by  dissolving  in 
another  phase  (as  will  be  discussed  in  chapter  7),  a  relatively 
large  amount  of  Cr  is  seen  far  away  from  the  Cr  source  after 
short  annealing  times. 

Mosu  of  the  previously  mentioned  junction  diffusion 
anneals  feere  done  with  no  source  or  sink  in  the  ampoule.  A  few 
were  done  with  J2  sink  (for  a  surface  concentration  of  -lO^-^ 
cm~3 ) t  wi  :h  no  observable  difference  in  the  results.  This  is 
because  tie  Cr  concentration  in  the  LPE  layer  was  at  that  level 
to  begin  with  and,  just  like  the  outdif fusions  done  without  a 


source  or  sink  (as  explained  in  terms  of  the  phase  diagram  in 
section  6.5.3),  it  would  remain  at  that  level.  That  is 


probably  one  reason  why  one  did  not  see  an  increase  in  the  bulk 
Cr  level  in  the  LPE  layer  (Case  III  diffusion),  as  was  seen  in 
the  indiffusion  experiments;  the  equilibrium  solubility  level 
was  already  achieved.  Another  possible  reason  is  that  the  LPE 
GaAs  might  not  be  able  to  generate  vacancies  as  effectively  due 
to  a  smaller  number  of  dislocations.  Or  perhaps  it  is  because 
the  initial  gallium  vacancy  concentration  in  the  LPE  layer, 
which  was  grown  from  a  Ga-rich  melt,  is  very  low.  Any  or  all 
of  these  factors  could  explain  the  bulk  Cr  level  in  the  LPE 
layer  not  rising. 

The  fact  that  the  diffusion  coefficient  is  evidently 
independent  of  concentration  is  further  proof  that  the 
Frank-Turnbull  mechanism  for  the  interstitial-to-substitutional 
conversion,  rather  than  the  kickout  mechanism,  is  operating. 

As  discussed  in  chapter  5,  the  kickout  mechanism  supposedly 
predicts  a  C"2  dependence  for  the  slow  diffusion  coefficient, 
while  the  F-T  mechanism  predicts  a  concentrat ion- independent 
Dsiow*  The  data  clearly  supports  the  F-T  model  over  the 
other  in  this  case.  This  mechanism  is  reasonable  for  Cr  in 
GaAs  because:  1)  the  Cr  solubility  is  very  close  to  the 
equilibrium  gallium  vacancy  concentration,  and  therefore 
relative  few,  if  any,  vacancies  must  be  generated;  and  2)  the 
effective  diffusion  coefficient  (Ds^ow)  measured  is  very 
close  to  the  vacancy  diffusion  coefficient  at  each 
temperature.  The  two  main  arguments  used  against  this 
mechanism  were  that  too  many  vacancies  were  required  to  be 
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generated  and  the  effective  diffusion  coefficients  were  much 
larger  than  the  reported  vacancy  diffusion  coefficients.  These 
arguments  are  not  valid  in  the  present  case,  and  the  F-T  model 
is  therefore  quite  reasonable. 

6 . 6  Chromium  as  a  Probe 

Because  of  its  great  sensitivity  to  its  environment,  as 
evidenced  in  the  results  so  far,  Cr  may  be  used  as  a  "probe"  to 
study  the  defect  chemistry  of  GaAs.  It  has  been  shown  in  the 
present  and  previous  studies  that  the  solubility  and  diffusion 
of  Cr  is  greatly  affected  by  other  impurities,  vacancies, 
precipitates,  and  dislocations,  and  thus  by  studying  the 
behavior  of  chromium  in  GaAs,  one  can  study  the  defect 
chemistry  of  GaAs.  In  the  next  chapter  it  is  shown  how  high 
levels  of  another  impurity  were  discovered  as  a  result  of 
studying  Cr. 

6 . 7  The  Analysis  of  Previous  Studies 

These  results  and  the  diffusion  model  can  now  be  used  to 
explain  some  of  the  results  and  observations  of  other 
investigators  concerning  Cr  in  GaAs. 

First,  the  results  and  model  explain  the  deep  penetration 
of  Cr  in  GaAs  that  some  have  observed  (5,  102),  penetration 
much  deeper  than  the  "reported  diffusion  coefficients"  could 
account  for.  This  of  course  is  due  to  the  fast  interstitials, 
which  can  convert  to  substitutionals  wherever,  and  whenever, 
vacancies  are  available. 

The  the  model  also  explains  why  a  simple  error  function 
diffusion  profile  fit  the  results  so  well,  which  led  to  the  low 
"reported  diffusion  coefficients";  most  of  these  experiments 
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were  done  for  short  times  and  were  outdif fusions ,  and  the 
rather  complex  model  reduces  to  a  simple  error  function 
diffusion  model  for  these  conditions.  It  also  explains  why  for 
longer  times  the  simple  model  no  longer  applied;  the  proposed 
model  predicts  such  different  behavior  for  longer  times. 

The  fast  interstitials  explain  the  gettering  behavior  of 
Cr,  specifically,  how  the  Cr  move  from  one  place  to  another 
so  quickly  (especially  from  the  bulk  to  the  surface).  Magee  et 
al.  (123)  recently  studied  the  gettering  action  of  Cr,  in  this 
case  to  a  mechanically  damaged  surface.  They  measured  the  flux 
of  the  gettered  Cr  (in  atoms /cm^ /sec. )  versus  temperature, 
from  300  to  400 3C.  Since  the  flux  due  to  a  stress  field  equals 
FDC/kT  (from  ref.  94),  where  F  is  the  force  exerted  on  the  Cr 
by  the  stress  field,  D  is  the  effective  diffusion  coefficient, 
and  C  is  the  concentration,  which  is  constant  with  temperature 
in  this  case,  a  plot  of  ln(T»flux)  versus  1/T  results  in  the 
activation  energy  for  the  diffusion  coefficient.  From  Magee's 
data,  a  value  of  0.9  eV  was  determined.  The  activation  energy 
for  the  interstitial  controlled  diffusion,  the  case  that 
allegedly  is  responsible  for  the  gettering,  was  determined  in 
this  work  to  be  0.83  eV,  which  is  within  experimental  error  of 
Magee's  result.  If  the  flux  seen  in  Magee's  results  is  not  due 
to  a  stress  field,  but  due  to  concentration  gradient  diffusion 
of  the  interstitials,  to  a  region  of  lower  Cr  activity  at  the 
surface,  then  plotting  ln(flux)  vs  1/T  results  in  a  value  for 
Ea.  In  that  case,  Ea  =  .88  eV,  even  closer  to  Magee's 
value.  In  either  case,  the  results  are  good  evidence  that  the 
Cr  species  that  accounts  for  gettering  (as  seen  by  others  and 


in  this  work)  and  the  Cr  species  which  accounts  for  the  deep 
penetration  seen  in  this  work  {modelled  by  an  interstitial 
species)  are  the  same. 

The  low  surface  concentration  of  Cr  observed  when 
annealing  Cr-doped  GaAs  (without  an  arsenic  overpressure)  is 
explained  in  terms  of  the  ternary  phase  diagram,  as  discussed 
in  section  6.5.3,  i.e.,  by  losing  arsenic  from  the  GaAs  surface 
and  moving  to  the  Ga-rich  side  of  the  phase  diagram.  This  is 
consistent  with  the  results  of  Kasahara  et  al.  (112)  in  which 
the  surface  concentration  of  Cr  changed  from  less  than  10^6 
to  «5  x  1()16,  an(fi  then  to  =>7  x  10^-6  cm“3  (at  850°C),  as 
the  arsenic  overpressure  was  increased  from  0  to  0.15,  and  then 
to  7.6  torr.  In  the  7.6  torr  case,  virtually  all  the  arsenic 
for  the  vapor  phase  is  externally  supplied;  therefore,  none 
comes  from  the  GaAs  wafer,  and  the  system  remains  in  the 
GaAs ( s )-CrAs ( s ) -Ga ( ^ )  tie  triangle  region  where  the  Cr 
solubility  in  GaAs(S)  is  =»6  x  lO-'-6  cm"3 .  As  discussed  in 
section  6.2,  Kasahara  (112)  attributed  the  different  profiles 
to  different  diffusion  coefficients;  however,  as  mentioned,  the 
curves  can  be  fit  better  to  profiles  with  the  same  D  value 
(»3  x  10“12  cm^/sec.),  but  with  different  surface 
concentrations. 

Annealing  Cr-doped  GaAs  with  an  encapsulent,  as  reported 
by  others  and  discussed  in  section  6.2,  also  results  in  Cr 
profiles  with  a  low  concentration  of  Cr  at  the  surface.  As 
suggested  by  others,  arsenic  may  diffuse  into  and  through  the 
encapsulents ,  leaving  a  Ga-rich  phase  behind  at  the  surface  of 
the  GaAs  (40).  According  to  the  results  of  this  work,  this 
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would  result  in  a  low  solubility  of  Cr  in  the  GaAs,  which  would 
explain  the  low  levels  of  Cr  at  the  surface. 


The  indiffusion  results  of  Tuck  et  al.  (5,  120),  as 
discussed  in  section  6.2,  can  be  explained  in  either  of  two 
ways:  1)  the  Cr  at  the  surface  is  "built-up  Cr" ,  to  be 
discussed  in  the  next  chapter,  and  the  maximum  bulk  level  is 
lower  than  seen  in  this  work  because  the  GaAs  is  in  equilibrium 
with  a  Ga-rich  liquid  rather  than  CrAS(S)  +  Ga(z)  (because 
no  arsenic  source  was  used);  or  2)  the  Cr  at  the  surface  is  the 
"vacancy  diffusion  controlled  case"  Cr,  and  the  bulk  Cr  level 
doesn't  rise  as  high  because  the  vacancy  generation  mechanism 
is  not  as  efficient  in  their  case  (due  to  less  dislocations, 
perhaps,  in  the  GaAs).  The  fact  that  the  Cr  levels  were  even 
lower  when  elemental  arsenic  was  added  is  probably  not  due  to 
moving  the  system  from  the  Ga-rich  side  to  the  tie  triangle 
region  (which  would  predict  higher  levels  of  Cr),  but  is 
probably  due  to  the  formation  of  other  Cr-As  compounds  from 
elemental  Cr  and  As,  which  probably  changes  the  chemical 
potential  of  Cr. 

The  phenomenon  of  the  high  build-up  of  Cr  at  the  surface, 
seen  in  this  and  in  others'  work  will  be  discussed  in  chapter 


6.8  Summary  of  Chromium  Diffusion  Studies 


In  this  chapter  the  studies  of  the  diffusion  of  Cr  in  GaAs 
were  presented.  The  diffusion  of  Cr  in  GaAs  was  studied  with 
well-defined  experiments,  that  used  relevant  phase  equilbria, 
and  using  well-defined  solubility  information,  all  of  which  was 


determined  in  this  research.  The  diffusion  behavior  was  shown 
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to  follow  the  substitutional/interstitial  dissociative 
diffusion  model,  uncharged  case.  Values  for  key  parameters  of 
this  model,  including  the  effective  diffusion  coefficients  for 
vacancy  diffusion  controlled  Cr  diffusion  and  interstitial 
controlled  diffusion,  and  the  equilibrium  time  constant  (which 
controlled  the  rate  at  which  the  latter  type  of  diffusion 
became  dominant  over  the  former  type),  were  determined  for  a 
large  section  of  the  Ga-As-Cr  ternary  phase  diagram  from  700  to 
1000 °C.  This  diffusion  model,  and  the  phase  equilibria  and 
solubility  information,  explain  much  of  the  observed  phenomena 
concerning  Cr  in  GaAs  during  heat  treatments,  including:  deep 
penetration  of  Cr  into  GaAs;  the  near-uniform  increases  or 
decreases  of  the  bulk  level  of  Cr  in  GaAs;  the  two  different 
types  of  diffusion  -  with  different  effective  diffusion 
coefficients  -  seen  in  various  situations;  changes  in  the 
diffusion  profiles  due  to  a  change  in  surface  concentration  of 
Cr  (and  not  to  a  change  in  diffusion  coefficient,  as  previously 
assumed);  and  the  gettering  behavior  of  Cr  in  GaAs,  especially 
from  a  Cr  doped  substrate  to  the  surface  of  the  LPE  layer  many 
microns  away.  Also,  it  was  shown  that  the  Frank-Turnbull ,  or 
vacancy  model,  for  the  interstitial-to-substitutional 
conversion  is  probably  operating  rather  than  the  kickout 
mechanism.  In  addition,  it  was  demonstrated  how  chromium  can 
be  used  as  a  probe  in  GaAs,  due  to  its  extreme  sensitivity  to 
its  environment  in  regards  to  its  diffusion  behavior  and 
solubility,  to  study  the  defect  chemistry  of  GaAs. 


7.  SURFACE  PHENOMENA 


7 . 1  Introduction 

Of  great  importance  in  the  study  of  the  behavior  of 
chromium  in  gallium  arsenide  is  the  subject  of  surface 
phenomena.  There  have  been  two  types  of  unusual  behavior 
exhibited  by  chromium  in  the  surface  region.  One  is  the  large 
surface  build-up,  which  is  most  evident  after  outdif fusion.  In 
normal  outdif fusion ,  the  concentration  level  of  the  impurity 
decreases  to  a  low  value  at  the  surface;  however,  with  Cr  in 
GaAs ,  the  concentration  level  often  increases  to  a  higher  value 
right  at  the  surface,  as  discussed  in  chapter  6.  As  was 
stated,  many  causes  for  this  have  been  proposed,  including: 
stress  due  to  encapsulents ,  a  high  dislocation  density  at  the 
surface,  and  arsenic  evaporation  or  outdif fusion  leading  to  the 
formation  of  gallium  droplets.  This  high  and  erratic  build-up 
was  seen  in  my  studies,  both  after  outdiffusion  and 
indiffusion.  In  addition,  the  observation  of  the  Cr  level 
above  the  solubility  level  (for  the  specific  conditions  of  the 
system)  was  sometimes  observed  in  the  LPE  layers  grown  for  the 
junction  diffusion. 

The  other  surface  phenomena  observed  was  the  depletion  of 
Cr  near  the  surface  after  indiffusion.  A  depletion  near  the 
surface  after  outdiffusion  is  normal;  a  depletion  after 
indiffusion  is  quite  unexpected  and  suggests  uphill  diffusion. 

In  this  chapter,  the  results  of  my  studies  on  these 
phenomena  are  presented.  First,  what  exactly  was  observed 
is  presented,  followed  by  a  discussion  on  several  topics 
concerning  these  phenomena:  SIMS  analysis  effects?  phase 


diagram  considerations;  and  solubility  effects,  which  might  be 


of  help  in  determining  their  cause  or  causes.  Next,  the 
experimental  testing  of  some  of  the  theories  is  discussed,  as 
well  as  the  probable  causes  based  on  the  results.  Next,  the 
implications  of  these  results  on  the  behavior  of  Cr  in  GaAs  is 
discussed,  especially  regarding  chromium  outdif fusion.  In 
addition,  the  Mgh  Cr  levels  in  some  of  the  LPE  layers  is 
explained  based  on  the  proposed  theories. 

7.2  Observations  Reaardina  Cr  Build-uc  and  Depletion 


Shown  in  figure  7-1  are  two  SIMS  profiles  of  Cr  m  GaAs 
following  800°C,  12  hour  Cr  indiffusion  runs  (using  a 
well-defined  Cr  source).  In  these  profiles,  as  in  virtually 
all  the  indiffusion  profiles  obtained  initially  in  this 
research,  two  characteristics  are  evident:  1)  a  high  and 
erratic  build-up  of  Cr  at  the  surface,  up  to  1000  times  higher 
than  the  bulk  solubility  level;  and  2)  a  depletion  of  Cr  1-4 
microns  deep  for  700-900 °C  anneals.  Both  of  these  phenomena 
increase  with  annealing  temperature  and  time.  In  addition, 
these  two  phenomena  appear  to  be  totally  unrelated;  there  is  no 
correlation  between  the  occurrence  and  magnitude  of  them. 
However,  there  does  seem  to  be  a  correlation  between  the  Cr 
build-up  and  the  occurence  of  certain  surface  features  observed 
by  optical  microscopy;  these  are  small  (10-200  microns  across) 
solidified  droplets  of  a  silver  appearing  material  sitting  in 
rectangular  recesses.  Examples  are  shown  in  figures  7-2  and 
7-3  (after  12  hour,  900°C  anneals).  Features  such  as  these 
have  been  reported  quite  often  and  are  believed  to  be  gallium 


droplets  which  form  when  arsenic  preferentially  vaporizes  from 
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CONCENTRATION  ( 


tne  surface  (108,  133). 

The  Cr  surface  build-up  is  also  observed  after 
indiffusion.  Sometimes  this  is  over  a  very  narrow  region  (0.01 
to  0.04  microns),  but  other  times  it  is  much  deeper  (up  to  0.5 
microns ) . 

7 . 3  Discussion  of  Possible  Causes  of  Cr  Anomalies 
7.3.1  SIMS  Eftacts 


One  must  consider  the  possibility  that  the  Cr  build-up 
and/or  depletion  seen  in  the  profiles  may  by  due  to  the 
profiling  process  itself  -  SIMS.  In  chapter  4,  several  SIMS 
phenomena  were  discussed,  including  enhanced  secondary  ion 
yield  due  to  the  presence  of  oxides,  and  the  "artificial  dip" 
near  the  surface  related  to  the  projected  range  of  the  primary 
ion.  As  discussd  previously,  these  effects  are  far  too  small 
to  explain  the  depletion  and  large  build-up  seen.  However,  the 
enhanced  yield  due  to  an  oxide  layer  (which  usually  is  present 
on  the  surface  of  GaAs )  may  be  the  cause  of  the  very  high,  yet 
extremely  shallow  (-0.01  -  0.04  microns  deep)  build-up  of  Cr 
usually  present  in  the  SIMS  profiles.  This  is  present  even  in 
the  Cr  profiles  of  unannealed  GaAs  (see  figure  7-4),  leading 
one  to  suspect  that  it  is  a  SIMS  artifact. 

The  other  possibility  for  the  large  build-ups  and/or 
depletions  is  that  uneven  sputtering  near  the  surface  could 
result  in  the  anomalous  profiles.  This  was  easily  checked  by 
profiling  the  arsenic  and  gallium  concentrations.  If  the 
sputtering  rate  changes,  then  the  As  and  Ga  concentrations 
would  change.  Figure  7-5  shows  the  As  and  Ga  profiles,  as  well 
as  the  Cr  profile,  for  an  indiffusion  sample  (12  hr.  @  800°C). 
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Figure  7-4.  Cr  concentration  profile  from  unannealed,  C 
GaAs. 
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Figure  7-5.  Cr,  Ga,  and  As  ion  intensity  profiles  from  800 
12  hr.  indiffusion  experiment. 
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As  one  can  see,  the  As  and  Ga  levels  are  uniform,  indicating 
that  the  Cr  depletion  is  not  a  SIMS  artifact.  In  figure  7-6, 
the  arsenic  profile  is  shown  with  the  Cr  profile  for  a 
different  sample.  In  this  case,  the  arsenic  level  is  slightly 
lower  in  the  region  where  the  Cr  build-up  occurs.  However,  the 
fact  that  the  percent  change  in  the  As  level  is  nowhere  close 
to  the  percent  change  in  the  Cr  level  means  that  the  Cr 
build-up  is  not  merely  due  to  a  change  in  sputtering  rate; 
rather,  whatever  is  causing  the  change  in  arsenic  rate  is 
probably  also  changing  the  Cr  level.  This  will  be  made  more 
clear  later.  In  conclusion  then,  the  phenomena  in  question  do 
not  appear  to  be  due  to  anomalous  profiling  by  SIMS.  However, 
a  clue  has  been  given  regarding  the  build-up  of  Cr  at  the 
surface,  i.e.,  a  slight  depletion  of  arsenic  there. 

7.3.2  Phase  Diagram  Considerations 

As  discussed  in  chapter  3,  one  of  the  results  of  the  phase 
diagram  determinations  is  that  a  Ga-rich  liquid  is  in 
equilibrium  with  GaAs  for  the  temperatures  of  interest.  It  is, 
in  fact,  one  of  the  phases  in  the  source  material  and  one  of 
the  phases  present  in  the  final,  equilibrium  state  of  the 
solubility  or  diffusion  system.  In  addition,  the  Ga(£)  phase 
has  a  high  Cr  solubility  (=>10  a/o  at  900°C),  and  the 
distribution  coefficient  between  Cr  in  GaAs ( s )  and  in  the 
Ga(ij  phase  is  quite  low  (in  the  10-3  range).  Therefore, 
if  some  Ga^)  should  be  present  in  the  system,  it  would 
dissolve  a  cons iderab1 n  amount  of  Cr,  and  at  a  much  higher 
concentration  than  in  GaAs(s).  If,  for  some  reason,  the 
Ga(Z)  phase  should  be  present  on  the  surface  of  the  GaAs 
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wafer  during  a  diffusion  or  solubility  anneal,  then  as 
equilibrium  is  approached,  more  and  more  Cr  would  be  dissolved 
in  the  liquid  phase  until  the  high  equilibrium  liquid 
solubility  level  is  reached. 

The  formation  of  Ga  droplets  on  the  surface  of  GaAs,  as 
previously  discussed,  is  commonly  observed.  As  discussed  in 
chapter  2  regarding  reaching  the  final  equilibrium  state, 
achieving  the  equilibrium  vapor  pressure  of  all  the  components 
must  be  taken  into  account.  In  the  case  of  annealing  pure  GaAs 
in  an  ampoule  at  900 °C,  the  equilibrium  partial  pressure  of 
arsenic  is  =>10“^  atm  (134),  while  that  of  gallium  is  «*10-7 
atm.  This  explains  the  formation  of  Ga  pools  on  GaAs;  much 
more  arsenic  leaves  the  surface  than  Ga  in  order  to  achieve  the 
equilibrium  vapor  pressure  and,  therefore,  regions  of  pure  Ga 
are  left  behind  (135).  This  behavior  is  illustrated  in  figure 
7-7  which  shows  the  surface  of  a  GaAs  wafer  which  was  annealed 
in  an  ampoule  at  1000°C  for  2  hours.  Also  shown,  in  figure 
7-8,  is  the  surface  of  a  GaAs  wafer  annealed  with  enough 
elemental  arsenic  added  to  the  ampoule  so  that  the  partial 
pressure  of  arsenic  is  =<1  atm.,  greater  than  the  equilibrium 
pressure  of  arsenic  over  GaAs,  and  hence  large  enough  to 
prevent  dissociation  of  the  GaAs  wafer.  (Annealing  in  an 
arsenic  overpressure  is  common  practice  in  GaAs  processing 
(13  6).  ) 

It  is  a  reasonable  to  assume  that  the  partial  pressure  of 
arsenic  in  the  three  condensed-phase  tie  triangle  region 
containing  GaAs(s),  CrAs(s),  and  Ga^j  is  high  relative 
to  the  partial  pressures  of  Cr  and  Ga.  Therefore,  one  would 
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igure  7-7.  Optical  micrograph  (5Q0X)  of  GaAs  surface  after 
1000°C,  2  hr.  anneal. 


Figure  7-8.  Optical  micrograph  (500X)  of  GaAs  surface  after 

1000°C,  2  hr.  anneal  -  with  elemental  arsenic  added 
to  ampoule. 
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expect  that  more  arsenic  than  gallium  would  leave  the  starting 
material  (both  the  GaAs  wafer  and  the  source  material)  in  order 
to  achieve  the  equilibrium  vapor  pressure.  This  would  leave  Ga 
droplets  on  the  GaAs  surface,  which  would  then  dissolve  Cr  to 
its  solubility  level. 

This  process  is  shown  schematically  using  a  "reaction 
path"  diagram  in  figure  7-9.  Initially,  the  GaAs  surface  is  at 
point  A.  Upon  heating,  arsenic  selectively  vaporizes  from  the 
surface,  leaving  Ga  droplets.  The  composition  of  these 
droplets  is  represented  by  point  B.  However,  in  attaining 
equilibrium,  the  composition  of  the  droplets  moves  towards 
point  C  with  an  increase  in  Cr  content.  Therefore,  the  total 
system  at  equilibrium  consists  of  GaAs(S),  CrAs(S),  Ga^) 

(at  point  C),  and  the  vapor  phase,  with  some  Ga^)  on  the 
GaAs(S)  surface,  as  well  as  in  the  source. 

Because  of  the  high  Cr  solubility  in  the  Ga^)  phase 
compared  to  in  GaAs  ( s )  ,  even  a  small  amount  of  the  Ga^j 
phase  on  the  surface  could  lead  to  a  large  Cr  level  as  measured 
by  SIMS.  Even  if  only  1/1000  of  the  area  analyzed  by  SIMS  was 
Ga(£),  the  Cr  level  could  increase  by  100  times  just  because 
of  the  relative  solubilities.  Because  of  these  considerations, 
and  because  of  the  correlation  between  the  occurrence  of  the  Ga 
droplets  and  the  Cr  build-up,  it  was  suspected  that  this  was 
the  main  cause  for  the  Cr  build-up  at  the  surface  seen  in  these 
studies. 

7.3.3  Discussion  on  Solubility  in  Regards  to  Cr  Depletion 

Another  consideration  regarding  solubility  is  the  effect 
that  other  impurities  have  on  the  solubility  of  chromium  in 
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Figure  7-9.  Schematic  diagram  of  reaction  path  of  GaAs  surface 
regions. 
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GaAs.  As  discussed  in  chapter  4,  the  presence  of  other 
acceptors  in  GaAs  would  decrease  the  solubility  of  an  acceptor, 
such  as  Cr,  and  the  presence  of  donors  would  increase  the 
solubility  of  an  acceptor.  This  latter  effect  was  shown  to  be 
true  for  the  case  of  chromium  in  Si  and  Te  doped  GaAs. 
Therefore,  a  build-up  of  Cr  in  GaAs  above  the  normal  case 
(i.e.,  in  undoped  GaAs)  might  be  due  to  the  presence  of  a  donor 
impurity;  likewise,  a  depletion  could  be  due  to  the  presence  of 
an  acceptor. 

As  discussed  in  chapter  4,  chromium  in  undoped  GaAs  is 
probably  uncharged  and  only  becomes  an  acceptor  when  n-type 
impurites  are  introduced.  Therefore,  it  might  not  act  as  a 
normal  acceptor  in  the  presence  of  another  acceptor  and 
decrease  in  solubility  due  to  electronic  effects.  However, 
since  most  acceptors  occupy  the  same  site  as  Cr  does  (Ga 
sites),  then  site  (or  vacancy)  competition  could  likewise  lead 
to  a  decrease  in  the  amount  of  each.  This  would  especially  be 
true  if  there  is  a  kinetic  problem  with  the  impurity  atoms 
incorporating  into  the  site,  such  as  is  the  case  for  Cr  and 
other  impurities  incorporating  into  a  Ga  site  in  GaAs  by 
finding  a  vacancy.  The  decrease  in  the  Cr  level  in  GaAs  due  to 
the  presence  of  other  impurity  atoms  has  been  reported  by 
Asbeck  (82),  who  observed  a  decrease  in  the  Cr  level  when  Zn 
and  Cu  atoms  were  introduced  into  the  GaAs. 

Another  cause  of  increased  or  decreased  solubility  is  an 
appreciable  change  in  defect  density  and/or  structure.  A 
region  that  has  an  abnormally  high  concentration  of 
dislocations,  vacancies,  stress  fields,  etc.,  could  change 
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the  basic  structure  of  the  crystal  enough  so  that  the 
"solubility"  of  an  impurity  in  it  could  be  different  than  that 
of  the  crystal  in  the  equilibrium  state.  By  definition,  the 
true  solubility  of  an  impurity  in  a  substance  is  the  solubility 
when  the  substance  is  in  the  equilibrium  state  (with  the 
equilibrium  density  of  vacancies,  etc.).  However,  depending  on 
the  history  of  the  material,  higher  or  lower  numbers  of  these 
defects  may  be  frozen  in.  There  are  three  ways  to  check  this: 
perform  the  solubility  anneal  for  a  "long"  time  to  see  if  the 
effect  decreases  or  not;  pre-anneal  the  substrate  at  a  high 
temperature;  or  remove  regions  such  as  surfaces  with  polishing 
or  slicing  damage. 

The  last  two  effects,  the  site  competition  and  defect 
effects,  are  kinetic  effects  and  will  therefore  disappear  in 
time.  In  contrast,  the  electrical  effect  applies  even  for 
equilibrium  conditions.  Another  way  one  impurity  can  change 
the  solubility  of  another  is  by  producing  stress  fields  in  the 
matrix  material  which  in  turn  influence  the  solubility  of  the 
other  impurity.  This,  like  the  electrical  effect,  is  an 
"equilibrium  effect"  and,  therefore,  will  not  anneal  out. 

An  additional  phenomena  is  that  of  impurity  segregation  to 
the  substrate  surface,  which  is  occasionally  observed.  It  is 
believed  to  occur  when  the  impurity  lowers  the  surface  tension 
and  hence  minimizes  the  energy  of  the  system  (137).  Therefore, 
for  some  impurities,  its  "solubility"  is  higher  at  the  surface 
than  in  the  bulk,  whereas,  some  impurities  have  lower 
solubilities  at  the  surface  because  they  raise  the  surface 
tension.  However,  either  case  would  result  in  reproducible 


levels  at  the  surface  (for  a  given  temperature),  which  is 
clearly  not  the  case  for  Cr  build  up  in  GaAs.  Also,  a  lowering 
of  the  surface  concentration  usually  occurs  for  cases  in  which 
the  impurity  increases  the  melting  point  of  the  substrate 
(137),  which  is  not  the  case  for  Cr  in  GaAs;  therefore  the 
depletion  is  probably  not  a  direct  result  of  this. 

7 . 4  Experimental  Testing  To  Find  Causes  of  Phenomena 
7.4.1  Chromium  Build-up 

First  to  be  considered  is  the  large  chromium  build-up  seen 
at  the  surface  following  Cr  solubility  and  diffusion  anneals. 

As  mentioned,  a  likely  explanation  for  this  is  the  formation  of 
gallium  droplets  which  have  a  high  Cr  solubility.  Two 
observations,  previously  discussed,  are  consistent  with  this 
theory:  1)  there  is  a  correlation  between  the  number  and  size 
of  the  droplets  and  the  magnitude  of  the  Cr  build-up  as 
measured  by  SIMS;  and  2)  the  build-up  increases  with  increasing 
annealing  time  and  temperature,  both  in  concentration  level  and 
in  depth  into  the  GaAs. 

To  test  the  theory,  two  indiffusion  runs  were  done  at 
800°C,  for  12  hours,  using  H  source.  Sample  E  was  in  an 
ampoule  with  a  small  volume  ( =  2  cm^  total),  with  a  large 
amount  of  source  material,  and  the  GaAs  wafer  was  large. 

Sample  L  was  a  small  GaAs  wafer,  had  a  small  amount  of  source 
material,  and  was  in  a  large  ampoule  (=<6  cm^  ) .  If  the 
build-up  is  due  to  Cr  in  Ga^),  which  forms  when  arsenic 
leaves  the  GaAs  surface  to  provide  for  the  equilibrium  arsenic 
vapor  pressure  in  the  ampoule,  then  sample  E  should  result  in 


less  Cr  build-up.  The  reason  for  this  is  that  less  Ga/ } \ 


will  form  on  the  surface  because  less  arsenic  is  needed  in  the 
ampoule  due  to  the  smaller  volume.  Also,  because  of  the  larger 
size  of  wafer  and  larger  amount  of  source  material,  less 
arsenic  will  leave  per  square  cm.  of  surface.  The  resulting 
SIMS  Cr  profiles  are  shown  in  figure  7-10.  As  predicted, 
sample  E  shows  a  much  smaller  Cr  build-up.  Shown  in  figures 
7-11  and  7-12  are  micrographs  of  the  two  surfaces.  Again  as 
predicted,  the  surface  of  sample  E  has  much  less  Ga/^ 
present. 

Another  experiment  was  done  to  confirm  the  fact  that  the 
Ga  droplets  on  the  GaAs  surface  do  contain  large  amounts  of 
Cr.  A  GaAs  wafer  was  annealed  in  an  ampoule,  with  H  source,  at 
1000 °C  for  12  hours,  and  quenched.  Many  Ga  droplets  were 
observed.  A  typical  one  was  analyzed  by  an  electron  microprobe 
(sampling  size  of  -1  urn2).  The  results  were:  68.7  a/o  Ga; 

14.9  a/o  As;  and  16.4  a/o  Cr.  According  to  the  1000 °C  ternary 
isotherm,  these  numbers  are  within  5%  of  what  the  Ga ( ^  ) 
composition  should  be.  In  addition,  a  spot  next  to  the  Ga 
droplet  (in  the  GaAs  matrix)  was  analyzed  to  check  the  accuracy 
of  the  analysis.  The  results  were:  50.01  a/o  Ga;  49.99  a/^ 

As,  and  0.000  a/o  Cr,  i.e.,  below  the  detection  limit  for  Cr; 
these  values  confirmed  the  accuracy  of  the  analysis. 

Finally,  a  GaAs  wafer  was  annealed,  again  at  1000°C  for  12 
hours,  with  H  source,  but  this  time  slow  cooled  to  room 
temperature.  According  to  the  phase  equilibria  studies,  the 
Ga(Z)  phase,  upon  cooling,  should  convert  to  CrGa^g)  and 
Ga(4)  (with  less  Cr),  plus  a  little  CrAS(S)  and  GaAS(S). 

The  sample  was  crushed  and  analyzed  by  XRD.  The  diffraction 
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Figure  7-11.  Optical  micrograph  (500X)  of  GaAs  surface  of 
experiment  L. 


Figure  7-12.  Optical  micrograph  (500X)  of  GaAs  surface  of 
experiment  E. 


pattern  indicated  that  there  was  some  CrGa4(s)  present. 

One  final  observation  was  that  evidence,  or  "artifacts", 
of  Ga  droplets  could  be  seen  in  the  interior  of  the  SIMS  crater 
in  many  cases  when  there  was  Cr  build-up  at  the  surface.  The 
actual  Ga  droplet  was  gone  by  the  time  the  entire  SIMS  profile 
was  completed  since  it  was  only  in  the  surface  region;  however, 
small  areas  of  uneven  morphology  resulting  from  slightly  uneven 
sputtering  of  that  phase  could  be  seen. 

The  SIMS  arsenic  profile,  which  shows  the  the  As  level  to 
be  slightly  less  in  the  region  where  the  Cr  build-up  occurs, 
can  now  be  explained:  the  sputtering  rate  may  be  a  little 
different  there  due  to  the  second  phase  present;  in  addition, 
that  phase  (Ga^j)  has  less  arsenic  than  the  matrix. 

In  conclusion,  the  major  cause  of  the  Cr  surface  build-up 
is  the  formation  of  the  Ga-rich  liquid  phase  which  forms  when 
arsenic  selectively  vaporizes  in  order  to  establish  the 
equilibrium  arsenic  vapor  pressure.  Even  with  small  volume 
ampoules,  the  surface  build-up  was  not  always  eliminated, 
especially  during  the  higher  temperature  anneals.  However,  it 
was  believed  to  have  had  a  negligible  influence  on  the  actual 
diffusion  process. 

7.4.2  Chromium  Depletion 

In  contrast  to  the  chromium  build-up  phenomenon,  the  cause 
of  the  Cr  depletion,  after  indiffusion,  was  not  quite  so  easy 
to  determine.  However,  after  many  experiments,  and  a  little 
luck,  a  cause  was  discovered. 

Initially,  a  few  theories  of  its  origin  were  considered, 
but  none  of  these  were  appropriate.  One  of  these  theories  was 


the  uneven  SIMS  sputtering  rate.  3ut,  as  previously  discussec, 
uniform  gallium  and  arsenic  concentration  levels  in  that  region 
disproved  that.  Another  theory  proposed  that  Cr  outdiff used 
during  the  quencn,  because  the  Cr  solubility  is  lower  at 
lower  temperatures.  This  was  the  reason  given  for  similar 
profiles  for  Au  in  InAs  (62)  and  Cu  in  I.nAs  (138).  Tuck  et 
al.  (59),  who  made  the  same  observation  with  Au  in  InP,  did  not 
accept  this  theory  for  all  those  systems,  saying  that  the 
results  were  not  consistent  with  one  another.  This  theory  was 
tested  with  Cr  in  GaAs  by  me  by  changing  the  quench  rate. 
However,  no  change  was  seen  in  the  profiles.  Besides,  the 
diffusion  coefficients  for  these  times  and  temperatures  are  far 
too  low  to  account  for  that  much  depletion. 

It  was  therefore  concluded  that  the  Cr  depletion  seen  is 
not  an  artifact  of  the  profiling,  but  is  real  and  occurs  during 
the  indiffusion.  This  suggests  "uphill  diffusion",  diffusion 
up  the  concentration  gradient.  However,  in  light  of  the 
proposed  diffusion  mechanism,  this  type  of  behavior  can  be 
explained,  as  discussed  in  section  5.3  of  chapter  5,  without 
using  normal  uphill  diffusion  arguments.  Cr  atoms  diffuse  in 
interstitially ,  down  the  interstitial  concentration  gradient. 
They  then  convert  to  substi tutional  atoms,  and  more 
interstitials  diffuse  in.  The  concentration  of  the 
substitutional  atoms  builds  up  to  the  substitutional  solubility 
level,  whatever  it  is  in  that  region.  If  a  region,  such  as  a 
region  near  the  surface,  has  a  lower  Cr  solubility  than 
adjacent  regions,  the  final  Cr  concentration  will  be  lower; 
this  situation  could  lead  to  profiles  such  as  we  have  seen. 
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Therefore,  it  was  presumed  that  the  cause  of  the  Cr  depletion 
near  the  surface  was  a  lower  Cr  solubility  near  the  surface. 

The  cause  of  the  lower  solubility  was  not  easy  to  determine. 

Initially  it  was  presumed  that  surface  damage  and/or 
stress  at  the  surface,  caused  by  slicing  and  polishing, 
might  cause  the  lower  Cr  solubility.  The  fact  that  the 
depletion  did  not  occur  after  very  high  temperature  anneals  (_> 
1000°C),  except  for  very  short  times,  suggests  that  the  damage 
was  annealed  out  at  the  high  temperatures.  To  test  this 
theory,  three  experiments  were  conducted.  One  was 
pre-annealing  the  wafer  before  the  actual  indiffusion  to  reduce 
the  damage  and/or  stresses.  This  was  done  at  1000°C  for  6 
hours,  with  either  elemental  arsenic  in  the  ampoule  to  prevent 
arsenic  vaporization,  or  with  J  source  (which  would  introduce 
only  a  very  low  level  of  Cr  into  the  GaAs,  and  yet  provide  a 
source  of  arsenic).  In  each  case,  after  doing  subsequent 
normal  indiffusions  at  800°C  for  12  hours  with  H  source,  the 
depletion  region  was  still  present  and,  if  anything,  worse. 
Another  experiment  was  pre-etching  the  wafer  to  physically 
remove  the  surface  region.  This  was  done  using  various  etching 
solutions,  e.g.,  16:1:1  H20: H 20 2 :H 2S0 4 ,  and  etching  off  from  1 

to  20  ums.  Subsequent  indiffusion  again  resulted  with 
depletion  regions.  The  third  experiment  was  simply  to  do  the 
indiffusion  for  long  times  and  see  if  the  depletion  is  reduced 
with  time.  However,  after  12,  18,  and  14  hours  at  900°C,  there 
was  still  a  depletion  in  each  case,  and  it  actually  increased 
in  depth  with  time.  From  these  experiments,  it  was  concluded 
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that  a  damaged  or  "non-equilibrium"  surface  region  was  not  the 
cause  of  the  Cr  depletion  there. 

Another  cause  for  a  change  in  solubility  or  expected 
concentration  level  locally,  is  the  presence  of  other 
impurities,  as  discussed  previously.  It  was  therefore 
considered  that  an  acceptor  or  an  impurity  on  a  Ga  site  might 
be  present  in  the  surface  region.  Using  SIMS,  the  various 
common  impurities  -  Si,  Cu,  Fe,  0,  and  C  -  were  searched  for. 
However,  none  showed  high  levels  near  the  surface.  It  was  then 
discovered  that  Klein  and  Nordquist  (139)  proposed  that  p-type 
surface  conversion  occurred  in  undoped  GaAs,  after  annealing  in 
a  H2  atmosphere,  due  to  high  levels  of  manganese  at  the 
surface.  The  Mn  levels  reached  the  mid  10^7  cm-3  range 
over  the  first  four  microns  after  a  20  minute  740°C  anneal. 
Virtually  all  of  it  was  electrically  active.  Klein  concluded 
that  the  Mn  was  gettered  to  the  surface  from  the  bulk,  of  the 
GaAs  (which  normally  has  levels  of  10^“*  to  lO^-3  cm~3  )  and 
didn't  come  from  outside  the  wafer.  Others  (41,  127,  140,  141) 
have  also  recently  studied  this  phenomena.  Palmateer  et  al. 
(127),  like  Klein,  found  that  the  Mn  accumulation  decreased 
after  successive  anneal/polishing  cycles,  consistent  with  the 
theory  that  the  Mn  comes  from  inside  the  GaAs.  Klein  suggested 
that  it  is  the  indiffusion  of  gallium  vacancies  that  causes 
this.  However,  when  annealing  in  an  argon  atomosphere  rather 
than  in  hydrogen,  Klein  found  that  no  Mn  accumulated  at  the 
surface,  the  reason  for  which  is  not  known.  It  appears  that 
there  is  a  complex  relationship  between  the  atmosphere  (and/or 
surface)  and  the  Mn  accumulation  phenomena.  For  whatever 
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reason  the  Mn  does  accumulate  at  the  surface,  if  it  _is  coming 
from  the  bulk,  it  strongly  suggests  that  the 
substitutional/interstitial  diffusion  mechanism  is  again  at 


work. 


In  any  case,  manganese,  a  known  acceptor  in  GaAs,  was 
analyzed  by  SIMS.  In  many  indiffusion  samples,  Cr  and  Mn  were 
simultaneously  profiled.  Three  examples  are  shown  in  figures 
7-13,  and  7-14.  As  one  can  see,  the  Cr  depletion  and  Mn 
build-up  complement  each  other  in  each  case.  This  was  true  for 
all  the  cases  that  Mn  was  profiled  (15  in  all).  In  addition, 

Cr  and  Mn  were  profiled  in  several  of  the  LPE  junction  samples 
after  various  anneals,  with  and  without  J  sources  in  the 
ampoules.  In  each  of  these  there  was  Mn  build-up  at  the 
surface,  but  not  in  the  interface  (junction)  region.  This 
explains  why  "normal"  Cr  diffusion  in  the  junction  region  can 
be  observed  (the  vacancy  diffusion  controlled  Cr  diffusion), 
whereas  this  cannot  usually  be  seen  at  the  surface  in  the 
indiffusion  samples.  From  all  these  observations  it  is 
concluded  that  the  decrease  in  Cr  concentration  near  the 
surface  after  indiffusion  is  due  to  the  build-up  of  Mn  there. 

As  stated  earlier,  it  appeared  that  the  depletion  did  not 
occur  at  higher  annealing  temperatures  (1000°C  and  above). 

The  SIMS  analysis  was  done  on  a  sample  annealed  at  1000 °C, 
again  profiling  both  Cr  and  Mn.  The  results  are  shown  in 
figure  7-15.  As  one  can  see,  the  Mn  is  still  there;  however, 
the  Cr  profile  does  not  appear  to  be  affected  by  the  presence 
of  the  Mn.  This  is  explainable  if  the  lowering  of  the  Cr  level 
is  due  to  an  electrical  effect.  At  900°C,  the  intrinsic 


carrier  concentration  is  =2  x  10^-7  cm-^  (13),  while  at 
1000°C  it  is  =5  x  1017  cm-^  .  Therefore,  as  the  temperature 
is  increased,  the  electrical  effect  of  the  Mn  is  swamped  out  by 
the  intrinsic  carrier  concentration,  and  the  Mn  has  little  or 
no  effect  on  the  Cr;  the  presence  of  Mn  would  not  change  the 
hole  or  electron  concentration  in  the  GaAs.  Also,  in  the  two 
Te-doped  GaAs  (n-type)  indiffusion  samples,  only  small 
depletions  are  observed.  This  is  consistent  with  the 
electrical  effect  theory,  since  the  electrons  from  the  Te  would 
compensate  the  holes  from  the  Mn. 

However,  one  can  not  rule  out  the  possibility  that  it  is 
not  entirely  an  electrical  effect  -  the  presence  of  the  Mn 
might  also  have  an  effect  on  the  Cr  level  just  because  it  is  on 
the  same  site  and  it  may  change  the  state  of  the  crystal  by 
mechanical  or  crystallographic  effects,  as  discussed  earlier. 
The  GaAs  might  only  be  able  to  accommodate  so  many  impurity 
atoms  on  the  gallium  sites,  whether  they  are  Mn  atoms  or  Cr 
atoms,  and  since  Mn  and  Cr  are  both  transition  elements  and 
right  next  to  each  other  in  the  periodic  table,  they  may  not  be 
very  distinguishable.  The  fact  the  Cr  level  is  affected  even 
when  the  Mn  level  is  below  the  n^  level  in  some  of  the 
profiles  would  suggest  that  it  is  not  entirely  an  electrical 
effect.  And  even  the  vacancy  or  site  competition  theory, 
closely  related  to  this,  is  not  impossible.  It  was  ruled  out 
because  it  was  believed  that  eventually  the  Cr  would  reach  its 
equilibrium  solubility  level  which  it  did  not  do  for  the  times 
measured.  However,  if  the  Mn  is  constantly  arriving  at  the 
surface,  then  the  Cr  level  might  be  kept  low  there  for  long 
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times.  The  observation  that  the  final  solubility  of  Cr  in 
Si-doped  (n-type)  GaAs  is  higher  than  in  undoped  GaAs,  but  that 
the  rise  to  the  solubility  level  is  slower  (Si  is  on  a  Ga  site) 
shows  that  both  electrical  and  "mechanical"  effects  are 
important.  At  this  point,  the  exact  reason  or  reasons  why  Mn 
causes  the  Cr  depletion  is  not  known;  however,  it  is  clear  that 
it  is  the  cause. 

An  attempt  was  made  to  get  rid  of  the  Mn  by  pre-annealing 
a  GaAs  wafer  in  0.1  atm  arsenic  overpressure  (using  elemental 
arsenic  in  the  ampoule)  for  24  hours  at  750 °C,  then  etching  the 
surface.  However,  the  arsenic  anneal  did  not  produce  any  Mn 
build-up  at  the  surface.  The  reason  for  this  is  not  known.  If 
the  Mn  builds  up  at  the  surface  by  gallium  vacancies  diffusing 
in,  one  would  expect  an  even  greater  tendency  for  this  to  occur 
in  an  arsenic  atmosphere  because  the  gallium  vacancy 
concentration  would  be  higher.  But  this  did  not  occur. 

7 . 5  Advantages  of  Having  Mn  Present 

Even  though  the  presence  of  manganese  at  the  surface  of 
the  GaAs  wafer  did  obscure  some  aspects  of  the  chromium 
diffusion  (especially  the  vacancy  diffusion  controlled  Cr 
diffusion  at  the  surface  for  indiffusions),  there  were  three 
benefits  of  having  it  there  and  affecting  the  chromium.  First, 
it  is  a  good  confirmation  of  the  substitutional/interstitial 
dissociate  diffusion  mechanism.  If  the  diffusion  was  the 
"normal"  type  (interstitial  and/or  substitutional  diffusion 
down  a  concentration  gradient),  and  the  surface  concentration 
was  lowered  (due  to  the  Mn),  then  that  would  be  the  maximum 
level  throughout.  However,  the  "uphill  diffusion"  seen  can 
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easily  be  explained  by  the  sub./int.  dissociative  mechanism,  as 
discussed  earlier.  And  as  Kendall  (77)  points  out,  this 
"uphill  diffusion"  is  a  good  indicator  of  this  type  of 
diffusion.  Secondly,  the  depletion  in  the  Cr  level  is  a  good 
indication  that  that  Cr  is  in  solution.  In  many  of  the  Cr 
profiles,  some  Cr  build-up  due  to  the  Ga  droplets  still 
occurred  to  some  extent,  and  it  might  have  been  difficult  to 
distinguish  the  "build-up  Cr"  from  the  Cr  in  solution. 

However,  the  Mn  only  affects  the  Cr  in  solution.  Therefore, 
the  "solution  Cr"  always  had  a  small  depletion  zone  near  the 
surface,  while  the  "build-up"  Cr  did  not,  and  the  two  were 
easily  distinguishable.  Thirdly,  because  the  Cr  level  in  the 
depletion  region  is  lower  than  the  bulk  level  (which  can  be 
quite  low  for  short  time/low  temperature  diffusion),  we  know 
that  the  level  in  the  bulk  (which  is  often  very  low  and  flat) 
is  above  the  background  level  of  the  SIMS. 

Therefore,  there  were  actually  some  benefits  to  having  the 
Mn  there  and  affecting  the  chromium.  Luckily,  neither  the  Mn 
or  the  Cr  build-up  had  any  apparent  affect  on  the  other  parts 
of  the  diffusion  process,  such  as  the  deep,  interstitial 
controlled  diffusion,  as  evidenced  by  how  well  the  model  fit 
the  profiles. 

7 . 6  Application  of  Surface  Phenomena  Results  to  Other  Cr 

Behavior 

The  results  of  these  surface  phenomena  studies  can  now  be 
used  to  explain  some  of  the  "anomalous"  behavior  of  Cr  in  GaAs 
seen  by  others  and  in  the  present  work.  For  outdif fusion ,  a 
high  Cr  build-up  at  the  surface  is  usually  observed,  as 
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mentioned  earlier.  Part  of  this  (0  to  .04  microns  deep)  can  be 
attributed  to  the  enhanced  yield  due  to  an  oxide  layer,  or  some 
other  SIMS  artifact.  However,  in  many  cases,  the  build-up  goes 
much  deeper  into  the  wafer,  and  this  can  be  explained,  in  part, 
by  gallium  droplet  formation.  As  was  the  case  with  the  Cr 
indiffusion,  this  build-up  could  be  reduced  by  using  a  smaller 
volume  ampoule  (see  figure  7-16). 

The  fact  that  others  have  shown  rather  conclusively  that 
the  presence  of  dislocations,  stress,  and  lattice  damage  (due 
to  ion  implantation,  for  example)  does  result  in  a 
redistribution  of  Cr  in  the  regions  where  the  defects  are 
present  (usually  with  increasing  Cr  level  concentrations  in 
these  regions)  suggests  that  these  factors  may  also  be 
responsible  for  some  of  the  Cr  build-up  at  the  surface. 

Perhaps  the  dislocations,  damage,  and/or  stresses,  or  whatever 
causes  them  ( encapsulents ,  etc.),  cause  the  formation  of 
gallium  droplets  or  inclusions,  which  then  results  in  the  high 
levels  of  Cr.  As  mentioned  earlier,  it  has  been  suggested  that 
arsenic  diffusion  into  or  through  the  encapsulents  may  result 
in  Ga  inclusions  at  the  GaAs  surface  (40,  142).  Gyulai  et  al. 
(143)  claim,  from  their  backscattering  measurements,  that 
gallium  does  accumulate  at  the  encapsulent/GaAs  interface  for 
both  Si02  and  Si3N4  encapsulents  following  700-800°C  anneals. 
Also,  Morkoc  (37)  reports  that  no  electrically  active  Cr  is 
detected  in  the  regions  of  high  Cr  build-up  at  the  surface. 

This  implies  that  the  Cr  is  not  in  solution  there,  consistent 
with  it  being  in  a  gallium  precipitate. 

An  extreme  example  of  Cr  surface  build-up  is  shown  in 
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figure  7-17.  This  is  a  result  of  annealing  Cr-doped  GaAs  (1 
hr.  3  800°C)  in  an  evacuated  ampoule  after  a  "different" 
cleaning  technique.  The  same  cleaning  schedule  as  outlined  in 
chapter  4  was  followed,  except  no  wafer  holder  was  used; 
rather,  the  GaAs  wafer  was  allowed  to  "slosh"  around  in  the 
bottom  of  the  cleaning  beaker.  As  one  can  see,  much  more  Cr 
build-up  occurs  than  when  the  GaAs  is  kept  in  the  wafer 
holder.  It  is  presumed  that  this  is  due  to  the  damage  that  is 
produced  on  the  surface  and/or  gallium  inclusions  that  form. 
This  illustrates  how  critical  the  handling  of  the  material  is 
and  how  sensitive  the  behavior  of  impurities  are  to  the  state 
of  the  crystal. 

It  is  not  clear  whether  or  not  the  presence  of  Mn  affects 
the  outdiffusion  of  Cr.  Mn  surface  accumulation  does  occur 
during  the  Cr  outdif fusions ,  as  evidenced  by  SIMS.  However,  a 
depletion  of  Cr  in  the  surface  region  is  expected  anyway  as  it 
is  a  normal  result  of  the  outdiffusion  as  discussed  earlier. 
The  manganese  might  be  expected  to  enhance  the  depletion  and 
result  in  higher  measured  effective  D  values.  The  fact  that 
most  of  the  Dsj_oy  values  for  the  outdiffusion  experiments  are 
slightly  higher  that  those  for  the  junction  diffusion 
experiments  (which  were  not  affected  by  Mn)  suggests  that  the 
Mn  is  affecting  the  outdiffusion  experiments. 

The  surface  build-up  of  Cr  in  the  junction  samples  after 
annealing  (see  figure  7-13  for  an  example  of  this)  can  also  be 
attributed  to  the  formation  of  gallium  droplets  and/or  damage. 
This  is  a  little  different  than  the  normal  outdiffusion  cases 


in  that  the  surface,  and  hence  Cr  build-up,  is  removed  from  the 
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Figure  7-17.  Cr  concentration  profile  from  800°C,  1  hr.  out 
diffusion  experiment,  using  different  cleaning 
technique. 


source  of  chromium,  sometimes  by  as  much  as  20  microns.  But 
the  substitutional/interstitial  dissociative  diffusion 
mechanism  can  easily  explain  this,  as  discussed  in  chapter  6; 
tht  fa  ,t  interstitial  Cr  species  can  quickly  diffuse  from  the 
substrate  to  the  surface.  (This  is  the  answer  to  the 
researchers  who  could  not  explain  how  Cr  got  to  the  surface  so 
fast  based  on  the  reported  diffusion  coefficients?  the 
"reported"  D  values  are  the  vacancy  diffusion  controlled  Cr  D 
values  and  only  apply  to  the  diffusion  seen  at  the  GaAs 
substrate/LPE  layer  interface. )  Once  the  Cr  interstials  get  to 
the  surface,  they  can  dissolve  in  the  gallium  phase, 
incorporate  into  the  stressed  or  damaged  surface  region,  or 
perhaps  become  subtitutional  by  reacting  with  excess  vacancies 
diffusing  in  from  the  surface.  Therefore,  why  the  chromium  can 
be  at  high  concentration  levels  at  the  surface  can  be  explained 
(dissolving  in  the  Ga  phase,  etc.)  as  well  as  how  it  can  get 
there  so  fast  (interstitial  diffusion). 

Another  use  of  the  results  of  this  section,  and  of  earlier 
chapters,  to  explain  other  observations  concerns  Cr  levels  in 
LPE  layers.  For  the  junction  diffusion  experiments,  as 
discussed  earlier,  undoped  LPE  GaAs  on  Cr-doped  GaAs  samples 
were  desired.  However,  the  first  two  samples  obtained  (from 
S.C.  Chang)  showed  relatively  high  levels  of  Cr  in  the  LPE 
layers,  as  measured  by  S IMS  (see  figures  7-19  and  7-20  for  the 
SIMS  profiles  of  such  samples).  In  these  samples,  the  meltback 
procedure  was  used  just  prior  to  growth  of  the  LPE  layer  in 
order  to  remove  any  residual  damage  and/or  impurities  at  the 
surface.  (This  is  done  by  slightly  raising  the  temperature 
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when  GaAs  (  s  )  and  Ga^)  (saturated  with  GaAs  (  s  )  )  are  first 
in  contact  with  each  other,  which  results  in  some  of  the  GaAs 
substrate  melting  and  becoming  part  of  the  melt. )  Because  of 
the  high  Cr  concentration  in  these  LPE  layers,  these  samples 
could  not  be  used  for  the  junction  diffusion  experiments. 
Samples  were  then  produced  without  the  meltback  step.  In 
these,  the  Cr  level  is  very  low  in  the  LPE  layer,  as  measured 
by  SIMS  (see  figure  7-21).  An  explanation  for  the  high  level 
of  Cr  in  the  LPE  in  the  first  samples  was  sought. 

Because  the  Cr  level  in  the  LPE  layer  is  fairly  uniform, 
it  was  initially  assumed  that  the  Cr  that  was  detected  was  in 
solution.  However,  based  on  the  solubility  results,  this  was 
not  possible.  For  a  GaAs  layer  to  be  grown  with  these  levels 
of  Cr,  it  must  be  in  equilibrium  with  a  gallium-rich  liquid  or 
with  ■  1-2  a/o  Cr  (based  on  distribution  coefficient  results). 
However,  the  gallium  melt  used  had  no  Cr  in  it  to  begin  with, 
and  diffusion  of  Cr  into  the  liquid  during  the  growth  procedure 
could  not  account  for  that  much  chromium  in  the  melt.  The  only 
explanation  was  that  the  Cr  was  not  in  solution  (i.e.,  not  in 
GaAs).  It  was  therefore  presumed  that  the  Cr  was  in  Ga (jj 
inclusions.  Optical  micrographs  were  taken  of  the  SIMS  craters 
(see  figures  7-22  ,  7-23,  7-24  ),  which  show  evidence  of  Ga^) 
inclusions  in  the  SIMS  craters  of  the  meltback  samples,  while 
the  craters  of  the  other  (low  Cr  level)  samples  show  no  such 
inclusions.  Electron  microprobe  analysis  revealed  that  the 
features  are  indeed  mostly  gallium  (91  a/o  Ga  and  9  a/o  As, 
which  is  exactly  the  composition  of  gallium,  saturated  with 
GaAS(s),  for  the  growth  temperature.  (No  chromium  was 
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Figure  7-24.  Optical  micrograph  (100X)  of  GaAs  surface  of  non 
meltback  sample,  showing  SIMS  crater. 


detected,  but  for  these  levels  in  the  LPE  layer,  none  would  be 
expected  based  on  the  detectability  of  the  analysis 
technique.)  Therefore,  it  was  concluded  that  during  the  LPE 
growth,  some  Ga^)  inclusions  remained  imbedded  in  the  GaAs 
layers,  with  substantial  Cr  dissolved  in  the  Ga (A)  which  was 
then  measured  in  the  SIMS  analysis.  To  test  this  theory,  a 
SIMS  Cr  analysis  was  done  on  another  area  of  one  of  the  samples 
which  showed  the  high  Cr  level.  If  the  Cr  was  in  solution  in 
the  GaAS(S),  then  all  SIMS  analysis  should  show  the  high  Cr 
level;  but  if  the  Cr  is  concentrated  in  the  few,  and  small,  Ga 
inclusions,  then  some  areas  (which  did  have  an  inclusion  in  it) 
should  show  low  levels  of  Cr  on  the  LPE  layer.  Figure  7-25  is 
such  a  SIMS  profile.  As  one  can  see,  the  Cr  level  in  the  LPE 
layer  is  very  low  (at  the  background  level  of  Cr  in  GaAs  by 
SIMS  analyses),  consistent  with  the  theory. 

The  results  suggest  that  perhaps  it  was  just  luck  that  the 
SIMS  analyses  of  the  two  meltback  samples  just  happened  to  hit 
areas  with  the  Ga  inclusion,  while  the  others  did  not. 

However,  it  does  appear  that  the  meltback  samples  do  have  a 
somewhat  larger  number  of  these  inclusions  than  the 
non-meltback  samples  (-30  cm”2  vs.  -  15  cm“^ ) .  The  reason 
why  the  Ga^)  inclusions  would  be  present  in  larger  numbers 
after  the  meltback  procedure  is  not  known;  however,  others  have 
observed  this  same  phenomena  (144). 

The  above  treatment  is  a  good  illustration  of  how  phase 
equilibria  considerations  and  information  are  used  in 
explaining  observed  phenomena. 


7 . 7  Summary  of  Surface  Phenomena  Studies 

In  this  chapter,  the  unusual  behavior  of  Cr  seen  in  this 
work  and  in  previous  studies  was  examined  and  explained  with 
the  help  of  information  learned  from  the  phase  equilibria, 
solubility,  and  diffusion  sections  of  this  work.  In 
particular,  the  large  Cr  surface  build-up  after  both 
indiffusion  and  outdif fusion  is  to  a  large  part  due  to  the 
formation  of  Ga (4)  droplets  or  inclusions  which  have  a  high 
Cr  solubility.  The  Cr  depletion  is  due  to  the  surface 
accumulation  of  manganese.  It  was  also  shown  how  substantial 
levels  of  Cr  in  LPE  layers  are  a  result  of  the  presence  of 
Ga^j  inclusions  and  not  Cr  in  solid  solution.  In  addition, 
the  proposed  diffusion  mechanism  for  Cr  in  GaAs  explains  how 
the  Cr  can  get  to  the  surface  region  so  quickly,  i.e.,  via  fast 


interstitials . 


8.  SUMMARY  AND  CONCLUSIONS 

In  this  chapter  a  review  of  this  research  is  given, 
followed  by  the  major  conclusions  of  the  work.  In  addition, 
suggestions  for  further  work  in  this  area  are  presented. 

The  present  research  concerns  the  phase  equilibria  of  the 
Ga-As-Cr  ternary  system  and  its  application  in  solubility  and 
diffusivity  studies  of  Cr  in  GaAs.  Using  the  anneal,  quench, 
and  x-ray  diffration  method  and  differential  thermal  analysis, 
the  Ga-As-Cr  ternary  phase  diagram  was  determined  from  room 
temperature  to  1300  °C  for  As  a/o  <  50.  This  information  was 
used  to  perform  well-defined  solubility  and  diffusion 
experiments,  which  were  done  by  annealing  GaAs  and  an 
appropriate  Ga-As-Cr  composition  -  as  the  Cr  source  or  sink  - 
for  various  times  and  temperatures  in  an  evacuated  closed 
quartz  ampoule,  and  measuring  the  Cr  concentration  vs.  depth  in 
the  GaAs  using  SIMS.  The  solubility  of  Cr  in  GaAs  was 
determined  as  a  function  of  position  in  the  ternary  phase 
diagram  from  700  to  1000 °C,  and  the  results  were  checked  by  C-V 
measurements.  Using  these  values  and  the  solubility  of  Cr  in 
Ga(£),  as  determined  in  the  phase  diagram  studies,  the 
distribution  coefficient  of  Cr  in  GaAs  was  determined  for  this 
same  temperature  range.  In  addition,  the  solubility  of  Cr  in 
doped  (n-type)  GaAs  was  studied.  The  diffusion  of  Cr  in  GaAs 
was  studied  by  performing  three  types  of  diffusion  experiments 
-  indiffusion,  outdif fusion ,  and  junction  diffusion  -  again 
under  well-defined  conditions.  The  diffusion  behavior  was 
modelled  by  a  complex  diffusion  mechanism,  and  values  for  key 
parameters  pertaining  to  the  model  were  determined,  including 
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two  types  of  diffusion  coefficients  corresponding  to  two 
different  regimes  of  the  diffusion.  In  addition,  unusual 
behavior  of  Cr  in  the  surface  region  of  GaAs  was  studied, 
specifically,  the  large  build-up  of  Cr  at  the  surface  after  all 
types  of  diffusion  and  depletion  of  Cr  near  the  surface  after 
indiffusion,  and  the  causes  for  such  behavior  were  determined. 
The  results  of  these  studies  were  applied  to  some  of  the 
results  reported  by  others,  which  explained  much  of  the 
observations  of  others  concerning  Cr  in  GaAs.  Finally  it  was 
shown  how  Cr  can  be  used  as  a  probe  to  study  the  defect 
chemistry  of  GaAs,  especially  concerning  other  impurities  and 
vacancies. 

The  major  conclusions  are  summarized  below: 

(1)  CrAs  and  Cr3Ga,  and  not  GaAs,  are  the  dominant  phases  in 
the  Ga-As-Cr  ternary  system,  thus,  most  other  phases  are 
in  equilibrium  with  them  and  not  with  GaAs. 

(2)  GaAs  and  Cr  are  not  in  equilibrium  with  each  other  under 
any  condition;  therefore,  when  heated  together  they  will 
react  with  each  other  to  form  Cr-Ga  and  Cr-As  phases. 

(3)  A  large  tie  triangle  region,  containing  GaAS(S), 

CrAs  (S),  and  Ga^j,  is  a  suitable  source  or  sink  for 
solubility  and  diffusion  studies  for  Cr  in  GaAs;  elemental 
Cr  is  not. 

(4)  The  solubility  of  Cr  in  GaAs  for  the  tie  triangle  region 
mentioned  above  increases  from  1.2  x  10^6  to  1.3  x 

1017  <^-3  as  the  temperature  increases  from  700  to 
1000°C.  The  solubility  of  Cr  at  a  given  temperature 
decreases  considerably  as  the  overall  composition  of  the 
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system  moves  toward  the  Ga-rich  corner. 


(5)  The  distribution  coefficient  increases  from  9  x  10”^  to 
3  x  10"5  from  700  to  1000°C,  but  is  virtually 
independent  of  Cr  concentration  at  a  given  temperature. 

(6)  The  Cr  solubility  increases  with  n-type  doping  in  GaAs  and 
decreases  with  p-type  doping,  consistent  with  Cr  being  a 
deep  acceptor  in  GaAs. 

(7)  The  diffusion  of  Cr  in  GaAs  follows  the  substitutional/ 
interstitial  dissociative  diffusion  mechanism.  In  this, 
there  are  three  regimes:  the  vacancy  diffusion,  the 
interstitial  diffusion,  and  the  vacancy  regeneration 
controlled  cases.  The  first  first  case  dominates  at 
shorter  annealing  times  and  lower  temperatures,  the  second 
case  dominates  at  longer  annealing  times  and  higher 
temperatures  (and  presumably  for  high  dislocation 
densities),  and  the  third  case  describes  the  behavior 

in  between.  The  effective  diffusion  coefficient  for  the 
first  case  ranges  from  4.0  x  10-13  to  1.2  x  10“-*--'- 
cm^/sec.  for  700  to  900°C,  with  an  activation  energy  of 
1.7  eV,  and  for  the  second  case,  ranges  from  1.5  x  lO-^ 
to  9.1  x  10~10  cm^/sec  from  750  to  1000°C,  with  an 
activation  energy  of  0.8  eV.  The  equilibrium  time 
constant  (or  the  vacancy  regeneration  rate),  which 
describes  the  transition  from  the  first  case  to  the 
second,  decreases  from  1.3  x  10^  to  3.0  x  10^  sec. 
from  750  to  900  °C,  for  LEC  GaAs.  It  is  considerably  less 
at  each  temperature  for  Bridgman  GaAs. 

(8)  The  Cr  diffusion  and  solubility  is  affected  by  defects  in 
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GaAs,  especially  other  imparities  and  vacancies.  For 
example,  donor  impurities  increase  the  Cr  solubility  but 
decrease  the  diffusivity  -  acceptors  do  the  opposite,  and 
increasing  the  gallium  vacancy  concentration  decreases  the 
time  it  takes  for  the  extrapolated  surface  concentration 
of  Cr  to  reach  the  solubility  level.  In  addition,  such 
defects  as  dislocations  introduced  by  polishing,  etc., 
probably  affect  Cr  behavior. 

(9)  The  anomalous  build-up  of  Cr  at  the  surface  of  GaAs  is 
mostly  due  to  the  presence  of  a  second  phase  -  Ga(Zj  - 
which  forms  upon  the  incongruent  vaporization  of  arsenic 
and  which  has  a  very  high  solubility  of  Cr  ( =10  a/o  at 

700 °C).  Similarly,  the  large  Cr  concentration  seen  in  LPE 
layers  after  growth  using  the  meltback  technique  is  due  to 
the  incorporation  of  Cr  in  Ga(£)  inclusions  that  remain 
in  the  layer. 

(10)  The  depletion  region  of  Cr  observed  near  the  surface  after 
indiffusion  is  due  to  uneven  Mn  distribution,  which  lowers 
the  solubility  of  Cr  in  that  region. 

(11)  The  gettering  of  Cr  to  the  surface  is  accomplished  via 
fast  Cr  interstials. 

The  following  are  suggestions  for  further  work  in  this 

area: 

(1)  Continue  the  use  of  Cr  as  a  probe  to  study  the  defect 
chemistry  of  GaAs.  As  discussed  earlier,  since  the 
behavior  of  Cr  is  very  sensitive  to  the  presence  and 
behavior  of  defects  in  GaAs,  the  defects  can  be  studied  by 
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studying  Cr.  This  would  be  particularly  useful  in 
investigating  vacancies  and  dislocations,  which  are 
inherently  difficult  to  study  directly.  Therefore,  more 
work  could  be  done  studying  the  effects  of  dislocations 
and  vacancies  on  the  behavior  of  Cr  in  GaAs. 

(2)  The  same  types  of  studies  done  in  this  research  on  the 
Ga-rich  side  of  the  phase  diagram  can  be  done  on  the 
arsenic-rich  side.  Since  many  GaAs  device  manufacturers 
are  going  to  the  use  of  arsenic  overpressure  rather  than 
encapsulents  to  prevent  surface  degradation,  this  would 
have  practical  applications. 

(3)  To  study  Cr  diffusion  in  the  absence  of  a  concentration  or 
activity  gradient,  one  can  do  isoconcentration  diffusion 
experiments.  This  also  eliminates  any  non-uniform 
electrical  effects  on  non-equilibrium  effects.  This  can 
easily  be  done  without  the  use  of  radioisotopes  by 
diffusing  Cr^3 ,  a  non-radioactive  isotope,  in 
Cr^^-doped  GaAs,  and  using  SIMS  to  measure  the  Cr^3 
concentrations  (SIMS  is  mass  specific).  However,  since: 

1)  the  concentration  of  Cr  in  GaAs  is  so  small  (ppb),  the 
activity  coefficient  is  probably  constant  with 
concentration  and  virtually  equal  to  unity  anyway;  2)  the 
diffusion  in  most  cases  follows  the  uncharged  model?  and 
3)  equilibrium  is  achieved  rather  quickly,  especially  at 
high  temperatures,  not  much  difference  is  expected  for  the 
behavior  of  Cr  in  isoconcentration  conditions  compared  to 
concentration  gradient  conditions,  except  at  low 


temperatures  (where  the  concentrations  of  Cr  and 
electrically  active  defects  are  less  than  nj_ )  and  for 
crystals  with  a  large  number  of  defects  present. 

Study  the  interaction  between  C_-  and  the  other  materials 
present  in  GaAs  devices,  especially  the  materials  used  to 
make  electrical  contacts  to  GaAs  These  include  Au,  Ge, 
Sn,  Se,  and  Zn.  Since  Cr  has  a  tendency  to  end  up  in  the 
surface  region  of  GaAs  and  to  react  with  other  materials, 
this  could  be  a  major  problem  in  GaAs  devices.  Not  only 
must  the  binary  phase  equilibria  between  Cr  and  the  other 
elements  be  considered,  but  also  the  multi-component 
phase  equilibria  between  all  the  elements  present  must  be 
considered.  The  same  concepts  used  in  the  present  study 
can  be,  and  should  be,  applied  to  these  cases. 

Now  that  the  mode  of  Cr  diffusion  in  GaAs  and  the  phase 
equilibria  are  known,  one  might  now  attempt  to  determine  a 
way  to  keep  Cr  from  adversely  affecting  the  performance  of 
GaAs  devices.  This  may  include  preannealing  the  Cr-doped 
substrate  under  conditions  such  that  it  is  saturated  with 
gallium  vacancies,  thus  reducing  the  tendency  for 
conversion  of  Cr  to  the  interstitial  state  and  thus 
reducing  the  fast  interstitial  diffusion.  Another 
possibility  is  using  backside  gettering  or  internal 
gettering,  such  as  are  used  in  silicon  I.C.  processing,  to 
keep  the  Cr  from  getting  into  the  active  region. 


APPENDIX  -  ERROR  ANALYSIS 
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Error  bars  indicating  the  maximum  expected  error  for  the 
particular  determinations  have  been  included  with  each  data 
point  in  figures  4-7,  4-9,  6-8,  and  6-9.  The  reasoning  behind 
such  determinations  for  each  case  is  as  follows: 

Cr  Solubility,  Figure  4-7 

As  discussed  in  Section  4.4.1,  the  reported  maximum 
expected  error  for  the  SIMS  analyses  of  Cr  in  GaAs  is  a  factor 
of  1.5  to  2 ,  and  the  standard  deviation  of  several  sample 
results  is  +  50%.  The  solubility  values  were  taken  directly 
from  SIMS  Cr  profiles  after  appropriate  anneals,  and  therefore 
the  error  in  the  solubility  is  just  the  error  in  SIMS  analysis 
(the  error  in  extrapolating  the  profile  to  the  surface  was 
trivial  in  comparison  to  the  error  in  SIMS),  Thus,  the  maximum 
expected  error  for  the  solubility  values  is  determined  to  be  "a 
factor  of  1.5".  This  means  that  for  a  reported  value  of  2  x 
lO-^,  for  example,  the  true  value  could  be  from  1.3  x  lO-*-^ 
to  3  x  10^-6.  This  type  of  expected  error,  rather  than  +  50 
percent,  is  much  more  suited  to  data  that  is  presented  on  a  log 
scale. 

Distribution  Coefficient,  Figure  4-9 

The  distribution  coefficient,  k,  was  calculated  by 
dividing  the  solid  solubility  of  Cr  in  GaAs  by  the  liquid 
solubility  of  Cr  in  the  Ga-rich  liquid.  The  error  in  the  solid 
solubility  is  determined  to  be  "a  factor  of  1.5",  as  discussed 
above.  The  error  in  the  liquid  solubility,  from  the  phase 
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diagram  determinations,  is  estimated  to  be  a  factor  of  1.3, 
resulting  in  an  overall  total  expected  error  of  a  factor  of  2 
Diffusion  Coefficients,  Figure  6-8 


The  diffusion  coefficients  were  determined  by  fitting  the 
diffusion  model  profiles  to  the  actual  profiles.  The  value  of 
the  diffusion  coefficient  depends  on  three  parameters:  1)  the 
relative  concentration;  2)  the  depth;  and  3)  the  annealing 
time.  The  relative  concentration  refers  to  the  ratio  of  the 
impurity  concentration  at  a  given  depth  to  either  the  surface 
concentration  (for  indiffusion  results)  or  the  bulk 
concentration  (for  outdif fusion  results).  In  SIMS  analysis, 
while  the  error  in  the  absolute  concentrations  are  somewhat 
large,  the  error  in  relative  concentration  are  regarded  as  very 
small.  This  is  especially  true  when  the  matrix  atom  profiles 
are  very  uniform,  indicating  uniform  sputtering,  as  is  the 
case  in  this  research.  Thus,  this  error  is  considered 
negligible.  The  depth  determination,  as  measured  by  a  DEKTAC 
instrument,  has  a  reported  error  of  15%.  Since  D  is 
proportional  to  the  square  of  the  depth,  this  results  in  an 
error  of  1  -  (1.15  x  1.15),  or  32%.  Finally,  the  error  in 
annealing  time,  which  also  takes  into  account  that  it  takes 
some  time  for  the  sample  to  reach  the  desired  temperature,  is 
estimated  to  be  10%.  D  is  inversely  proportional  to  t,  so  the 
net  error  is  determined  to  be  1— (1.32  x  1.1)  -  50%.  In  order 
to  place  this  on  a  log  scale,  error  bars  representing  a  factor 


of  1.5  error  are  used. 


The  equilibrium  time  constant,  0,  for  each  temperature  was 
determined  from  the  ratio  of  the  extrapolated  surface 
concentration  to  the  solubility  concentration  at  that 
temperature,  and  the  annealing  time  through  the  equation: 

0  =  -t/ln ( 1-Cq/C^ ) .  The  best  way  to  determine  the  maximum 
error  for  this  is  by  taking  the  standard  deviation  of  the 
results  at  each  temperature  for  the  various  annealing  times 
(the  value  for  0  should  be  the  same  for  all  times  at  each 
temperature).  The  maximum  expected  error  is  thus  determined  to 
be  a  factor  of  1.5. 

Error  in  Temperature 

In  all  of  the  above  mentioned  graphs,  the  values  are 
plotted  against  temperature.  The  absolute  error  in  the 
temperature  is  estimated  to  be  +  5°C,  or  »0.5%,  while  the 
relative  error,  from  one  experiment  to  the  next,  is  only 
+  1°C,  or  *  0.1%.  The  vertical  error  bars  for  this  were  not 
included  in  the  graphs;  however,  one  should  be  aware  of  this 


error  in  addition  to  the  others. 
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